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20 Background of the Invention 

Field of the Invention 

The present invention is related to digital signal processing and digital 
conmiunications. More specifically the present invention is related to interpolation, 
25 angle rotation, rectangular-to-polar conversion, and carrier and symbol timing 

recovery for digital processing and digital communications applications. 
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RelatedArt 



Advances in technology have enabled high-quality, low-cost 
communications with global coverage, and provide the possibility for fast Internet 
access and multimedia to be added to existing services. Exemplary emerging 
5 technologies include cellular mobile radio and digital video broadcasting, both of 

which are described briefly as follows. 

In recent years, cellular mobile radio has experienced rapid grov^h due to 
the desire for mobility while enjoying the two-way voice services it provides. 
GSM, IS-136 and personal digital cellular (PDC) are among the most successfiil 

10 second-generation personal communications (PCS) technologies in the world 

today, and are responsible for providing cellular and PCS services globally. As the 
technology advances, customers will certainly demand more from their wireless 
services. For example, with the explosive growth of the world wide web over the 
wired networks, it is desirable to provide Internet services over mobile radio 

1 5 networks. One effort to specify the future global wireless access system is known 

as IMT-2000 (Buchanan, K., etal., lEEEPers, Comm. 4:S-13 (1997)). The goal 
of IMT-2000 is to provide not only traditional mobile voice communications, but 
also a variety of voice and data services Avith a wide range of applications such as 
multimedia capabilities, Internet access, imaging and video conferencing. It is also 

20 an aim to unify many existing diverse systems (paging, cordless, cellular, mobile 

satellite, etc.) into a seamless radio structure offering a wide range of services. 
Another principle is to integrate mobile and fixed networks in order to provide 
fixed network services over the wireless infrastructure. Such systems might well 
utilize broadband transport technologies such as asynchronous transfer mode 

25 (ATM). 

For the applications of IMT-2000, a high-bit-rate service is needed. 
Moreover, for multimedia applications, the system should provide a multitude of 
services each requiring 1) a different rate, and 2) a different quality-of-service 
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parameter. Thus, a flexible, variable-rate access with data rates approaching 2Mb/s 
is proposed for IMT-2000. 

The advent of digital television systems has transformed the classical TV 
channel into a fast and reliable data transmission medium. According to the 
specifications of the DVB project (Reimers, U., IEEE Comm. Magazine 36: 104- 
110 (1998)), digital TV is no longer restricted to transmitting sound and images 
but instead has become a data broadcasting mechanism which is fiilly transpsirent 
to all contents. Digital TV broadcasting by satellite, cable and terrestrial networks 
is currently under intensive development. A typical system looks like this: a DVB 
signal is received fi-om a satellite dish, fi-om cable, or fi-om an antenna (terrestrial 
reception). A modem buih into an integrated receiver/decoder (IRD) will 
demodulate and decode the signal. The information received will be displayed on 
a digital TV or a multimedia PC. In addition to being used as a digital TV, DVB 
can receive data streams fi-om companies who wish to transmit large amounts of 
data to many reception sites. These organizations may be banks, chains of retail 
stores, or information brokers who wish to offer access to selected Internet sites 
at high data rates. One such system is MultiMedia Mobile (M% which has a data 
rate of 16 Mb/s. 



S^OT proper operation, these third generation systems require proper 
synchronizatton between the transmitter and the receiver. More specifically, the 
frequency and pha^esof the receiver local oscillator should substantially match that 
of the transmitter local osc^lator. When there is a mismatch, then an undesirable 
rotation of the symbol constelfet^on will occur at the receiver, which will seriously 
degrade system performance. Whm^he carrier frequency offset is much smaller 
than the symbol rate, the phase and fre^u^cy mismatches can be corrected at 
baseband by using a phase rotator. It is also necfe^ary to synchronize the sampling 
clock such that it extracts symbols at the correci^limes. This can be achieved 
digitally by preforming appropriate digital resamples. 

The digital resampler and the direct digital frequency synthesizer (DDS) 
used by the phase rotator are among the most complex components in a receiver 




1904.0140002 



(Cho, K., "A frequency-agile single-chip QAM modulator with beamforming 
diversity," Ph.D. dissertation. University of California, Los Angeles (1999)). Their 
performance is significant in the overedl design of a communications modem. For 
multimedia communications, the high-data-rate requirement would impose a 
demand for high computational power. However, for mobile personal 
communication systems, low cost, small size and long battery life are desirable. 
Therefore, it would be desirable to have an efficient implementation of the phase 
rotator, re-sampler, and DDS in order to perform fast signal processing that 
operates within the available resources. Furthermore, it would be desirable to have 
an efficient synchronization mechanism that uses a unified approach to timing and 
carrier phase corrections. 

For Internet services it is important to provide instantaneous throughput 
intermittently. Packet data systems allow the multiplexing of a number of users on 
a single channel, providing access to users only when they need it. This way the 
service can be made more cost-effective. However, the user data content of such 
a transmission is usually very short. Therefore, it is essential to acquire the 
synchronization parameters rapidly from the observation of a short signal-segment. 

For applications where low power and low complexity are the major 
requirements, such as in personal communications, it is desirable to sample the 
signal at the lowest possible rate, and to have a synchronizer that is as simple as 
possible. Therefore, it is also desirable to have an efficient synchronizer 
architecture that achieves these goals. 

For applications utilizing Orthogonal Frequency Division Multiplexing 
(OFDM), sampling phase shift error produces a rotation of the Fast Fourier 
Transform (FFT) outputs (PoUet T,, and Peters, M., IEEE Comm. Magazine 
37:80-86 (1999)). A phase correction can be achieved at the receiver by rotating 
the FFT outputs. Therefore, it is also desirable to have an eflBcient implementation 
structure to perform rotations of complex numbers. 



1904.0140002 



-5- 



Summary of the Invention 



The present invention is directed at an apparatus and method for 
preforming angle rotation. More specifically, the angle rotator rotates an input 
complex signal in the complex plane according to the angle 6, to produce a 
5 rotated output signal. In doing so, the angle rotator preforms a two stage rotation 

that includes a coarse stage rotation and a fine stage rotation. The two specific 
amounts of rotation are obtained directly fi^om the original angle, without 
performing iterations as performed by known CORDIC-type methods. Therefore, 
the critical path of the present invention is significantly shorter than that of the 

10 CORDIC-type methods. 

The coarse stage rotation is performed using truncated approximations for 
the cosine and the sine 6^. where 6^ is a radian angle that corresponds to a 
most significant word (MSW) of the input angle 6. The resulting intermediate 
complex signal is rotated by a second amount during the fine stage rotation. The 

1 5 fine stage rotation is performed using one or more error values that compensate 

for approximations and quantization errors associated with the coarse stage 
rotation. 

By partitioning the rotation into coarse and fine rotation stages, a two 
stage structure is obtained that requires much less hardware than a single stage 

20 rotator, without sacrificing angle precision. This can occur because the two stage 

rotator stores pre-computed cosine and the sine 6^ values in a small lookup 
table (e.g. memory device) for fast retrieval. This same lookup table can be reused 
by a phase rotator for carrier recovery, in systems wherein both subsystems are 
present. Furthermore, the angle rotator consolidates all operations into a small 

25 number of reduced-size multipliers, enabling the use of efficient multiplier 

implementations, such as Booth encoding, thereby yielding a smaller and faster 
overall circuit. When higher precision is desired, more accurate results can be 
attained simply by increasing the wordlength and the multiplier size, v^thout 
increasing overall circuit latency. Further efficiency in implementation can be 
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achieved by exploiting certain opportunities presented by Booth encoding methods 
in relation to the angle rotator's special "two-butterfly" structure. 

Although the angle-rotation processor is initially intended for the 
trigonometric-based resampler and the phase rotator for carrier recovery, its 
application well exceeds that of synchronization. It can be applied to efficiently 
implement many communication functions and non-communications functions, 
wherever the operations can be characterized as multiplying complex numbers by 
a complex exponential. Such applications include but are not limited to Quadrature 
Direct Digital Frequency Synthesizer/Mixer (QDDFSM) applications. Direct 
Digital Frequency Synthesizer (DDFS) applications. Fast Fourier Transform (FFT) 
and Inverse FFT applications, and, more generally, Discrete Fourier Transform 
applications. 

Further features and advantages of the invention, as well as the structure 
and operation of various embodiments of the invention, are described in detail 
below with reference to the accompanying drawings. The drawing in which an 
element first appears is typically indicated by the leftmost character(s) and/or 
digit(s) in the corresponding reference number. 

Brief Description of the Figures 

FIG. 1 A illustrates a PSK transmitter. 
FIG. IB illustrates a PSK receiver. 
FIG. IC illustrates a block diagram of an OFDM system. 
FIG. ID illustrates a PSK receiver with carrier and timing recovery. 
FIG. 2 illustrates an Interpolation environment. 
^^IGr^-ittngtrates a l^agrange baiiiii pofyneHuals, — 
FIG. 4 illustrates a Farrow structure that implements (2.5) and (2.6). 
FIG. 5 illustrates a flowchart 500 representing trigonometric interpolation 
according to embodiments of the present invention. 

FIG. 6A. illustrates an impulse response of a Lagrange interpolator. 
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FIG. 6B illustrates an impulse response of a Trigonometric interpolator 
according to embodiments of the present invention. 

FIG.7A illustrates a frequency response for N-4 according to 
embodiments of the present invention. 

FIG.7B illustrates a frequency response for A^=32 according to 
embodiments of the present invention. 

FIG. 8 A illustrates a signal with two samples/symbol and 1 00% excess BW 
according to embodiments of the present invention. 

FIG. 8B illustrates an NMSE of the interpolated signal according to 
embodiments of the present invention. 

FIG. 9 illustrates the critical path of the Lagrange cubic interpolator. 

FIG. 10 illustrates a trigonometric interpolator with N=^ according to 
embodiments of the present invention. 

FIG. 1 1 illustrates a trigonometric interpolator with i\^=8 according to 
embodiments of the present invention. 

FIG. 12 illustrates a conceptual modification of input samples according 
to embodiments of the present invention. 

FIG. 13 illustrates correcting the offset due to modification of original 
samples according to embodiments of the present invention. 

FIG. 14 illustrates the modified trigonometric interpolator for A^=4 
according to embodiments of the present invention. 

FIG. 15 illustrates the modified trigonometric interpolator for A^=8 
according to embodiments of the present invention. 

FiGJ^A-D illustrates a comparision of the amount of interpolation error 
using (A) Language^^Cnbio.,,^) the trigonometric interpolator 1000, (C) the 
trigonometric interpolator, (D) tRB^-'-epjimal structure (to be discussed in 
Section 4). 

FIG. 17 illustrates a flowchart 1700 representing trigonometric 
interpolation according to embodiments of the present invention. 
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FIG. 1 8 illustrates trigonometric interpolation using a table lookup for 
angle rotation according to embodiments of the present invention, 

FIG. 19 illustrates trigonometric interpolation using modified samples 
according to embodiments of the present invention. 
5 FIG. 20 illustrates normalized impulse responses / of the interpolation 

filters according to embodiments of the present invention. 

FIG. 21 illustrates normalized fi-equency responses F of the interpolation 
filters according to embodiments of the present invention. 

FIG. 22 illustrates analysis of the fi-equency responses according to 
10 embodiments of the present invention. 

FIG. 23 illustrates the effect of a more gradual transition at the band edge 
according to embodiments of the present invention. 

FIG. 24 illustrates reducing the transition bandwidth by increasing 
according to embodiments of the present invention. 
1 5 FIGs. 25 A-B illustrate (A) impulse response of the original filter and the 

modified filter; (B) The equivalent v^ndow, according to embodiments of the 
present invention. 

FIG. 26 illustrates forming the fi-equency response of the discrete-time 
fi-actional-delay filter according to embodiments of the present invention. 
20 FIGs. 27 A-B illustrate a fractional-delay filter with (A) ^=0.12 and 

(B) n=0.5, using the preliminary A^=8 interpolator according to embodiments of 
the present invention. 

FIGs. 28 A-D illustrate modification to F(Q), and the corresponding F^(o>) 
according to embodiments of the present invention. 
25 FIGs. 29 A-B illustrate F/g)), with n=0.5, N=S, (A) before and (B) after 

optimization according to embodiments of the present invention. 

FIGs. 30A-30B illustrate F^(o>) for ji=0.5, A^=4, (A) before and (B) after 
modification according to embodiments of the present invention. 
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FIGs. 31A-31B illustrate F^(o)), ^==0.5, simplified A^=4 structure, 
(A) before and (B) after modification according to embodiments of the present 
invention. ^ 

FIG. 32 illustrates real and imaginary components of the ^ value 

according to embodiments of the present invention. 

ates signal with two samples/symbol and 40% excess 
bandwidth embodiments of thepresent invention. 

FIG. 34 illustrates a flowchart 3400 for optimizing trigonometric 
interpolation according to embodiments of the present invention. 
10 FIG. 35 illustrates a flowchart 3500 for optimizing trigonometric 

interpolation according to embodiments of the present invention. 
5 FIG. 36 illustrates an optimized interpolator 3600 according to 

embodiments of the present invention. 
Cn FIG. 37 illustrates an optimized interpolator 3700 according to 

U 15 embodiments of the present invention. 

FIG. 3 8 illustrates an angle rotator 3 800 according to embodiments of the 
□ present invention. 

5 FIG. 39 illustrates an angle rotator 3900 according to embodiments of the 

present invention. 

20 FIG. 40 illustrates an angle rotator 3900 and example multiplier sizes 

according to embodiments of the present invention. 

FIG. 41 illustrates a flowchart 4100 for angle rotation according to 
embodiments of the present invention. 

FIG. 42 illustrates an angle rotator 3900 and multiplier sizes to achieve 
25 90.36 dB SFDR according to embodiments of the present invention. 

FIG. 43 illustrates an output spectrum showing 90.36 dB SFDR according 
to embodiments of the present invention. 

FIG. 44 illustrates a modified angle rotator 4400 when only one output is 
needed according to embodiments of the present invention. 
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FIG. 45 illustrates a flowchart 4500 for angle rotation when only one 
output is needed according to embodiments of the present invention. 
FIG. 46 illustrates a phase accumulator 4600. 

FIG. 47 illustrates a Quadrature Direct Digital Frequency 
Synthesizer/Mixer (QDDFSM) 4700 according to embodiments of the present 
invention. 

FIG. 48 illustrates an angle rotator 4800 according to embodiments of the 
present invention, 

v^*^^^' FJ<T 49 illustrates a Booth multiplier according to embodiments of the 
present invention. 

FIG. 50 illustf^s an original Booth table 5000. 

FIG, 5 1 illustrates annfcgating booth table 5 1 00 according to embodiments 
of the present invention. 

FIG. 52 illustrates an negating Boothsqiultiplier 5200. 

FIG. 53 illustrates a conditionally negatingBq^h decoder 5300. 

FIG. 54 illustrates a conditionally negating multipri'^400. 

FIG. 55 illustrates as an angle rotator configured as quadrature direct 
digital synthesizer (QDDS) 5500 according to embodiments of the present 
invention. 

FIG. 56 illustrates an angle rotator as a cosine only direct digital 
synthesizer based on angle rotator 3900 according to embodiments of the present 
invention. 

FIG, 57 illustrates an angle rotator as a cosine only direct digital 
synthesizer based on angle rotator 4400 according to embodiments of the present 
invention. 

FIG. 58 illustrates a common packet format for packet based 
communications. 

FIG. 59 illustrates a system model for packet based communications. 
FIG. 60 illustrates mean values of a preamble correlator output, for 0 = 0, 
according to embodiments of the present invention. 
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FIG. 61 illustrates a synchronizer 6100 according to embodiments of the 
present invention. 

FIG. 62 illustrates a flowchart 6200 associated with the synchronizer 6100 
according to embodiments of the present invention. 
5 FIG. 63 illustrates bias due to truncation. 

^^^^^FIGr-64illustrates a synchronizer 6400 according to embodiments of the 
present invention. 

FIG. 65A-65B illustrafea^owchart 6200 associated with the synchronizer 
6100 according to embodiments of m&wesent invention. 
□ 10 FIG. 66 illustrates timing variance, a= 0.4. 

fS FIG. 67 illustrates timing jitter variance, a = 0. 1 . 

FIG. 68 illustrates phase jitter variance, a = 0. 1 . 
rU FIG. 69 illustrates Cartesian-to-polar conversion. 

S FIGs. 70A-70B illustrate using Newton-Raphson iteration to find 1/^, 

15 FIG. 71 illustrates a rectangular-to-polar converter 7100 according to 

B embodiments of the present invention. 

Q FIG. 72 illustrates angle rotation associated v^th the rectangular-to-polar 

;S converter 7100 according to embodiments of the present invention. 

FIG. 73 illustrates a flowchart 7300 associated with the synchronizer 7 1 00 
20 according to embodiments of the present invention. 

FIG. 74 illustrates Interpolation in a non-center interval according to 
embodiments of the present invention. 
VT^^ ^tes774A;^Ulustrate impulse responses of the non-center-interval 
interpolation filter (^"^Ire^Qve and (b) after optimization, according to 
25 embodiments of the present inventiSt 

FIGs. 76A-B illustrate firequency responses of the non-center-interval 
interpolator (A) before optimization and (B) after optimization, according to 
embodiments of the present invention. 

FIG. 77 illustrates an exemplary computer system 7702, according to 
30 embodiments of the present invention. 
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FIG. 78 illustrates a data rate expansion circuit 7800 according to 
embodiments of the present invention. 

Detailed Description of the Preferred Embodiments 
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L Introduction 

As discussed herein, third generation and other cutting edge 
communications systems require proper synchronization between the transmitter 
and the receiver. More specifically, the fi-equency and phase of the receiver local 
oscillator should substantially match that of the transmitter local oscillator and 
accurate symbol timing must be achieved. The following section discuss some 
exemplary modulation schemes and configurations, and their related 
synchronization issues. These example configurations are not meant to be limiting, 
and are provided for example purposes only. After which, an overview of the 
present invention is provided. 

i. 1 Exemplary Modulation Schemes and Synchronization Issues 

A key to the evolution of third-generation PCS is the ability to provide 
higher data rates via increased spectral efficiency of the access scheme. The IS- 
136 conmiunity intends to add a 200-KHz carrier bandwidth and adopt 8PSK 
modulation. This allows for data rates up to 384 Kb/s. 

A simplified 8PSK transmitter 102 and receiver 104 are shovra in FIG. 1 A 
and FIG. IB, respectively. The receiver 104, as shovm, performs baseband 
sampling. Alternatively, the received signal could be sampled at an EF fi"equency, 
where the down-conversion to baseband is performed digitally. However, since 
it does not alter the main subject in the present invention, the baseband-sampled 
system is used as an example. 

Referring to FIG. IB, PSK receiver 104 down-converts an IF input signal 
106 to baseband by multiplication with a local oscillator signal 108 using mixers 
110. After filtering 111, A/D converters 112 sample the dovm-converted signal 
according to a sampling clock 114 in preparation for logic examination. After 
fiirther filtering 116 and equalization 118, the logic decision devices 120 examine 
the sampled signal to determine a logic output for the two channels. 

1904.0140002 



-16- 



During down-conversion, an undesirable rotation of the symbol 
constellation will occur if the frequency and phase of the oscillator signal 108 does 
not match the oscillator signal of the transmitter 102. This symbol rotation can 
seriously degrade system performance. When the carrier frequency offset is much 
smaller than the symbol rate, the phase and frequency mismatches can be corrected 
at baseband, using a phase rotator 124, as shown in FIG. ID. 

The sampling clock 1 14 is generated locally in the receiver 104. The logic 
decision devices 120 make more accurate decisions when the sampling instant is 
optimal, i.e., sjmchronous to the incoming symbols. 

If the timing information can be extracted from the signal 106, it can be 
used to adjust the phase of the sampling clock 1 14. This adjustment would require 
a voltage controlled oscillator (VCO) to drive the A/D converters 112 . In this 
scenario, the digital portion of the circuit 1 04 needs to keep in synchronization with 
the A/D converters 112, which places strict requirements on the VCO. Moreover, 
changing the phase of the sampling clock 1 14 would cause jitter. High data-rate 
receivers are more sensitive to such jitter when used in multimedia 
communications. 

Another solution to timing errors is to correct them entirely in the digital 
domain, with the equivalent of A/D sampling adjustment performed by a digital 
resampler 122, as shown in FIG. ID. This resampler 122 is controlled by a timing 
recovery circuit (not shown) and it attempts to supply the optimal samples (i.e. 
synchronous) to the decision circuits 120. Using the digital resampler 122, the 
timing recovery loop is closed entirely in the digital domain. This allows the 
complete separation of digital components from analog components. 

The digital resampler 122 and a direct digital frequency synthesizer (not 
shown) used by the phase rotator 124 are among the most complex components in 
a receiver (Cho, K., "A frequency-agile single-chip QAM modulator with 
beamforming diversity," Ph.D. dissertation, University of California, Los Angeles 
(1999)). Their performance is significant in the overall design of the modem. For 
multimedia communications, the high-data-rate requirement imposes a demand for 
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high computational power. However, for mobile personal communication systems, 
low cost, small size, and long battery life are desirable. Therefore, efficient 
implementation is the key to implementing fast signal processing within the 
available resources. It is also desirable to provide an eflRcient synchronization 
mechanism by using a unified approach to timing and carrier phase corrections. 
This can be accomplished by sharing resources between the resampler 122 and the 
phase rotator 124. 

As for the digital video broadcasting system (DVB) systems, the most 
challenging of all DVB transmissions is the one used in terrestrial channels (DVB- 
T) due to the presence of strong echoes which characterize the propagation 
medium. A common approach for DVB-T is based on Coded-OFDM (orthogonal 
frequency division multiplexing). The major benefit of OFDM is that the serial 
baseband bitstream which needs to be transmitted is distributed over many 
individual subcarriers. Such spreading makes the signal robust against the effects 
of multipath and narrowband interference. The simplified block diagram of an 
OFDM modem 108 is shown in FIG. IC. 

FIG. IC illustrates an orthogonal frequency division multiplexing system 
(OFDM) 126 having an OFDM transmitter 128 and an OFDM receiver 130. For 
the OFDM system 126, synchronization errors produce a rotation of the fast 
Fourier Transform (FFT) outputs of the OFDM receiver 130. (PoUet T., and 
Peeters, M., IEEE Comm. Magazine 37:80-86 (1999)). A sampling phase 
correction for the received signals can be achieved by rotating the FFT outputs at 
the receiver. For FFT rotation, it is desirable to have an efficient implementation 
structure to perform rotations of complex numbers. 

The example applications and modulation schemes described above in this 
section were provided for illustrative purposes only, and are not meant to be 
limiting. Other applications and combinations of such applications will be apparent 
to persons skilled in the relevant art(s) based on the teachings contained herein. 
These other applications are within the scope and spirit of the present invention. 
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1.2 Overview 



LU 



The following is an overview of the sections that follow. 
Sections 2, 3 and 4, discussed herein, present a novel interpolation method 
for digital resampling using a trigonometric polynomial. In Section 2, after a brief 
5 review of interpolation methods, particularly those using a conventional 

polynomial, the trigonometric interpolation method is introduced. Efficient 
implementation structures for trigometric interpolation are given. The performance, 
hardware complexity and computational delay are compared with conventional 
J polynomial interpolators. The trigonometric- polynomial based resampler can use 

p 1 0 the same hardware as is employed in the phase rotator for carrier synchronization, 

thus further reducing the total complexity in the synchronization circuitry. 

In Section 3, a signal processing approach is used to analyze the 
interpolation method devised in Section 2. It shows how an arbitrary frequency 
response is achieved by applying a simple modification to the original interpolation 
1 5 algorithm. This enables the interpolator to also perform matched filtering of the 

received signal. 

The approaches in Section 3 can be employed to improve the interpolator 
performance by optimizing the frequency response of the continuous-time 
interpolation filter. This method is based on optimizing the performance by 
20 conceptually reconstructing the continuous-time signal from existing samples. 

From the point of view of designing digital resamplers, however, what we are 
actually interested in are new samples corresponding to the new sampling instants. 
In Section 4, we optimize the interpolation filter such that the error in producing 
a new sample corresponding to every resampling instant is minimized, hence further 
25 improving the overall interpolation accuracy. 

Section 5 presents an angle rotation processor that can be used to efficiently 
implement the trigonometric resampler and the carrier phase rotator. This structure 
can also implement the resampler for an OFDM receiver, which rotates the FFT 
outputs. It has many other practical applications. 
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The discussions in the previous Sections have assumed that the sampling 
mismatch that is supplied to the resampier is known. The problem of obtaining the 
synchronization parameters is studied in Section 6. For burst mode transmissions 
in packet data systems, we present an eflBcient architecture for feedforward symbol- 
5 timing and carrier-phase estimation. 

Section 7 presents an efficient implementation of a key component in the 
feedforward synchronizer of Section 6, as well as in many other such synchronizers. 
This involves computing the angle from the real and imaginary components of a 
complex number. The discussion, however, extends to a general problem of 
10 Cartesian-to-polar conversion, which is encountered in many communication 

applications. An architecture that efficiently accomplishes this conversion is 
presented. 

Section 8 presents an exemplary computer system in which the invention 
can be operated. 
1 5 Section 9 includes various appendices. 

Further discussions related to materials in Sections 2-7 are included in 
Dengwei Fu, "Efficient Synchronization for Multimedia Communications," Ph.D 
dissertation, University of California, Los Angeles, 2000, which is incorporated-by- 
reference, in its entirety. 
20 Additionally, the following articles are herein incorporated by reference: 

D. Fu and A. N. Willson, Jr., "A high-speed processor for digital sine/cosine 
generation and angle rotation," in Conf. Record 32nd Annual Asilomar Conference 
on Signals, Systems and Computers^ vol. 1, pp, 177-181, Nov. 1998; 

D. Fu and A. N. Willson, Jr., "Interpolation in timing recovery using a 
25 trigonometric polynomial and its implementation," in Proc. GLOBECOM 1998, 

Comm. Theory Mini-Conference Record, pp. 173-178, Nov. 1998; 

D. Fu and A. N. Willson, Jr., "Design of an improved interpolation filter 
using a trigonometric polynomial," \nProc, Int Symp. Circuits & Systems, vol. 4, 
pp. 363-366, May 30- June 3, 1999; 
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D. Fu and A. N. Willson, Jr., "A high-speed processor for rectangular-to- 
polar conversion with applications in digital communications," in Proc, 
GLOBECOM 1999, vol. 4, pp. 2172-2176, Dec. 1999; 

D. Fu and A. N. Willson, Jr., "A fast synchronizer for burst modems with 
simultaneous symbol timing and carrier phase estimations," in Proa Int Symp, 
Circuits & Systems, vol. 3, pp. 379-382, May 28-31, 2000; and 

D. Fu and A. N. Willson, Jr., "Optimal interpolator using a trigonometric 
polynomial," inProc, 4 3rd Midwest Symp, Circuits & Systems, 4 pages, Aug. 8-11, 
2000. 

2. Interpolation Using a Trigonometric Polynomial 

As discussed in Section 1, when an analog-to-digital converter (ADC) is 
clocked at a fixed rate, the resampler must provide the receiver with correct 
samples, as if the sampling is synchronized to the incoming symbols. How can the 
resampler recover the synchronized samples by digital means without altering the 
sampling clock? Since the input analog signal to the ADC is bandlimited, as long 
as the sampling rate is at least twice the signal bandwidth, according to the 
sampling theorem, the sampled signal carries as much information as the 
continuous-time signal. Therefore, the value of the original continuous-time signal 
at an arbitrary point can be evaluated by applying an interpolation filter (e.g., sine 
interpolation) to the samples. Hence, the design of the resampler has been 
transformed to the design of effective interpolation filters or, in other words, 
fi-actional-delay filters with variable delay. 

There are numerous methods for designing fi-actional-delay filters. These 
filters have diflFerent coeflScients for different delays. Thus, to implement variable- 
delay interpolation filters, one could compute one set of coefficients for each 
quantized delay value and store them in a memory. Then, in real-time, depending 
on the fractional delay extracted from the incoming signal, the corresponding 
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coeflBcients could be loaded. However, this method is likely to result in a large 
coefficient memory. 

To design low-cost modems, a large coefficient memory is undesirable. 
Gardner, et al, have shown that polynomials can be incorporated to compute the 
desired samples that are synchronous with the transmitted samples (Gardner, F.M., 
IEEE Trans. Comm, ^7:502-508 (1993); Erup, L., et aL, IEEE Trans. Comm, 
^7:998-1008 (1993)). In this case, an extensive coeflBcient memory is not needed. 
Moreover, the polynomial-based structure can be implemented efficiently with a so- 
called Farrow structure (Farrow, C, "A continuously variable digital delay 
element," in Proc, IEEE Int Symp. Circuits Syst (June 1988), pp. 2641-2645). 
This method is reviewed in Section 2.1. Although this approach achieves 
reasonable performance, the hardware complexity grows rapidly as the number of 
samples used to calculate each new sample is increased for better accuracy. In 
addition, given a fractional delay n, to produce a new sample using a degree A^-1 
polynomial, there will be A^-1 sequential multiplications that involve \i since we 
must compute \i raised to the (A^-l)-th power times a data value. Thus, the critical 
data path gets longer as A'^increases, thereby creating a limitation on the achievable 
data rate. 

Starting in Section 2.2, a new approach to interpolation is introduced. 
Instead of approximating the continuous-time signal with a conventional (i.e., 
algebraic) polynomial, a trigonometric polynomial is used according to the present 
invention. First, some background information is given. Next, the detailed 
implementation is discussed. We then evaluate and compare the performance and 
computational complexity of the algebraic polynomial interpolation to that of our 
method, giving numerical results. 

2. 1 Interpolation Using an Algebraic Polynomial 

To interpolate using function valuesX^) at A^equally-spaced sample points, 
also referred to as a "base point set," we can fit an algebraic polynomial of degree 
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NA to the data, as in FIG. 2. As explained in (Gardner, F.M., IEEE Trans, Comm. 

502-508 (1993)), there should be an even number of samples in the base point 
set, and the interpolation should be performed only in the center interval of the base 
point set. That is, Nis restricted to be even. In other words, given 4 samples points 
in FIG. 2 including y(-l), y(0), y(l), and y(2), the interpolation is performed at 
offset ji between the points y(0) and y(l) to determine the point 202 on the curve 



It seems that one would have to solve for the coeflBcients of the (A^-l)-th 
degree polynomial from these available samples before the synchronized (i.e., 
interpolated) samples can be computed. However, a method devised by Farrow 
(Farrow, C, "A continuously variable digital delay element," in Proc. IEEE Int 
Symp. Circuits Syst (June 1988), pp. 2641-2645) can compute the synchronized 
sample from the available samples efficiently with use of an algorithm that is well 
suited for VLSI implementation. To illustrate, we consider the following example 
of interpolation using a cubic Lagrange polynomial. Without loss of generality, let 
us assume the sampling interval is = 1 . Using the Lagrange formula for = 4, 
the synchronized samples can be computed as 



M=y{-\)C_X^i)^y{0)cXM)^y{\)cXpi)+y{2)cX^) (2.1) 



where C„(ji), n = -1, 0, 1, 2, are the third degree polynomials that are shown in 
FIG. 3. 

Obviously, 



Thus, y(n) in (2. 1), the sum of polynomials C„(n) weighted by the y(n) values, must 
be a third degree polynomial and must go through the samples y(- 1), y(0), y(l) and 
y(2). Writing C„(^) as 



P(t). 




1 ^ = nT,nan integer 
0 all other integers. 



(2.2) 



(2.3) 



k = 0 
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the coefficients of C„(n) are fixed numbers. They are independent of |i. We 
can re- write (2. 1) as 



3 ( 2 



n= - 1 Jt= 0 



Jt=0 ^n=-l 



Ar=0 



where 



10 



v{k)= Y.yi^)c^' (2.5) 

To minimize the number of multipliers, we can use a nested evaluation of (2.4), as 

y{ju) = ((v(3)// + v(2))// + v{\))jLi + v(0). (2.6) 
A Farrow structure 400 (for A'^= 4) that implements equations (2.5) and 
(2.6) is shown in FIG. 4. It consists of multiplications of data by fixed coeflBcients, 
and data multipliers, as well as addition operations. 



2,2 The Trigonometric Polynomial Method 



15 



20 



To solve the problems discussed in the section 2.1, the present invention 
utilizes a trigonometric polynomial to fit the asynchronous samples in FIG. 2. 
Using = e-'^"^ notation, for / e [-N/2 + 1, N/2], the polynomial may be written 
as: 



1 N/2~l 1 1 



. (2.7) 

The polynomial must cross the N samples. Thus, the coeflScients can be 
determined by solving the linear equations in TV unknowns: 



k=-N/2+\ ^ ^ 



(2.8) 



whose solution is 
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= z Anw:r^ k = -—+i,...,—. (2.9) 

The expression in (2.9) is simply the A^-point discrete Fourier transform 
(DFT). This suggests that, given iV^ equally-spaced samples, we can compute the 
DFT of these samples as in (2.9) to obtain the coeflBcients of the interpolating 
trigonometric polynomial in (2.7). Then, for a given offset \i, the synchronized 
sample >'(n) can be computed using that polynomial as: 



I N/2-l ^ 

Since Cj^ and c,j^ are conjugates, this equation can be simplified as 

1 r N/2-i > 



(2.10) 



(2.11) 



Flowchart 500 in FIG. 5 sunmiarizes the interpolation between two sample 
points at an offset \x. using a trigonometric polynomial, where the two data samples 
that are to be interpolated are part of a set of N-data samples (see FIG.2). The 
flowchart is described as follows. 
15 In step 502, a set of N-data samples are received having the two data 

samples that are to be interpolated. 

In step 504, coefficients of a trigonometric polynomial are determined based 
on the set of N data samples, according to equation (2.9). In doing so, the N data 
samples are multiplied by a complex scaling factor W^^^to generate a coefficient 
20 for the trigonometric polynomial, wherein Wn= e'^^"^, and wherein n represents an 

n* data sample of said N data samples. 

In step 506, the trigonometric polynomial is evaluated at the offset \i based 
on equation (2. 10) or (2. 1 1). 

In step 508, the real part of the evaluated trigonometric polynomial is 
25 determined and represents the desired interpolation value 202 in FIG. 2. 
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There are three issues that are to be considered in evaluating an 
interpolation scheme: 1) accuracy of interpolation, 2) complexity of implementation 
and 3) latency. In following sections, the trigonometric interpolation method is 
compared to the traditional polynomial method, particularly the Lagrange 
5 interpolator, in all these categories. 



2.3 Performance Comparisons 



10 



15 



20 



Let us first derive the impulse responses of the interpolators. With N 
samples, A^an even integer, the Lagrange formula (2. 1) is 

M/^)= X yin)cX^). (2.12) 

In addition, the interpolation is performed in the center interval. Thus 0^ |i ^ 1 . Let 
us define a new functiony(') such that 



/(//- w)= C„(//), 0<//<l, - N/2+l<n< N/2. 
Using the example of FIG. 3, defining t= \i-n,vfe have 



(2.13) 



( C,{t + 2) 

cJt-D 

. 0 

Thus, the Lagrange formula becomes 



-2<t<-l 
-l<t<0 
0<t<l 
l<t<2 
otherwise. 



(2.14) 



y{t)= Z y(rdC„(t)= S yin)f{t-ny 



(2.15) 



Therefore, the approach to reconstruct the continuous signal using the Lagrange 
polynomial is in fact equivalent to applying an interpolation filter to the 
available samples, withXO being a piecewise polynomial. The interpolator's 
impulse response obtained fi"om (2.14) is shown in FIG.6A. 
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Taking the Fourier transform offit), we obtain its frequency response. This 
allows us to evaluate the interpolation accuracy by examining the frequency 
response of the interpolation filter. The frequency response 702 of the Lagrange 
cubic interpolator (N = 4) is shown in FIG. 7A. The horizontal axis is the 
normalized frequency // Fy, with = 1 / T^. An ideal frequency response should 
have value one in the passband (0 ^ // < 0. 5) and be zero in the stopband (f/ Fs 
^ 0.5). 

For the interpolator using a trigonometric polynomial, we can express XO 
in terms ofX^) by substituting (2.9) into (2.10): 

1 Nn-\ f Nil \ 



1 



( Nil \ 
^n='Nll+\ 



COS TTJU 



(2.16) 



1 Nil ( Nll-l \ 

^ n=-NIM ^k=-Nll+l ^ 



— Z M 1+ 2 Z cos— -(//- /i) + cos4//- n) 



Jt=l 



Defining 



f NI2-\ 27ik 

l+2y COS— -r+ cos;jf k|<A''/2 
V 0 \t\>NI2 



(2.17) 



we have 



(2.18) 
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The impulse response (2. 17) is shown in FIG. 6B. The corresponding 

frequency response 704 of the trigonometric interpolator (for N=4) is shown in 
FIG. 7 A in thin lines. 

By comparing the frequency responses of the two interpolators, we can see 
that the trigonometric interpolation response 704 has a sharper roU-oflf in the 
transition band and more rapid attenuation in the stopband than the Lagrange 
response 702. These traits are enhanced as N increases, as demonstrated in 
FIG. 7B. For N=32, the trigonometric response 708 has a sharper roUofF than the 
Lagrange response 706 as shoAvn in FIG. 7B. 

Next we verify these observations by interpolating the samples of a practical 
time-domain signal. As an example, we interpolate a baseband signal with raised 
cosine spectrum and roll-off factor a = 1.0, sampled at two samples per symbol 
period, as shown in FIG. 8A. 

The interpolation accuracy here is measured as the normalized mean- 
squared difference between the signal interpolated with an ideal interpolator and 
the signal interpolated with the practical interpolator. The normalized mean- 
squared error (NMSE), discussed above, is calculated for both the Lagrange 
interpolator and the trigonometric interpolator for a range of typical values of N. 
The results are plotted in FIG. SB. 

Our test results demonstrate that the performance is improved with the 
trigonometric method. Using the same number of samples to interpolate, the 
proposed method produces a smaller NMSE, and the performance gain becomes 
greater as the number of samples increases. 

2, 4 Efficient Implementation Structures 

Recalling from Section 2.2, the trigonometric interpolation algorithm 
includes substantially two steps: 
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to 



10 



15 



Step I . Given a number of data samples N, calculate the Fourier coefficients 
Cf^k = 0, N/2 using (2.9). In a preferred embodiment, an even number of N data 
samples are used. In other embodiments, an odd number of data samples are used. 

Step 2. Compute the synchronized sample >'(|i) for any given according 
to (2.11). 

The first step involves multiplying the data samples by complex scaling 
factors W^^, Since these factors lie on the unit circle, the computation in Step 1 
can be simplified. Let us examine the case when 4: 

Example 2.1: ¥orN= 4, the Fourier coefficients are obtained as: 

Co = yi- 1) + ;^(0) + yd) + y{2) 

= b(0) - 3/(2)]+ y[->;(l) + 1)] (2.19) 
c^ = y{0)-y{\) + y{2)-y{-\). 
As seen in (2. 19), there is no nontrivial scaling multiplier required {orN= 4. 

Example 2. 2: 'We now compute the coefficients in (2.9) for N= 8. Using 
the trigonometric identities, we can obtain the following simple form for Ci,k = 0, 
...,4: 

Co = ^(-3) + y(-2) + y{- 1) + y(0) + y(\) + y(2) + >.(3) + y(4) 
q = {>^(0) - y{4) + [ -j/(-3) + >^(l) + >^(- 1) - j/(3) ] cos(;r/ 4)| 

+ 2) - y(2) + [ >^(-3)-;^l) + ;/(-!) -><3) ] cos(;r/ 4)} 
[-yi-^) + M - y(2) + y(4)} + j{-y(-3) + y(- 1) - ><l) + ;;(3)} (2.20) 
;^(0) - y{4) + [ >/(-3)-;<l) + (-;^(-l) + ;^3)) ]cos(;r/ 4) j 



Co = 



+ j\-y(-2) + y(2) + 



y{-^)-y{i) + {-y(-i)-y{3)) ]cosU/ 4) 



c, = -y{-3) + y(-2) - y(- 1) + y{0) - y(\) + y(2) - y(3) + y(A) 
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The only non-trivial scaling multiplications are those multiplications by cos 
(ttM). It appears that four such multiplications are needed to compute all the 
complex coefficients Cj^ A: = 0, 4. However, if we examine the data being 
multiplied by cos (n/4) (those terms embraced by the [] brackets), we observe that 
they are either the sums or differences of the [>(-3)-Xl)] and \y('\)-yO)] values. 
Therefore, we can compute {y('3)-y(l)] cos(7r/4) and I>(-l)-X3)]cos(7i/4), then 
use these two products to generate the c„ coefficients. Thus, only two scaling 
multiplications are needed in computing all the coefficients. 

Having observed the simplicity in the first step, let us focus on the second 
step. The second step may look complicated because of the complex 
multiplications cJV^'^^ and c^/^e'^''. However, since \W^'^^\ = le'^^l = 1, these 
products are just rotations of points Cj^ and c^^/j in the complex plane. Furthermore, 
this is the same type of operation performed in the phase rotation for carrier 
recovery by the phase rotator 124 that is shown in FIG. ID. This suggests that we 
can reduce the total complexity of the synchronization circuitry by sharing some 
resources needed by the digital resampler 122 and the carrier phase rotator 124. 
In one embodiment, a lookup table is utilized to determine the angle rotation 
associated with the angle for rotation of the c^ coefficients. In another 

embodiment, an angle rotator processor is utiHzed. Both embodiments are 
discussed further below, and the angle rotator processor is discussed in detail in 
section 5. 

\J^^^ PJQ;^^ illustrates an trigonometric interpolator 1000 that is one circuit 
configuration thatSfo^ements the trigonometric interpolator equations (2.9)- 
(2. 11), where the numberot^tet^^amples is N=4. The interpolator 1000 is not 
meant to be limiting, as those skilledTit^tii^^^s may recognize other circuit 
configurations that implement the equations (278)^(341) These other circuit 
configurations are within the scope and spirit of the present mention. 

The trigonometric interpolator 1000 receives input data samples having 
two data samples that are to be interpolated at an offset \i (see FIG. 2). The 
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resulting interpolated value y(n) represents the interpolated point 202 in FIG. 2. 
The interpolator 1000 includes a delay module 1004, an adder/subtractor module 
1006, and an angle rotator module 1008, and an adder 1012. 

The delay module 1004 includes one or more delay elements 1012. The 
delay elements 1012 can be configured using known components. 
\o\ Thri — frddnr/nubtrnrtormoflulr 1006 includes multiple adders (or 




subtractors) 1014, where subtraction is indicafea"te5ra--f-)..5igns. 

The angle rotator module includes two angle rotators 1010. The angle 
rotators 1010 can be configured using an angle rotator processor or a table lookup 
10 (e.g. read only memory ) as discussed below. 

The operation of the trigonometric interpolator 1000 is discussed fijrther 
in reference to the flowchart 1700 in FIG. 17, which is discussed below. 

In step 1702, the interpolator 1000 receives a set of N-input data samples. 
H The N-data samples include the two data samples that are to be interpolated at the 

1 5 offset |i relative to one of the data samples, as shown in the FIG. 2. In FIG. 2, the 

interpolation is to be performed between y(0) and y(l) at the offset |i to determine 
the interpolation value 202. 

In step 1704, the delay module 1004 delays the input data samples. 
\>^^>- Irrstep4706^tfie adder/subtractor module 1006 generates one or more 
20 trigonometric coefficients""ac5tJrdiagto the equation (2.9). In FIG. 10, the 

coefficients are represented by Cq, Ci, andC^^jF^I^, where the coefficients Cj 
is a complex coefficients. 

In step 1708, the angle rotators 1010 rotate appropriate real and complex 
coefficients in a complex plane according the offset \i, resulting in rotated complex 
25 coefficients. More specifically, the angle rotator 1 0 1 Oa rotates the real coefficient 

C2, and the angle rotator 1010b rotates the complex coefficient Ci in the complex 
plane. 

In embodiments, as discussed herein, the angle rotators 1010 are table 
look-ups. In which case, a rotation factor is retrieved fi-om the table lookup based 
30 on the offset |i, where the rotation factor includes the evaluated cosine and sine 
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functions of ^kfi that are shown in the equations (2.21) below. The rotation 

factor is then muhiplied by the corresponding real or complex coefficient, to 
generate the respective rotated complex coefficient. An interpolator 1800 having 
a table lookup ROM 1802 and a complex multiplier 1804 are shown in FIG. 18 
5 for illustration. 

In step 1710, the adder 1012 adds together Q, a real part of the rotated 
coefficient Ci, and a real part of the rotated coefficient Cj. The adder 1012 also 
scales the output of the adder as necessary according to equation 2.10. The 
resulting value is the desired interpolation value at the offset ^i, as represented by 
10 point 202 inFIG.2. 

Nj^^^^Tl^e trigonometric interpolator is not limited to the 4th degree embodiment 
that is shownhi^^IG. 10. The trigonometric interpolator can be configured as an 
jp N* degree interpolafe«;^ased on N-data points, as represented by the equations 

(2.9)-(2. 11). These other N^egree interpolators are within the scope and spirit 
t^" 15 of the present invention. For exmnple and without limitation, FIG. 1 1 illustrates 

ty an interpolotor 1 100 having N=8. Th^H^onometric interpolator 1 100 includes: 

% a delay module 1102, an adder/subtractoNmodule 1104 (having two scaling 

multipliers having coefficients cos (7i/4)), an an^e rotator module 1 106, and an 
adder 1 108(having an 1/8 scale factor that is not showi). The operation of the 
20 interpolator 1100 will be understood by those skilled in the arts based on the 

discussion herein. 

2. 4. 1 Using a Lookup Table 

For carrier recovery, the phase correction is generally accomplished by 
25 looking up the sine and cosine values corresponding to the phase, then by 

multiplying these values with the complex data. This requires the same operations 
as the rotation of by an angle ^ kju , that is: 
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Re[c^W-'"')= Re[c,)cos—kM- In{c^)s\n-^kM 

2 21) 
/m(c^)cos— ^//. 

The sine and cosine table can be used for both the resampler, as in (2.21), and the 
phase rotator for carrier recovery. In embodiments, a read only memory (ROM) 
is utilized as the lookup table. However, other embodiments could be utilized in 
including other types of memories. An interpolator 1800 utilizing a table lookup 
ROM 1 802 and complex multiplier 1 804 are shown in FIG. 1 8 for illustration. The 
ROM table access time is insignificant as compared to the computation time in 
other parts of the interpolator. Therefore, this method results in low hardware 
complexity and low computational delay. This implementation will be referred to 
as the table-lookup method. 



2. 4. 2 Using an Angle Rotation Processor 



When a very low complexity implementation is desired at the expense of 
a slight increase in computational delay, we propose to use an efficient structure 
for angle rotation, which is described in Section 5. Based on this structure, each 
angle rotator has a hardware complexity slightly greater than that of two 
multipliers. In addition, a very small ROM is needed. We will subsequently refer 
to this particular implementation of our algorithm as the angle-rotation method. 

Thus, there are at least two choices for implementing (2.21) as well as the 
phase rotator for carrier recovery. The trade-off is between complexity and speed. 
In a base-station where computation power is more affordable, the table lookup 
method might be used. In a hand-set, where cost is a major concern, an angle 
rotation processor might be used for both the resampler and the phase rotator, 
multiplexing the operations. 

Now let us compare the complexity of the trigonometric resampler with 
that of the Lagrange method. Table 2. 1 summarizes the comparisons for several 
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typical A'^ values. The numbers are based on the table-lookup method. It indicates 
that, for the same filter length A^the trigonometric interpolation method needs less 
hardware. 



Table 2.1 Complexity and latency comparisons. 
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2.5 Delays in the Computation 

'^"^^The-critical path ofthe Farrow structure 400 (FIG. 2) are now compared 
to that of the trigonometriG^jjtmjolator. The Farrow structure implements the 
Lagrange interpolator as discussed aBSve^^Jlie Farrow structure 400 is shown in 
FIG. 9(or FIG. 4), with the critical path 902 indicated,.;nie c^^^ path 902 for 
this structure includes one scaling multiplier 904 and - 1 datanitii^iers 906. 

In contrast, the critical path for the trigonometric interpolator 1000 is path 
1002 and it contains just one angle rotation 1010, or only one multiplier if the 
table-lookup method is employed to replace the angle rotator 1010. Since the 
angle rotations for various angles can be carried out in parallel, the critical path 
does not lengthen as increases. 

Table 2.1 compares the computational delays for the trigonometric 
interpolator with that of the Lagrange interpolator for various values of N. The 
delay data for the trigonometric interpolator 1000 are based on the table-lookup 
method. As shown in FIG. 10, the trigonometric interpolator (for N=4) has only 
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one multiplier in the critical path, whereas the Lagrange interpolator has 4 
multipliers in the critical path. Therefore, the trigonometric interpolator has less 
latency than the Lagrange interpolator, which is important for voice 
communications. 

2. 6 Simplifications of the Preliminary Structures 

As mentioned in Section 2.4, to produce each^^i) we first calculate the 
Fourier coefficients Cj^ using existing samples, according to (2.9). We then 
compute Re(c^I^^^) (2 11) This is accomplished either by 

retrieving W^^^ from a lookup table, followed by two real muhiplications, or by 
an angle-rotation processor. 

2. 6. 1 The Simplified Structure for N=4 

Let us examine the trigonometric interpolator 1000 having N==4. To 
compute Re(2c^^) and Re{c.^tV^^^) the system requires either two angle 
rotators 1004 or two accesses to a lookup table. 

If the input samples would happen to be such that 0^=0 then one fewer 
angle rotator, or one fewer ROM access, would be required. Of course, the 
original data samples X-l)^ XO), yO), andX2) are not likely to have the special 
property that c^=>'(0)->'(l)+y(2)->'(-l)=0. However, if the data samples are 
changed, then the modified samples can be determined to satisfy ^^=0. If the 
modified samples for interpolation, then the angle rotator 1010a can be 
eliminated. However, the interpolation result will not then correspond to the 
original data samples. 

It seems that the data samples can be changed to attain savings in 
hardware, as long as the interpolation result is fixed so that it corresponds to the 
original data samples. Of course, it must also cost less in hardware to "fix" the 
interpolation result than is saved in hardware by using the modified samples. 
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If the ssmplesyfk) are modified to form yik) according to: 
yi-l)=y{-l)-K 

y(i2)=yi2)+2K 
then the^ value can be adjusted to force the y{k) samples to satisfy 0^=0, where 

a: is the slope of a straight line 1202 in FIG. 12. 

To find K, the value that corresponds to the modified samples is 
determined according to: 

C2 = K0)-Ki) + K2)-K-1) 

= y(0) - W) +K)+ iy(2) + 2K) - (>(- 1) - K) (2.23) 
= 2K-{y(il}+y{-l)-y(0)-yi2)). 
To force c^=0, requires: 

K = ^iydl) + yi- 1) - J(0) - yi2)) (2.24) 

Therefore, the angle-rotator can be eliminated, and Cq and Ci are determined 
accordingly as, 

Co=2(Kl)+y(-i)) 

c,=[2X0)->'(l)-><-l)]+y[-2Xl)+>'(0)+><2)]. (2.25) 
Then, the interpolated sample is 

1 I f -^ 

y(M)=^c, + -K^[c/^'j. (2.26) 

However, >'(//) should be adjusted so that it corresponds to the original 
samples. From FIG. 12, the values expressing the difference between the original 
and the modified samples lie on the straight line 1202. From FIG. 13, it follows 
that the offset due to the modification of samples is ATji. Therefore, the 
37(//) value can be compensated by: 
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y{fi)^y{M)-K^. (2,27) 

Using equations (2.25), (2.26) and (2.27) leads to an interpolator 1400 
as shown in FIG. 14. This simplified interpolator structure is not limited to N=4 
configurations. In fact, this simplification technique can be applied to an 
interpolator with an arbitrary TV^ value. To eliminate the angle-rotation needed by 

Re(c^/2^''^^ ) in (2. 1 1), the samples are modified according to 

y{n) = y{n) + nK, n = integer. (2.28) 

Using (2.9), and then applying (2.28), results in 

A ( N/2 \ ^ f N/2 \ i 

<^«. = -^ I Z (-!)>(«) +V. (2.29) 

If we choose 

2 

^=-T7 Z (2.30) 

iV „=-N/2+\ 

we can force c^/2=0. 

Referring to FIG. 14, the interpolator 1400 includes the delay module 
1004, an adder/subtractor module 1402, the angle rotator 1010b, a multiplier 
1404, and an adder 1406, The interpolator 1400 has a simplified structure when 
compared the interpolator 1 000 (in FIG. 10), as the interpolator 1400 replaces the 
angle rotator 1010b with a multiplier 1404. As discussed above, this can be done 
because the coefficient Cn/2=0 (Cj^O for N=4) by modification of the data samples, 
and therefore there is no need to have an angle rotator for C^^. The operation of 
the interpolator 1400 is fiirther discussed in reference to the flowchart 1900 that 
follows. 

In step 1902, the interpolator 1400 receives a set of N input data samples. 
The N data samples include two of the data samples that are to be interpolated at 
the offset \x, as shown in the FIG. 2. 
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In step 1904, the delay module 1004 delays the input data samples. 
In step 1906, the adder/subtractor module 1402 modifies one or more of 
the data samples so that a coefficient c^^ is 0, In embodiments the data samples 
are modified according to y(n)„^^ = y(n) + n-K, wherein K is determined by the 
equation (2.30) above so that Cn/j is 0, and wherein y(n) represents the n* data 
sample of the N data sample set. 
^'^^^^^ 14 step 1908, the adder/subtractor module 1402 generates one or more 
tngonomefric coefficients according modifications to the equation (2.8). In the 
N=4 case, equations 2.25 are implemented by the module 1402. In FIG. 14, for 
N=4, the coe^cients are represented by Co and Cj, where the coefficient Ci is a 
complex coeffic^nts. By comparing with FIG. 10, it is noted that the Cj 
coefficient is 0. Additionally, the adder/subtractor module 1402 outputs the K 
value for further processing. Notice also that in FIG. 14, the output scaling factor 
has been changed fi-orti Vi to 1/2. This reflects several other straightforward 
simplfications that have be^en made to module 1 402 and in the angle rotator 1 0 1 Ob. 
In embodiments, the steps r906 and 1908 are be preformed simultaneously by the 
adder/subtractor module 1402, as will be understood by those skilled in the 
relevant arts. \ 

In step 1910, the angle rotator 1010b rotates the complex coefficient 
in a complex plane according the\. offset \x, resulting in a rotated complex 
coefficient. In embodiments, as discussed herein, the angle rotators 1 0 1 Ob is table 
look-up. In which case, a complex rotation factor is retrieve fi-om the table lookup 
based on the offset ji, and the resulting rotation factor is then multiplied by the 
corresponding complex coefficient, to generate the respective rotated complex 
coefficient. The rotation factor includes the eWluation of the cosine and sine 
factors that are shown in equations 2.21 . Note thaKsince C2 = 0, the angle rotator 
1010a is replaced v^th the muhiplier 1404. ^ 

In the step 1912, the multiplier 1404 multiplies the K factor by the offset 
H, to produce a K\i factor. 
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In step 1914, the adder 1406 adds together Cq and K^i and a real part of 
the rotated complex coefRcient Cj, and scales the sum by the trivial factor Vi, to 
produce the desired interpolation value. The addition of the K|i factor 
compensates the desired interpolation value for the modification that was done to 
the data samples in order to force 0^/2 to zero in step 1906. 
V/^O^ TliOs^mplified trigonometric interpolator is not limited to the four sample 
embodiment thatiS^-sljown in FIG. 14. The simplified trigonometric interpolator 
can be configured as ajT^S-s^ple interpolator based on N-data points, as 
represented by the equations (2.28H!2r3iJ). These other N-sample interpolators 
are within the scope and spirit of the present in^«ijtion. For example and without 
limitation, an interpolator with A^=8 is discussed belo^Sv,^ 

2. 6.2 The Simplified Structure for N^8 

According to (2.30), we choose 

K=^ (yil) + X3) + X- 3) + y(~ 1) - ^0) - y(2) - yi4) - y(-2)l (2.3 1) 

The coefficient values can be computed following this procedure: 
pi = (4^ + XI) -X-3)) cos 71/4 

^ (2.32) 
p2 = (4K + y(3) -y(-l)) cos n/4 

Co-2(Xl) + X3) + X-3)+M-l)) 

= [yiO) - yi4) -4K+pl- p2] + j[y(-2) - yi2) -4K-pl- p2] 
c, = [4K + >;(0) + X4) - y(2) - y(-2)] 4- j[4K + y(3) + y(- 1) - y(l) - y(-3)] 
c, = [y(0) - y(4) -4K-pl^ p2] + j[- y(- 2) + X2) +4K-pl- p2] 
C4 = 0. 

(2.33) 

The resulting modified structure for N=8 is shown in FIG. 15 as 
interpolator 1500. Similar to the interpolator 1400, the interpolator 1500 includes 
a delay module 1504, an adder/subtractor module 1506, an angle rotator module 
1 508, a multiplier 1510, and an output scaling adder 1512. As in the interpolator 
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1 400, the multiplier 1510 substantially replaces an angle rotator module. As in the 
interpolator 1 100 for A^^8 (FIG. 1 1), only two non-trivial scaling multiplications 
are needed for the modified structure 1500. 



2. 6.3 Performance Comparisons with Other Structures 

^^S^HQwdoes the simplified interpolator 1400 (FIG. 14) perform as compared 
to the interpolatoMQOO (FIG. 10)? FIGs. 1 6 A-C show the fi-equency responses, 
in solid lines, of the La^^Qge cubic interpolator 400 (FIG. #4), the interpolator 
1000 (FIG. 10) and the simplmelMnterpolator 1400 (FIG. 14), respectively. For 
an input signal whose spectrum is a rais^d>^sine with a=0.4, as shown in dashed 
1 0 lines, the amount of interpolation error correspohd.§to the gray areas. Clearly, the 

interpolator 1400 produces less error than the Lagran^i&^^iubic interpolator 400 
and the interpolator 1000 . (FIG. 16D will be discussed in Sec^teq^.) 

Next, let us verify this performance improvement by interpolating two 
signals: Signal 1, which is the same as signal 802 in FIG. 8 A (a=1.0), and Signal 
15 2, which is signal 1 602 in FIG. 1 6 (a=0.4). Then, the NMSE values are compared. 

We use three interpolators: 1) Lagrange cubic interpolator 400; 2) the 
trigonometric interpolator 1000; and 3) and the trigonometric interpolator 1400. 
Also compared is the number of multipliers. 

The results in Table 2.2 show that the modified structure for A'^ = 4 not 
20 only requires less hardware, it also obtains the highest accuracy among the three 

methods for these practical signals used in our simulation. 
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Table 2.2 Comparison of interpolators for A^=4. 



CO 
y ^ 



10 



15 



20 



25 



N=4 


Lagrange 
cubic 


Structure in 
FIG. 10 


Structure in 
Fig. 14 


NMSE for Signal 1 in dB 


-25.80 


-28.45 


-29.41 


NMSE for Signal 2 in dB 


-31.08 


-31.21 


-33.51 


Nontrivial scaling 
multipliers 


2 


0 


0 


Data multipliers 


3 




* 


Multipliers in critical path 


4 


2 


2 



1 lie UllgUllUillCUll^ UiLCl I^UidlUl l\J\jyj I^llipiU^a IW\J «J1X^-Vrui.put CUX^iW-xv/kCl.kVrii?, 

having the hardware equivalent of, slightly more than two multipliers. The 
trigonometric interpolator 1400 employs one such angle-rotator and one 
multiplier yielding an equivalent of slightly more than three multipliers. 



2. 7 Trigonometric Interpolator Application 

An important application of the interpolation method and apparatus 
described in this patent is the foUowdng. It is often necessary to increase the 
sampling rate of a given signal by a fixed integer factor. For example, a signal 
received at a rate of 1000 samples per second might need to be converted to one 
at a rate of 4000 samples per second, which represents an increase of the sampling 
rate by the integer factor four. There are methods in common practice for doing 
such a conversion. One method is a very popular two-step process wherein the 
first step creates a signal at the desired higher samphng rate but one where simple 
zero-valued samples (three of them in the example situation just mentioned) are 
inserted after each input data value. The second step in the interpolation scheme 
is to pass this "up-sampled" or "data-rate expanded" signal through an 
appropriately designed lowpass digital filter which, in effect, smoothes the signal 
by "filling in" data values at the previously inserted zero-valued samples. In the 
process of doing this filtering operation it might or might not be important that the 
original data samples remain intact, and when this is important there exist certain 
special lowpass filters that will not alter those samples. 
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e can easily adapt the trigonometric interpolator described herein to 
Efficiently c^ate such a sampling rate conversion system, but one that does not 
require such filtering operations. If we denote the integer factor by which we 
desire to increase the data rate as L (in the above example, L = 4) we proceed as 
follows. We build ^ system 7800 shovm below in Fig. 78A. System 7800 
includes a Delay Module\7802 and Add/Subtract Module 7804 (that are similar 
to that in FIG. 10), and suck that it can accommodate incoming data at a rate r. 
We now build L copies of the Ai^le-Rotation Module 7806 (similar to that in FIG. 
10), with each one being fed by t^e same outputs of the Add/Subtract Module. 
Within each of these L Angle-Rotaticm Modules 7806 we fix the value; that is, 
each one has one has a different one of tn^yalues: 1/L, 2/L, . . . , (1 -L)/L. With such 
fixed ]! values, each Angle-Rotation Modufe 7806 can be constructed as a set of 
fixed multipliers (a very special case of the table^Iookup method), although any of 
the Angle-Rotation Module implementations previously discussed can be 
employed. 

As shown in Fig. 78, the L-1 outputs, i.e., the intVpolated samples that are 
offset by the values 1/L, 2/L,..., (L-1)/L from the first\)f the two data points 
(indicated as p. = 0 and pi = 1 in the Delay Module of Fig\A) are routed to a 
multiplexer 7808, along with the input data point from whicRvall interpolated 
samples are offset. The multiplexer 7808 simply selects theste samples, in 
sequence, and provides them to the output at the expanded data rateSL x r. 

A major advantage of the system 7800 is that all of the system's 
components are operated at the (slow) input data rate except the output 
multiplexer 7808. If desired, it would also be possible to employ fewer Angle- 
Rotation Modules 7806, but operating them at a higher data rate, and using 
several |i values, sequentially, for each. This would resuh in a system that 
employed less hardware but one that traded off the hardware savings for a higher 
data rate operation of such modules. 
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2. 8 Trigonometric Interpolator Summary 

In^his Section we have described an interpolation method that we have 
i that os^s trigonometric series for interpolation. Comparing the 
interpolations using thetrigQnometric polynomial and the Lagrange polynomial of 
the same degree, the trigononielric^ased method achieves higher interpolation 
accuracy while simultaneously reducingth^omputation time and the amount of 
required hardware. Moreover, the trigonmrt©t|TC-based method preforms 
operations that are similar to that of a phase rotator fohs^ier phase recovery. 
This allows a reduction in the overall synchronization circuit corhplexity by sharing 
resources. 

This trigonometric interpolator yields less computational delay, as 
compared to algebraic interpolators. To achieve the same throughput rate, this 
translates into more savings in hardware using the proposed structure, because the 
data registers that are required by algebraic interpolators to pipeline the 
computation for a faster rate would not be needed by our structure. 

We have also introduced two implementations of the trigonometric 
interpolation method: one using a lookup table, and one using an angle rotation 
processor (to be discussed in Section 5). 

After introducing a first interpolation method, we have shown that we can 
trade one angle rotator for a muhiplier by conceptually modifying the input 
samples, then by "correcting" the interpolated value obtained from the "modified" 
samples. Through this modification, we have obtained a simpler implementation 
structure while simultaneously improving the performance when interpolating 
most practical signals. This performance improvement has been demonstrated by 
comparing the frequency responses of the interpolators and the mean-squared 
interpolation errors using these interpolators. Our discussion of the optimal digital 
resampler in Section 4 will be based on this simplified interpolator. 
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3. Interpolation Filters with Arbitrary Frequency Response 

In Section 2, an interpolation method using a trigonometric polynomial 
was introduced, along with an example of such an interpolation structure of length 
N=4. In addition to being a very simple structure, our analyses and simulations 
5 also demonstrated that the trigonometric interpolator outperformed the 

interpolator of the same length using a Lagrange polynomial. In this Section, a 
more systematic approach will be taken to analyze this method from the digital 
signal processing point of view. A formula giving its impulse response allows us 
j!; to analyze the frequency response of the interpolation filter. We then show how 

£^10 to modify the algorithm to achieve arbitrary fi*equency response. The discussions 

wis? 

fy in this Section will provide the framework for the optimal interpolator of 

Section 4. 



3,1 Formulating the Trigonometric Interpolator as an 
Interpolation Filter 



15 We have shoAvn that, given N equally-spaced samples y(n), a 

continuous-time signal can be reconstructed as 

N/2 

At)= Z (3.1) 

n=-N/2+l 

where f(t) is the impulse response of a continuous-time interpolation filter. As in 
Section 2, it is assumed that the sampling period is . This assumption makes 
20 the notation simpler and the results can easily be extended for an arbitrary T^. In 

other words, the invention is not limited to a sampling period of T=l, as other 
sampling periods could be utilized. In Section 2 we have shown that f(i) can be 
expressed as: 
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( 



NI2-\ 

\+2y\ cos— kt + cosTa \t\<NI2 



\t\> N/2. 



(3.2) 



0 



FIG. 20 illustrates / for the trigonometric interpolation filters for A^==8 and N=16. 
By computing the Fourier transform of /, we obtain the frequency response of the 
interpolation filter. The fi'equency responses for the A^=8 and A^=16 cases are 
plotted in FIG. 21. Since f(t) is real and symmetric around ^0, its fi-equency 
response has zero phase. In FIG. 21, the oscillatory behavior near the band edge 
is quite obvious. In addition, by comparing FIGs. 21a and b, we can see that as the 
filter length is increased fi-om N=S to N=16 the amount of ripple does not 
decrease. Well knovra as the Gibbs phenomenon, the magnitude of the ripples 
does not diminish as the duration of the impulse response is increased. 

It is apparent that the amount of oscillation cannot be reduced using the 
method discussed thus far while only increasing the filter length N. Moreover, it 
seems that an arbitrary fi-equency response cannot be achieved using this method. 
To address these problems, let us examine how the fi-equency response of this 
method is determined. 

3.2 Analysis of the Frequency Response 

Let us examine how the frequency response of / for the trigonometric 
interpolator is obtained using the example in FIGs. 22 with A^=8. According to 
(3 .2), the interpolation filter's impulse response f(t) on the interval - — </< — 
is a weighted sum of cosine fiinctions. We can view the finite-length filter fin (3.2) 
as having been obtained by applying a window 2206 according to the following: 




|f|< N/2 
\t\> N/2 



(3.3) 
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to an infinitely-long, periodic function 2204 with period N: 



f^(t)= 1+2^ COS — kt+COS7Ct -oo<r<QO 



N/2-1 



.k=l 



(3.4) 



such that 

fit)=M)w(,t\ -oo<r<oo. (3.5) 
Thus F the frequency response of f can be obtained by convolving and fTthe 
Fourier transforms of /, and respectively. 

The Fourier transform of the periodic function f/t)^ ^ </<<», is 

NI2-\ f OtT \ \ 1 

FA^)= E A^--^k]+-S{^-7t)^-S(Ci + 7u) (3.6) 

which consists of a sequence of impulses 2208. We will subsequently refer to the 
weights of these impulses as frequency samples. Denoting the weight of 



4° - f *) 



Q - by F{k) , we have 



M ( 271 ^ 



(3.7) 



N N 
where ^ - is an integer. In the case of (3.6), M=— . For our particular 

interpolation fiher, according to (3.6), all in-band frequency samples F(k) = 1 

for \k\ <NI2. For \h\ >NI2, the out-of-band samples F{k) = 0 . The two samples 

in the transition band are F(N /2)= F(- N 12)- 1 / 2 . The transition 

bandwidth is determined by the distance between the last in-band, and the first 

out-of-band frequency samples. 

-^m cg-w is ^ TfifM^n^nl^^ must be a sine function 2210. 

Convolving and the sine function iTsimplyTnteFpolatesJh frequency samples 
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^btain F(Q),-oo < Q < oo , shown as response 2212. (Here we have plotted 
the sj^innetric F only on the positive half of the Q axis.) We thus have 

lk)= F{Q)\^^^„„^, "T-^^y- ^^^^ 

From response 2212/tl^e continuous-frequency response F(Q) is uniquely 
determined by an infinite numb^of equally spaced frequency samples F{k) . If we 
modify the frequency samples 22P4ssnear the passband edge to let the transition 
between the passband and stopband be mere gradual, as depicted in FIG. 23, then 
the ripple is decreased. FIG. 23 demonstrates*^;adually reduced samples 2302, 
and the reduction of ripples in the overall response^2^4, as compared to the 
response 22 12 in FIG. 22. The cost of this improvement is, an increased transition 
bandwidth in the response 2304, as compared to the response z2J2. 

If a narrower transition band is desired, we can increase the duration of the 
filter /(1r/ This can be seen by comparing FIG. 24, where A^=16, with response 
2304 in FIG. 23, in which N=%. 

3. 3 Implementing the Modified Algorithm 

By comparing (3 . 6) and (3 . 4) we can see that the frequency sample values, 
i.e., the weights of the impulses in (3.6), are determined by the weights in the sum 
in (3.4). 

We can modify our original interpolation filter in (3.2) for |/| <, Nil as 

fm(O=H0)+2^ F(k)cos—kt. (3.9) 
By expressing (3.9) using the notation, for \t\ ^ N/2 we have 

Ut) = F(0)W°' + 2 F(k)(W^ + W-'^) = £ F(k)W-''. (3.10) 

andfjt)=0 for |/|>M2. Substituting this result into (3.1), and re-ordering terms, 
we have 
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K0=T7 Z ^.W-" (3.11) 



where 



= I y(rtw;^\. (3.12) 

N 

By comparing (3. 12) to (2.9), we can see that, for k= -— +1,..., — , 

= ^(^K- (3.13) 
Thus, a modified algorithm can be implemented as the following steps: 

Step I (Same as Step 1 in Section 2.4): Given an even number of samples 
N, calculate the Fourier coefficients Cf^ A=0,,..,A72 using (2.9). 

Step 2': Multiply the coeflBcients by scale factors F(k) using (3.13). 

Step i': Given a fractional delay ji, compute the synchronized samples 
using (3.11), which, due to = , can be simplified as: 



1 f ^ 

^(//)=-Re C0 + 2Z a.^*^ 



f M \ 

_ ■ (314) 

It seems that, in Step 3*, we need coefficients c{k) (hence not only for k^NIl 

but also for k>NI2 while only Cj^ values for k:^NI2 are computed in Step V. 
However, c^^ values for k>NI2 can be obtained using 

~ ^k^mN (3- 15) 

where /w is an integer such that 0:^k-mN:^NI2. We have (3.15) because Cj^ is 
periodic in k with period N, because Cf^ is obtained fi^om the Fourier transform of 
the discrete-time signal X'')^ 'N/2+1 ^n^N/2. 

At this point, we have shown that the continuous-time fi-equency response 
of the interpolation filter having impulse response f(t) can be improved by 
modifying the weights F(k) in (3 .10). Now a question arises: the modification of 
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the weights would alter the shape of the impulse response of the f(t) filter. How 
do we guarantee that the resulting filter does not change the original samples? 

3. 4 Conditions for Zero ISI 

In order for f(i) not to alter the original samples when used for 
5 interpolation as in (3 .1), it must have zero-crossings at integer multiples of the 

sampling period: 

•^^">={i «'/o. « an integer. ^^ '^^ 

^2 The well-known Nyquist condition for zero ISI (Proakis, J.G., Digital 

CO Communications, McGraw-Hill, New York, NY ( 1 993)) states that the necessary 

10 and sufficient condition for (3.16) is 



2 -2mi)= 1 -00 < Q < 00. (3.17) 

Since the filter's impulse response f(t) has a finite duration, i.e.f(t)=0 for |/|>M2, 
(3.16) holds if and only if the fi'equency samples F(k) satisfy 

CO 

^F{k-Nn)=l, A: = integer. (3.18) 

n=—oo 

1 5 The proof is given in Appendix A. 

In summary, we can still guarantee that the modified interpolation filter 
/ does not alter the original samples as long as the modified weights F(k) 
(fi'equency samples) satisfy (3.18). Using this constraint, we can design the 
weights F{k) to meet given fi-equency response requirements. 

20 3.5 Optimization Algorithm 

Using the approach discussed, one can approximate an arbitrary fi-equency 
response by choosing appropriate weights F(k) . For example, a desired fi-equency 
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10 



15 



response Fj{Cl) for an interpolation filter should be unity in the passband and 
be zero in the stopband, as 



71 ^ TV 

1 -— <n<y 
0 |qT> ^• 



(3.19) 



The interpolation error using our interpolation filter is defined as 

e(Q)= ^^(Q)|F^(Q)- F(n)| (3.20) 
where is a weighting fianction. 

From Section 3.2 we have F(Q ) = (Q ) ® 5/>2c(Q ). Thus, we can 



express F(Q) in terms of FQc) , using (3 .7), as 



® sinc(Q) 




(3.21) 



= ^ F(>t)sinc| 



k — M 



An optimal interpolation filter can be designed by choosing F{k) to 
minimize the peak interpolation error, as 



L„ = max{e(Q)} 



(3.22) 



or the mean-squiired interpolation error 



= r e^(Q)c?n (3.23) 

subject to the constraint described by (3.18). 

By examining FIGS. 23 and 24, we can see that, by modifying only two 
frequency samples, those nearest the band edge, a significant improvement is 
achieved. In these cases we haveF(A:) = 0 for |^1>A72+1 . 
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3. 6 Conclusion 

In this Section, an interpolation method was presented that achieves 
arbitrary frequency response by modifying the trigonometric interpolator discussed 
in Section 2. Using this approach, the performance of a trigonometric interpolation 
filter can be further improved. 

It is interesting to note that this procedure is equivalent to the well-known 
filter design method using windows. FIG. 25a depicts the impulse responses of the 
original filter (3 .2) as the dashed line, and the modified filter (3.9) as the solid line. 
By comparing the two impulse responses, we have found a function illustrated in 
FIG. 25b. If we multiply the original impulse response by this function, we get the 
impulse response that we obtained by modifying the frequency samples. Therefore, 
this function is equivalent to a window. According to this interpretation, our 
frequency domain design method is equivalent to designing a better window than 
the rectangular window (3.3) in the time domain. 

4. Design of Optimal Resamplers 

4.1 Motivation 

We have thus far discussed digital resampling using interpolation methods. 
To accomplish this, we conceptually reconstruct the continuous-time signal by 
fitting a trigonometric polynomial to the existing samples and then re-sample the 
reconstructed signal by evaluating this polynomial for a given sampling mismatch 
(or offset) \i. The reconstruction of the continuous-time bandlimited signal y{t) 
from existing samples >'(w) using interpolation filter according to (3 .1), is 



M0= S y{m)f{t-m). 



(4.1) 



m=-Nn+\ 
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Theny is resampled at / = as 

y(M)= Z yMf(M-m)= (4.2) 

where/^(w) =J(m - ^), 

In the previous sections we approached the problem only from the point 
of view of reconstructing the continuous-time signal, as in (4.1), since we have 
only examined the frequency response of the continuous-time filtery(/). However, 
what we actually are interested in is the new sample ^(^i) that is obtained by 
resampling the continuous-time signal at a new sampling instant, / = ji. 

Now, a question arises: Even when the frequency response F(Q) of the 
continuous-time filter is optimized as in Section 3 .5, do we necessarily obtain the 
minimum error in producing a new sample Xm) ^or each ^ value? 

According to (4.2), the new sample X^) is actually obtained by filtering 
the original samplesX'w) using a discrete-time filter whose impulse responsey|^(/w) 
depends on a particular delay \i. What is the desired frequency response of the 
f^(m) filter? 

A digital resampler that is used in timing recovery simply compensates for 
the timing offset in sampling the received signal. Ideally, the /^(/w) filter should 
not alter the signal spectrum as it simply delays the existing samples by [i. 
Obviously, the desired frequency response of the discrete-time f^(m) filter is 

where co is the normalized emgular frequency. Let us define the frequency 
response of f Jim) as Fj,((o). The error in approximating the ideal frequency 
response i\i((o, \i) by FJ^io) for a given \i value, is 

e{co) = WXa?)\FXco.ju)- ^^(^)| (4.4) 
where W^icD) is a weighting fiinction. 

We now examine how the discrete-time fractional-delay frequency 
response FJ<d) is obtained. We denote by F(Q) the Fourier transform of the 
continuous-time filter^/). Hence, the Fourier transform ofy(/-n) must be 
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e'^^F(Q). We know that f^(n) =f{n-\i) are just samples ofJ{t-\x\ where -oo < / 
< oo. Therefore, according to the sampling theorem (Proakis, J.G., Digital Signal 
Processing, Macmillan, New York, NY (1992)), the Fourier transform of fjin) is 



after we replace Q on the right-hand-side expression by the normalized angular 
frequency variable co (co = Q since = 1). This relationship is shown in FIG. 26, 
where F(Q) corresponds to the N = S interpolator of FIG. 11. As discussed in 
Section 3 .2,^(0 is symmetric around / = 0. This implies that F(Q) has zero phase. 
To make the^O filter physically reedizable, of course, we must introduce a delay 
ofN/2, where N corresponds to the length of the filter. However, this delay 
simply "shifts" all input samples hy N/2, which is an integer because is even, 
and it does not change the characteristic of the input signal. Thus, it does not 
influence the interpolation accuracy. Therefore, to simplify our notation, we just 
use F(Q) as a real function. Hence, the phase of the complex function e'^^F(Q) 
is -Q|i if F(Q) ^ 0, or -Q^ + tt if F(Q) < 0 — the phase depends on |a, FIG. 26 
shows that the frequency response of the discrete-time filter F^(co) is obtained by 
first making an infinite number of copies of e'^^F(Q) by shifting this function 
uniformly in successive amounts of 2n, then by adding these shifted versions to 
e'^^F(Ci). As a sum of complex functions, the shape of i^^(o>) depends not only 
on the shape of the continuous-time frequency response F(Q) but also on the value 
li. The dependence of F^(io) on ^ is illustrated in FIG. 27, where F^J[o:>) is 
obtained from the function F(Q) in FIG. 26, using = 0.12 and = 0.5. The 
magnitude of the ideal fractional-delay frequency response, defined in (4.3), is 
shown in both FIG. 27A and FIG. 27B as the dashed lines. It is evident that the 
frequency response 2706 is worse for ^i = 0.5 than the frequency response 2704 
is for = 0. 12, since the response 2706 deviates more from the ideal frequency 
response 2702 than does the frequency response 2704. Hence, the interpolation 
error is larger for = 0.5 than for ^ = 0. 12. We have observed in our simulations 
that the largest interpolation error occurs when fi = 0.5, i.e., when the interpolated 



CO 




(4.5) 
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sample is exactly in the middle of the two nearest existing samples. As \i 
approaches 0 or 1 (i.e., as the interpolated sample gets closer to an existing 
sample), the interpolation error becomes smaller. Moreover, the interpolation 
errors obtained for and 1 - are the same. 



5 4.2 Resampler Optimizations 

In Section 3, we analyzed the relationship between the weights p (k) in 
(3.9) and the frequency response of the interpolator. We have shown that we can 
enhance F(Q) by adjusting the p (k) values. In FIG. 23 {N^ 8), for example, the 
p (3) and p (5) values correspond to the magnitude of the pulses 2302 near the 
10 band edge. If we adjust p (3) and p (5) such that the transition between the 

passband and stopband is more gradual, we can achieve a better frequency 
response. 

To flirther improve the interpolation performance, we could take ji into 
account, by optimizing i%,(o>) for each \i value. As in Section 3, we could adjust 

1 5 F{k) near the band edge to change the shape of F(C1), for each value, such that 

the discrete-time frequency response FJ^iS), which is obtained from (4.5), best 
approximates the desired response of (4.3). 

As discussed in Section 3, to guarantee that the original samples are not 
altered using the modified interpolator, the weights P{k) should satisfy (3.18). 

20 WhenA^=8, for example, we modify F(3) and F{5) together in order to satisfy 

(3.18). Here, however, our goal of optimization is to obtain the best 
approximation for the sample corresponding to a specific \i. Hence we need not 
be concerned with the issue of altering the original samples in Section 3, where 
there is only one set of optimized weights for all values. 

25 Let us demonstrate this method using the example of A^=8. We chose to 

modify F(3) and F(4) . For each given ji, we search for the F(3) and F(4) 

values suchthat (4.4) is minimized. Wedenotesuch P(3) and F(4) valuesby /;(3) 
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and ^(4) respectively, since they are now also dependent on ^. FIGs, 28 A and 

B show the modifications to F(Q) for ji=0.12 and ^=0.5, respectively. The 
corresponding optimized F^,(io) functions are illustrated in FIGrs. 28C and 28D, 
respectively. 

To demonstrate the performance improvement, let us use this example: 
Given n=0, 5, we optimize F^((x>) for the signal whose spectrum is shovm in dashed 
lines 2902 in FIG. 29. Comparing the un-optimized fi-equency response 2904 with 
the optimized fi-equency response 2906, the modification clearly produces a better 
fi-equency response. More specifically, the response 2906 is flat in the fi-equency 
band where the power of the signal 2902 is concentrated, and the deviation from 
the ideal response mostly falls in the "don't care" band. 

Similar to Section 3, to implement this improved method, we first compute 
the coefficients Cj, from the existing samples as in (2.9). Then, given the \i value, 
wemultiply, e.g., for A^=8, the and values by ^(3) and (4) , respectively. 
Finally, we compute the synchronized sample y(^) using (2.11), where C3 and C4 
are replaced by ^(3) Cj and ^(4)c4, respectively. 

We can apply similar modifications to the interpolator v^th iS^=4. FIG. 
30A show the frequency response of the interpolator 1000, for ^1=0. 5, while 
FIG. 3 OB displays the resuhs of a modified interpolator 1000, where parameters 

F(l) and F(2) are optimized, for ^1=0.5, to maximize the interpolation accuracy 

for the signal whose spectrum is shown in dashed lines. As can be seen the 
optimized response 3006 is flatter in the part of the spectrum of the signal 3002 
where most of its energy is concentrated than is the un-optimized response 3004. 

The flowchart 3400 in FIG. 34 generalizes the optimization of the 
trigonometric optimization procedure. The flowchart 3400 is similar to the 
flowchart 1700, but includes the additional steps of 3402 and 3404 that are 
described as follows. 
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In step 3402, a factor is determined to adjust the frequency response 
of the trigonometric interpolator so that it is consistent with the frequency 
response of N-data samples and the offset |i. 

In step 3404, one or more of the complex coefficients are multiplied by the 
to modify the frequency response of the interpolator so that it is consistent 
with the input signal and the offset |i. 

The optimization routine can also be used with K-modified data samples 
that leads to the simplified interpolator structures of FIGs. 14 and 15. The 
flowchart 3500 illustrates the F^ factor modification in the context of the 
flowchart 1900. 

As will be shown in the section that follows, the steps 3402, 3404, and 
1708 can be combined into a single step if a table lookup is used to determine the 
rotation factor. In other words, the sine and cosine values can be multipUed by the 
F^ factor before they are stored in the ROM. 

In Section 2.6, we have presented an efficient algorithm that eliminates one 
angle-rotation. For example, for A^=4, we can "modify" the input samples 
according to (2.22). With this modification, we can treat the new samples as if the 
input signal satisfies C2=0. The remaining non-zero coefficients are Cq and c^. In the 

example for A^=4 in the previous Section, two parameters, (1) and (2) , are 

optimized to achieve the best approximation of a desired fractional-delay 
frequency response described by (4.3). Now, v^th C2=0, we have only one 

parameter, F^ (1) , to choose. 

The impulse response of the simplified interpolation filter is derived in 
Appendix B . From the mathematical expression of the impulse response (B . 5), we 
can obtain the corresponding frequency response. The frequency responses of the 

interpolator 1400 (FIG. 14) before and after appljang the F^ (1) modification are 

shown in FIG. 3 1 A-B, respectively. We can see an improved frequency response 
3106 over the response 3104, as the response 3106 is flatter in the part of the 
signal 3102 where its energy is concentrated. Furthermore, it seems that the 
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frequency response 3106, where only is modified (c2=0!), is as good as the 
modified response 3006 in FIG, 3 OB where both and are modified. 

4. 3 Implementations 

It may appear that additional hardware is needed to implement the 
multiplication by, for example, ^^(1) for the simplified N=4 structure. Let us 

examine the corresponding computations. As we know, we first compute 
coeflBcients Cq and according to (2.25) (c2=0, of course). We then compute y(\i) 
using 

yiM) = ^RQiF^(l)c/^')- Kju (4.6) 

according to (2.26) and (2.27), where K is defined in (2,24). As discussed in 
Section 2.4, the computation 

ReCF^COc/^") = Re(c,)(Re(F/l)/^)- lm(c,)Im(F^(l)e'^') (4.7) 

^ J — M 

can be accomplished by retrieving the F^(V)e ^ value fi-om a ROM lookup 

table and then multiplying Re(Cj)+j Im(Cj) by the retrieved value, since both 
FAV) and ^^2^ can be pre-determined for all values. 

If the angle-rotation method is used, we can use a lookup table to store the 
F^iX) values. In real-time, afl:er computing c^e^'^^ using an angle-rotation 

processor, we can multiply the result by (1) . In this case, if F^iS) is allowed 
to be a complex number in optimizing performance, we then need two real 
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multiplications to obtain RG(F^(l)cy ) in (4.6). However, if we restrict 
i^O) to be a real number, we can use just one real multiplication as 

Re(F^(l)qe'?^) = Re(qe'?^) (48) 

According to Table 4, 1, the NMSE using complex and real (1) values 

are -37.41 dB and -37.08 dB, respectively. Therefore, the performance 

degradation caused by restricting to be a real number is insignificant. 

When the table-lookup method is employed, the implementation structure 
for the optimal interpolator is the same as that for the interpolator 1400, except 
for the coefficient which is now multiplied by (l)e^^^ instead of ^-^f ^ 

^ J — 

The table should therefore contain the Ke(F^ (l)e ^ ) and 

^ ) values, rather than the sin— ju and cos— // values used by 

the interpolator 1400. We now show that the size of the table is the same as the 
one storing the sine and cosine values. 

Let us examine the contents of the lookup table. FIG. 32 displays the 
^Q(F^(V)e^^^) and lm(F^(l)e^ values, used by (4.7), where the real 
values are represented by curve 3202, and the imaginary values are represented 
by the curve 3204. These values are monotonic with respect to \i, just like the 
Sin Y and cos— // values for 0:^ ^ ^ 1 . Moreover, simulations show that, when 

optimal values of (l) are reached, the real and imaginary components of 

^ J — /i 

F^(l)e ^ display the following complementary relationship: 

Im(/;(l)e'?^) = Ro(F,_^(l)e-'^''-''). (49) 
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Therefore, we need only store one of the Re(/^(l)e ^ ) and 
Im(i^(l)e ^ ) values. The other can be obtained by looking-up the value 
corresponding to 1 This is the same technique used in storing and retrieving the 
sine and cosine values with the purpose of reducing the table size. 

Various circuit implementations of optimized interpolators having N=4 are 
illustrated in FIGs. 36-37, These circuit implementations are presented for 
example purposes only and are not meant to be limiting, as those skilled in the arts 
will recognize other circuit implementation based on the discussion given herein, 
including interpolator configurations having different N values. 

FIG. 36 illustrates an optimized interpolator 3600 that is based on the 
simplified interpolator 1400 ( FIG. 14). The interpolator 3600 includes an 

ROM 3602 and a multiplier 3604. The ROM 3600 stores the appropriate 
value indexed by ^. The multiplier 3604 multiples the complex coefficient 

by the appropriate F^ value, and therefore optimizes the frequency response of 

the interpolator 3600. As discussed above, the order of the angle rotator 1010b 
and the multiplier 3604 can be interchanged so that the rotated complex 

coefficient is modified by the F^ value. 

FIG. 37 illustrates an optimized interpolator 3700 that is similar to the 

optimized interpolator 3600, except that the F^ ROM 3602, the multiplier 

3604, and the angle rotator 1010b are combined into to a single ROM 3702, that 
Stores the Re(/^(l)e ^ ) and Im(^(l)e'^^'") values. Therefore, coefficient 
optimization and angle rotation are performed in a simultaneous and efficient 
manner. 

It will be apparent to those skilled in the arts that the combined F^ and 

angle rotator ROM 3702 can be implemented for interpolator configurations that 
include more than N=4 elements, based on the discussions herein. 
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4,4 Simulation Results 

We have verified the new design with the following simulation. A baseband 
signal, shown in FIG. 33, with raised cosine spectrum, two samples per symbol 
and 40% excess bandwidth was generated. Table 4. 1 compares the result for N=4 
using four interpolation structures: 1) the Lagrange cubic interpolator, 2) the 
interpolator 1000, 3) the interpolator 1400, 4) the optimal resampler using a 
complex ^ (1) value, and 5) the optimal resampler employing a real ^ (1) 
value. 

Using the optimal structure, the NMSE is reduced by 4 dB over the 
method without optimization (FIG. 14 structure). The performance is improved 
by more than 6 dB compared to the Lagrange cubic interpolator, while the 
hardware is reduced. Comparing the optimal structure to the FIG. 14 structure, 
a 4 dB performance gain was obtained without increasing the amount of hardware. 



Table 4.1 Comparison of interpolators for N=4. 



N=4 


Lagr. 
cubic 
struct. 


Struct, in 
FIG. 10 


Struct, in 
FIG. 14 


Optimal 
table- 
lookup** 


Optimal 
angle- 
rotat.** 


NMSE in dB 


-31.08 


-31.21 


-33.51 


-37.41 


-37.08 


Scaling multipliers 


2 


0 


0 


0 


0 


Data multipliers 


3 


3 


3 


3 


4 


Multipliers in critical path 


4 


2 


2 


2 


3 



a. Complex /^(l) values are used. The table stores the ^(l)e''^^ values, 
multiplied by ^^XlX^T^JusT^ne'raorej; is needed. 



IS 
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The frequency response of an optimized interpolator 1400 (FIG. 14) using 
a lookup table is shown in FIG. 16D. Also shown in FIG. 16A-D are the 
frequency responses of the Lagrange cubic interpolator 400, the interpolator 1 000, 
and the interpolator 1400 without optimization. The signal spectrum of FIG. 33 
is shown in FIG. 16D in dashed lines 1604. The interpolation error corresponds 
to the gray area 1606. From FIG. 16D, the performance improvement achieved 
by the optimal interpolator is evident because the gray area 1606d is a lower 
amplitude than the corresponding gray areas 1606a-c. In addition, these 
improvements are accomplished without increasing the amount of hardware. 

For a high-performance interpolator, we now turn to the structure 
described in Section 2.6.2, for ^=8. Applying a similar approach for A^=4, as just 
discussed, to the A^=8 interpolator of Section 2.6.2, we can multiply the 

coefficient by (3) , whose value optimizes the fi-equency response FJo^) of a 
fractional-delay filter with delay |i. 

In designing the proposed N=S interpolator, only one parameter ^(3) 

was adjusted to minimize the MSE in (3.22). 

Next, three interpolators of length N=8 were used to interpolate the signal 
in FIG. 33 : 1) a Lagrange polynomial interpolator, 2) a Vesma-Saramaki optimal 
polynomial interpolator (Vesma, J., and Saramaki, T., "Interpolation filters with 
arbitrary fi-equency response for all-digital receivers," in Proc. J 996 IEEE Int. 
Symp, Circuits Syst (May 1996), pp. 568-571) (with length 8 but a third degree 
polynomial) and, 3) the proposed interpolator. 

Table 4.2 shows the simulation results. These results demonstrate that our 
method has an NMSE more than 1 6 dB lower than the Lagrange interpolator, and 
more than 4 dB lower than the Vesma-Saramaki polynomial interpolator in 
(Vesma, J., and Saramaki, T., "Interpolation filters with arbitrary frequency 
response for all-digital receivers," in Proc, 1996 IEEE Int, Symp. Circuits Syst. 
(May 1996), pp. 568-571). 
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Table 4.2 Performance comparison. 





Lagrange 
interpolator 


Vesma- Saramaki 
interpolator 


Proposed 
interpolator 


NMSE in dB 


-45.29 


-57.34 


-62.17 


Scaling multipliers 


25 


16 


2 


Multipliers 


7 


3 


7 



4,5 Conclusion 

Instead of optimizing F(Q), the frequency response of the continuous-time 
interpolation filter, we could optimize F^(co) of the fi-actional-delay filter for each 
\i value. By doing this, better interpolation performance can be achieved, as 
demonstrated by the simulations. 

As for the implementation complexity, when a table-lookup method is 
employed, the optimal interpolator does not require additional hardware, just table 
values that implement the coefficient optimization and angle rotation. When the 
angle rotation method is used, one additional real multiplier is needed. 

For ^=4, the optimal interpolator attained a 6 dB lower NMSE than the 
Lagrange cubic interpolator, while requiring less hardware. 

5. A High-Speed Angle Rotation Processor 

In previous Sections, an interpolation method and apparatus for timing 
recovery using a trigonometric polynomial has been discussed. The major 
computation in this method is the angle rotation, such as angle rotator 1010 (in 
FIG. 10 and FIG. 14). As mentioned in Section 2.4, these operations, together 
with the phase rotator for carrier recovery, can be implemented by table-lookup 
in a ROM containing pre-computed sine and cosine values, followed by four real 
multipliers to perform the angle rotation (see FIG. 1 8). Herein, going forward, this 
approach will be referred to as the single-stage angle rotation. Although fast 
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angle rotation can be achieved with eflBcient multiplier design techniques, for 
practical precision requirements, the ROM table can be quite large. For 
applications where low complexity and low power are the major concern, can we 
further reduce the amount of hardware for angle rotation with slightly more 
computational delay? 

There are various hardware designs that accompUsh angle rotations, 
notably the CORDIC processors (Ahn, Y., et al, "VLSI design of a CORDIC- 
based derotator," in Proc. 1998 IEEE Int Symp, Circuits Syst, Vol II (May 
1998), pp. 449-452; Wang, S., etal, "Hybrid CORDIC algorithms," lEEETrans, 
Comp, ^6:1202-1207 (1997)), and, recently, an angle-rotation processor 
(Madisetti, A.,e/a/., "A 100-MHz, 16-b, direct digital frequency synthesizer with 
a 100-dBc spurious-free dynamic range," IEEE J. Solid-State Circuits 3^:1034- 
1043 (1999)). These algorithms accomplish the rotation through a sequence of 
subrotations, with the input to each subrotation stage depending on the output of 
the previous stage. In these cases, the latency is proportional to the precision of 
the angle. 

We now propose a diflferent approach for angle rotation. Here the rotation 
is partitioned into just two cascaded rotation stages: a coarse rotation and a fine 
rotation. The two specific amounts of rotation are obtained directly fi-om the 
original angle without performing iterations as does CORDIC. The critical path 
is therefore made significantly shorter than that of the CORDIC-type methods. In 
addition, only a small lookup table is needed. 

In this Section, methods and apparatus for two-stage angle rotation will 
be described. These method and apparatus are meant for example purposes only, 
and are not meant to be limiting. Those skilled in the arts will recognize other 
methods and apparatus for two stage angle rotation based on the discussion given 
herein. These other methods and apparatus for angle rotation are within the scope 
and spirit of the present invention. 

It will be shown that more precision and less hardware can be obtained 
using the two stage angle rotator compared to the single-stage angle rotator, with 
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slightly more computational delay. We will then show that, given an overall output 
precision requirement, various simplifications can be applied to the computations 
within the two stages to reduce the total hardware, 

5,1 The angle rotation problem 

If we rotate a point in the X-Y plane having coordinates (Xq, Yq) 
counterclockwise, around the origin, by the angle (|), a new point having 
coordinates (X, Y) is obtained. It is related to the original point (X^, Yq) as: 

X=XqCOS^ -Y^sin^ 

Y = YqCOS^ -fXosini (5. 1) 

5. 1, 1 Single-Stage Angle Rotation 

The operation in (5.1) is found in many communication applications, 
notably in digital mixers which translate a baseband signal to some intermediate 
frequency and vice versa. In addition to accomplishing (5.1) with CORDIC, a 
very common implementation is to store pre-computed sine/cosine values in a 
ROM (Tan, L. and Samueli, H., "A 200-MHz quadrature frequency 
synthesizer/mixer in O.S-^im CMOS," IEEE J, Solid-State Circuits 50:193-200 
( 1 995)). Then, in real-time, the computation in (5 , 1 ) is accomplished v^th a ROM 
access for each given 4> followed by four real multipUcations. This method avoids 
the excessive latency of the iterations performed by CORDIC and can jdeld lower 
latency than the angle-rotation method (Madisetti, A. , "VLSI architectures and IC 
implementations for bandwidth efficient communications," Ph.D. dissertation. 
University of California, Los Angeles (1996)). Furthermore, a very fast circuit can 
be built, based on efficient multiplier design techniques. However, since the size 
of the ROM grows exponentially with the precision of (|), a rather large ROM is 
required to achieve accurate results. 

ROM compression can be achieved by exploiting the quarter-wave 
symmetry of the sine/cosine functions and such trigonometric identities as 
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sinO = cos(7i/2-0). The angle (|) in the full range [0,27r] can be mapped into an 
angle 6e[0,7i/4]. This is accomplished by conditionally interchanging the input 
values and a;, and Yq, and conditionally interchanging and negating the outputs 
and rvalues (Madisetti, A., "VLSI architectures and IC implementations for 
5 bandwidth efficient communications," Ph.D. dissertation. University of California, 

Los Angeles (1996)). Thus, we will focus only on 0e[O,7r/4] and replace <|) by 6 

in (5,1). Defining 0 = (tc /4)0 , we must have 6 e [0,1]. 

Obviously, the sine/cosine ROM samples must be quantized because of the 
limited storage space for sine/cosine samples. This produces an error in the ROM 
10 output when compared to the true (unquantized) sine/cosine value, which will 

yO subsequently be referred to as the ROM quantization error. Next we examine how 

fn 

f y this error affects the output. Let cos 6 and sin 6 be quantized to N bits, to 

become [cos6] and [sinB], respectively. We have 

U cose = [COSe]+ A^,e 

n 15 where A^^q and A^q are the ROM quantization errors, which satisfy | A \ < 2 

O and I A gi^Q I < 2'^ . The error in X due to the ROM quantization is the difference 

between Jfcalculated using infinite-precision sine/cosine values and the quantized 
values, that is 

A;, = (X^cosQ - 7o sine)- (;^o[cos9]- 7o[sine]) 
20 Its upper bound is 

|A;,|<(|A^oK|i^ol)2-^- (5-4) 
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5. 1. 2 Rotation by a Small Angle 
If the rotation angle happens to be so small that 

then its sine/cosine values can be approximated as 

sin0«0 

COS0 « 1-(0V2), 
For such 0 no table is needed. Next, we show how accurate (5.6) and (5,7) 
estimating their approximation errors. 

The Taylor expansion of sin0 near 0=0 yields 

sin' t ^ , 

sine = e-— ^9' 

o 

where l, = hQ,0 <. h <. L Thus, since 

lsin'^| = |cos^|< 1 
and in view of (5.5), an error bound on (5.6) is 

lA3i„el=|sine-e|<|0^/61< 2-^/6. 
Similarly, the Taylor expansion of cos0 yields 

1 , cos^^4 

cosG = 1-t6 +-:^7^9 • 

2 24 

Thus, an error bound on (5.7) is 

|A^3el = |cos0-(l-eV2)|<|eV24| 
which is negligible in comparison to the bound on 

5. 1. 3 Partitioning into Coarse and Fine Rotations 
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While it is unlikely that (5.5) is satisfied for a given eG[0,7i/4], i.e., 
0 G [0,1] , if we let B>N/3 be the number of bits in 0 , then we can express 

9 = 0 A/ + 6 z, where 

0^ = c/i2-^...+ c/^/32-^'' (5,13) 
5 ^, = c/^,3,,2-^^^V...^^2-^ (5.14) 

with d,e{0,l}. Next we define = (tc /4)Qj^ and 0^ = (ti /4)0^ . Clearly, 
with 7i:/4<l and fi-om (5.14), 6^ satisfies (5.5). 

If we substitute 6 = 0^ + 6^ for <> in (5 . 1) and expand 005(6^^^ + 6 J and 
1 0 sin(6j^^ + e J, we obtain: 

X= A^, cos^j, - I^sin^^ 

r= l^cos^^ + Xjsin^^ ^^^^^ 

and 

JCi = cosO^ - Xy sin0^ 

yi= yocos0^+ ^oSin^A^. ^^^^^ 
Now the rotation (5.1) is decomposed into two stages: a coarse rotation 
15 (5 . 1 6) by 0ji^ followed by a fine rotation (5 . 1 5) by 6^. With this partitioning (5.5) 

and (5.6) can be applied to the fine stage: 

X= X,{\-Ql/2)-Y,Q^ 

7= 7i(l-0j/2)+^i0^. (^-^^^ 

A benefit of this partitioning is that the fiinctions cosB^^^ and sin0j^ in (5 . 1 6) 
depend only on the A73 most significant bits of the angle 0 , where 
20 0 = (tc / 4)0 . They can be stored in a small lookup table. This results in a 

significant ROM size reduction. However, the approximation (5.6) introduces 
additional error. We now seek to achieve an overall precision comparable to that 
in the implementation having one stage and a large ROM table. 
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Defining the approximation errors A ^j^jq^ = sin9 ^ - 0 ^ and 
^ cosG J- = <^os6 i,-(l-0^/2), and neglecting terms that are products of error 

terms or products of an error term and sin0^, which is always small, we calculate 
the total error in X as the difference between X calculated using (5 .15) and (5.16) 
and X calculated using quantized sinOj^^ and cosOj^^, and in (5.16) and using 
(5.17) instead of (5. 15). We obtain: 

= ^o(Acos^?^ cos^^ + A, ,3^^ cos^^ - A^^^^^ sin^^^^) 
- J'o(Asin^?^ cos^^ + A,,3^^ sin^^ - A3i„^^ cos^^). ^^^^^ 
Comparing this error estimate with (5.3) and (5.4) it is evident that, so 
long as the errors due to A ^^^^ and ts.^^^ are sufBciently small, the error ts^^ 

in (5.18) can be made comparable to that of (5.4) by reducing the A^^^g^ and 

A ^^0^ values, i.e., by increasing the number of bits in the sine/cosine samples 

stored in the ROM. For example, if we add one more bit to the sine/cosine 
samples, then lA.^^e^K 2""^'' and |A,i„0^|< 2-'^"V Therefore, from (5. 18), 

we have 

I A ^ I < I A^o I (2-^-' + 2-^^^^ / 24 + 2-^ / 6) 

+ 1 1^1(2-^"^ + 2-^^^V24+ 2-^ /6) (5.19) 
= (|A-ol+|5^ol)2-'^(l/2+ 1/6+ (1/24)2-^'^) 

which is smaller than (5.4). A similar relationship can be found for A^. This 
demonstrates that, if we add one more bit of precision to the ROM for the coarse 
stage, we can achieve the same precision as that in the one-stage case, but with a 
significantly smaller ROM. 

A straightforward implementation of this method is illustrated by the angle 
rotator 3800 in FIG. 38. The angle rotator 3800 includes a ROM 3802, butterfly 
circuits 3806 and 3810, and fine adjustment circuit 3804. 
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The ROM 3802 stores the cos Qj^ and sin 0^^^ values, where Q^^ is the most 
significant part of the input angle 6. In embodiments the input angle 6 is 
normalized and represented by a binary number, so that 0^^^ is the most significant 
word of the binary number, and 0jr is the least significant word of the binary 
number. 

The first butterfly circuit 3806 multiplies the input complex number 3812 
by the (cosOj^^)^ and the (sin Q^Jf to perform a coarse rotation, where the ( ) ^ 
denotes that the appropriate ROM quantization errors have been added to the cos 
0j^ and sin Qj^, by the adders 3814. 
v^I^> The,^neadjustment circuit 3804 generates a fine adjust value (1- 0^'^, 
where 0^ is the leaSt-significant word of the input angle 0. 

The second butterflySh^cyit 3810 multiples the output of circuit 3806 by 
0^^ and the fine adjustment value fi*om cu'lHut3804, to perform a fine rotation that 
results in the rotated complex number 3814. The^rfth^O^^ denotes the an error 
value A ^ has been added to improve the accuracy of the fin^Tet^on. 

The three error sources A « , A and A -r, are shown. The much 

smaller error source A ^^^g^ has been neglected. The thick line depicts the path 

along which the ROM quantization error ti^^^^^ propagates to X. The error 

A is multiplied by and then by cos0^ as it propagates along this path to 

become A cosej^^ cos 6 ^ when it reaches the output. This matches the error term 

in (5 . 1 8) obtained fi"om our calculation. In subsequent discussions we will use this 
graphical approach to find the error at the output due to various error sources. 

The ROM table 3802 in the rotator 3800 contains many fewer sine/ cosine 
samples in comparison to the number of samples needed to implement (5 .1) using 
a conventional (single stage) table-lookup approach. Although the approximation 
(5.6) introduces additional error, so long as that error is smaller thsm the 
conventional ROM quantization error, we can increase the precision of the 
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samples in our small ROM table such that, overall, precision is not sacrificed. In 
principle, we can reduce the hardware complexity significantly in one block of our 
structure, with the corresponding accuracy loss compensated by higher precision 
from another block, and at the cost of a slight increase in the complexity of that 
block. As a result, the complexity of the overall structure is reduced without loss 
of accuracy. We will now exploit this idea again to further reduce the 
computational complexity. 



5.2 Simplification of the Coarse Stage 



The coarse stage, according to (5 .16), involves muhiplications of input 
dataJi^o and by the cos0j^^ and sin6j^ values. Writing %mQj^ as the binary number 

sine^ = 0,b,,..b^r,b^^^^.,. (5.20) 
where 6„ e {0,1}, we now round sinS^^^ upward, to obtain an (M3H-l)-bit value 

[sinOJ, as 

[sin0^] = 0.6,...6^;3V.i ^2<^^^'\ (5.21) 
where J\r represents the number of bits in the real part and the imaginary part of the 
input complex number. In other words, the real part has A'^bits, and the imaginary 
part has Ambits. 

Letting 0i be the angle for which 

sinGi = [sinGj (5.22) 
we must have Q^^Q^^. Next we can compute the corresponding cos 0i value. 
Using sin0i = [sinO^^^] and cosOi values, we actually rotate the point having 
coordinate (Xq, I^o) by 0i, instead of 0^^^, as 

= cose 1 - 7o sine , 
7i = Fo cose 1 + ^0 sine p ^^'^•'^ 

Since 0i = arcsin([sin0^]) and, of course, 0^^^ = arcsin(sin 0^^), applying the mean 
value theorem we have 

, . f ' ~ ^ = asin' ^ = -r=i== (5.24) 
[sin^^]-sin^^ 
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where ^ satisfies sin 6^ ^ [sin 0^^^] . Since sin 9 ^ 1/ (V2 ) , according to 
(5.21) we must have 

^ < [sine^]< sinG^ + 2-^^'^^> < 2"^^'^^\ (5.25) 

For most applications, A^^9. Thus, according to (5.25), we have 0.7696. 
Applying this value to (5.24), 

9i-9a^^ Vl- 07696^ ' ([sine^]- sinG^) < 1.566 x 2''"''^'\ (5.26) 
Because 0 ^ 0^^^ ^ttM, we have, for ^ 12, that 

0 ^ 01 :e0,0978 + tz/4 = 0.8832. (5.27) 

The resulting fine-stage angle is 0-0i, instead of 0^ = 0-0^ Thus, as in 
(Madisetti, A., "VLSI architectures and IC implementations for bandwidth 
eflBcient communications," Ph.D. dissertation. University of California, Los 
Angeles (1996)), a modified fine-stage angle compensates for a simplified coarse- 
stage angle. Since sin0i = [sin0ji^], by rotating by 0i, the (M3+l)-bit number sin0i 
decreases the number of partial products needed in computing Xq sin0i and Yq 
sin0i to just over a third of those needed for sin ^ and Yq sin This 
simplifies the computation in (5.23). However, if we can also reduce the multiplier 
size in computing Xq cos 0^ and FoCOsOi, we can fiirther simplify (5.23). 
Certainly, truncating the cos0i value would reduce the number of partial products 
in computing Xq cos 0i and Yq cos 0i. Let us truncate cos Oj to 2N/3 bits to 
obtain [cos 0i]. Then, 

0< A,„,e. = cose, - [COS0,] < 2-^^'' (5.28) 

We now have 

X, = ^o[cosej-7oSine, 

y; = yo[cos0,]+ Xosme,. (5.29) 
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Apparently, by truncating cos0i, smaller multipliers are needed. But the amount 
of rotation is no longer Gj. We now examine the effect on of using the 
truncated value [cos6i] instead of cosOj as 

First, we determine how 0„ is different from 6^ due to the truncation of cos0i. 
Letting cosGi and [cos0i] denote specific values of a variable z we consider the 
function 

^ sin9i 
0(z) = atan ^ . (5.31) 

Hence, 0^ and 0„ are the 0(z) values corresponding to Zj = cosGj and 
Z2 = [cos0i], i.e,, 01 = 0(zi) and 0„ = ©(zj) . According to the mean value 
theorem, we have 

0(Z,)-0 (Z,) 

-=Q'(i) (5.32) 



or 



9i-e^ ,sine, 

1 — — = atan' — - — (5.33) 

COS01 - [cosGj] ^ ^ 

where [cosej s ? :s cose,. The negation of the derivative atan' — satisfies 



-atan'^=- ^^^i 



1 + 



' sinOA ' (sin^,)' + ~ (.sinO^Y + [cos^*]^ ' ^^^"^^ 

According to (5.27), for ^ 9 we have 

0 £ sine 1 < 0.7727 
0.6347 < COS0 , ^ 1 (5 35) 

0.6406 <[cosei]< 1 



1904.0140002 



-72- 



1:1 \ 

£0 



Since [cosG J = COS0 ^ - A ^^^^ from (5.28) we also have 
(sm0,f + [cos0,f = (sm0,y + (cos6^, - A 



2 ^ 1 ^~2N/M (5.36) 



= 1- 2cos^iA,^^ + A^coso, > 1- 2 
Combining (5,35) and (5.36), one can verify, fovN^ 9, that (5 .34) satisfies 

- atan' < 0.7976 . (5.37) 
Thus, according to (5.28) and (5.33), we have 

e^-d,< (cos^, - [cos^,]) X 0.7976 < 0.7976 x (5.38) 
Combining (5.26) and (5.38), and ioxN^ 9, we have 



ry 0 < 6>„ - 6^ < 1.566 X 2-^"^''"^^ + 0.7976 x 2 '^'' < 0,8827 x 2 ''''. (5.39) 

ffi This is the amount of coarse rotation error, due to coarse-stage simplifications, 

10 that a modified fine stage must compensate. Let us examine the bound on the fine- 

O stage angle. 

p Now, the fine rotation angle is 6/ = 6 - 0„ instead of 6^. If 0, satisfies 

^ \d\<2-'"^ (5.40) 

then we have |sin 6, - 0,| < 2"^/6. That is, the approximations sin 6, « 0, and 
15 cos0,« 1-0,^2 can be applied as discussed in Section 5.1. Let us now examine 

the bound on 6;. By definition, 

0< = < 0.7854 x 2 "^^'. (5.41) 
Therefore, subtracting (5.39) fi-om (5.41) yields 

- 0.8827 X 2-^^' <0,- {O^ - 0^) < 0.7854 x 2"^'^ (5.42) 
20 which implies (5.40) because 

0, = ^-^. = 0M + O.-O„=0,-(O„- d^). (5.43) 
Hence, no lookup table is needed for the fine stage. 
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Next, we examine the magnitude of the complex input sample after 
rotation. One can verify from (5.29) that 

+ Y,' = {X', + y:)[[cos0,J + (sin^,)'). (5.44) 

Since [cos6j] is obtained by truncating cos 61, we must have 0 i [cos SJ ^ cos 
0i, thus 

[cos (9, ]' + (sin 0, y < (cos 0, )' + (sin 0, J = 1. (5.45) 

Therefore, 

X^^Y,'<X',+Y,\ (5.46) 
To maintain the magnitude, the result and must then be multiplied by 

1 / ^[cos^jp + (sin^j ) . We define a new variable ^5|^^^ j such that 



2 COS^i] 



Since -^(cos 0^ )^ + (sin ^, ) = 1 , and [cos0i] is very close to cosSi because of 



(5.28), we have that ^[cos^J^ + (sin^J^ is very close to 1. Thus, the 
^[cos^ ] ^® small. We now examine the bound on S^^^^^ j . We can 



write 1 / -^[cos^i]^ + (sin^j) 



2 

as 



(cos ft y + (sinfty 

— (5.48) 
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Substituting (5.28) into (5.48), we have 



^ (5.49) 



,2 

' ecn«| 



cos^,]%(sm^,)^ 



y[cos^]%-(sin^,) 
Because (5.28) and (5.35) imply that A^^,^ << [cos^J, we have 

A^^^ « [cos^,]a^j„^ , hence we can omit A^„,^ in (5.49). Defining 

fcosfllA - 
[cos^J +(sin^J 



then (5.49) becomes Vl+ 2S . From (5.28) and (5.35) we must have 8^0. 
Applying the mean-value theorem to V1+ 2S , we have 



Vl+2^- yf[+0 1 
= / ^ 1 (5.51) 

^-0 71+24- 

where 0 ^ C ^8 Hence, 

Vl+2^< l+S (5.52) 

According to (5.35), 

fcos^.l 

^ <1 (5.53) 



[cos^i]^ + (sin^j) 
and therefore, from (5.28) and (5.50), we have 0 ^ 6 < 2■^'^ 

By definition, in (5.47), we have VI + 2^ =1 + , . Thus 

3 \ /J ' [cosOi] [cos6l] 

S, Hence . is bounded by 

0<^, .,<2-'^'^ (5.54) 

According to (5 .40) and (5 .54), instead of storing the sin0^ and cos0j^^ 
values in ROM, we may store sinG^, which has N/3 + 1 bits for each sample, and 
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[cosei], which has 2A^ / 3 bits. Given 6^ , the sinOj and [cosej values are 
retrieved from the ROM to be used in performing the coarse rotation. Since the 
actual angle 0„ differs from the 0^^^ = (7i/4) Oj^ , we must also store the - 0^ 

values, so that the fine stage can compensate for the difference. The 
approximations (5 .6) and (5.7) still apply to 0;, in view of (5 .40). In addition, the 
change of magnitude in the rotation using the sin0i and [cos0i] values, as seen in 
(5.45), must also be compensated. Therefore we store the S, , values in order 

[cos^J 

to scale the coarse-stage output by 1 +<X i 

[cos61 J 

We can now implement the coarse rotation stage as in (5.29). Later we 
will show that the scale factor 1 + ^ can be combined with the scaling that will 

be done for the second stage (i.e., the fine stage) at its output. 

To compute 0; we must first convert the normalized 0^ value to the radian 
value 0^, which involves a multiplication by tt / 4. Since 7i / 4 = 2"^ + 2"^ + 2"^ + 
2-8 + 2-13 ^ ^ jf multiply 0 ^ ^ < 2"^^^ by (2-* + + 2"^ + X% this product 

and (71 / 4) ^ differ by no more than 2-'" x 2"^^ = 2^^^ ^^2) ^YAch is sufficiently 

small for a 12-bit system. (And two more bits would suffice for building a 16-bit 
system.) 

5. 3 Reduction of Multiplier Size in the Fine Stage 

In the fine rotation stage, the computations involved in generating are 

X, = xX\-0:i2)-YA 

/ , \ (5 55) 

Since 16/1 < ^'^'^ it follows that 0, can be expressed as 



1904.0140002 



-76- 



6^ - S, S ... S • • • ^2N/3 ^2Nn+i • • • (5.56) 

where 5 is the sign bit. The A'^/ 3 MSBs do not influence the result. This property 
helps to reduce the size of the multipliers that implement (5.55). Even more 
savings in hardware can be achieved by further reducing multiplier size, with just 
a small loss of accuracy. 

Let [7i] represent the 2N /3 MSBs of J^i as in 

Y, = s,y^... y^^^^ 3/3^,3^, . . . = [i;] + A . (5.57) 



Then we must have 



Ay\< 2 ^^^^ . The error contributed to the product Y^dj by 



^1 

using instead of is 

|i;^,-[};]^J=|a,^,|<2-\ (s.ss) 

Therefore, for A^-bit precision, the multiplication Y^d^ can be accomplished with 
a (2N/ 3) X (2N/3) multiplier. 

This method can be applied to the computation Defining [6 J as 

the 2Ar/ 3 MSBs of 0;, and letting A ^ denote the remaining LSBs, we have 

[0]-s.s... 

^ ^N/2+\ • • • ^2N/3 (5-59) 

^0. = ^2N/3^l'-' (5.60) 
The error in calculating 6,^2 using [6 J instead of 0, is 



[0,]+ A,)' - [9,j\/2« < 2-". (5.61) 

Thus 0;^ can be implemented with an(N/3) ^ (Nf^) multiplier, since the A^/ 3 
MSBs of [0;] are just sign-bits. 
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5. 4 Scaling Multiplier Simplification 



10 



As mentioned in Section 5.2, the scale factor \+ S. ^ , can be applied at 

Icos^ I 

the output of the fine stage. A straightforward implementation would use the fiill 
wordlength of ^^^^^^^ ] the product X= ^|^^^^ ] ), which would require 

a multiplier of size N ^ N. But this multiplier's size can be reduced as follows: 
According to (5,54), 0 s S^^^^^^ < 2•^^^ Moving the factor 1+ ^[^^^^] into the 

fine stage, we have 

X, = jr,(i- ^ /2)(i+ ■^™,|)- '^^„) (5 62) 

= ^. + ^.(^„.,|-^/2)-j;^, (5.63) 
The only significant error in approximating (5.62) by (5.63) is the absence of the 
^j^l^^^^j term in the factor multiplying Y^. But this is tolerable since, according 

to (5.54) and (5.40), 



OA 



< 2- 



(5.64) 



In view of (5.40) we have 0 ^ 0/ < 2'^" which, combined with (5.54), yields 



15 



cos^, I * 



< 2 



-2N/3 



(5.65) 



Thus, if we truncate d^^g - to Ambits, only the least significant N I bits in 



the truncated result will be non-sign bits. Therefore, in our computation of 
Xy{ S^^g - 0/^/2) in (5.63), if we truncate XytoN/3 bits, we can use an (iV/ 3) 

X (N/3) multiplier, with the product's error bound being 



20 



<2-^. 



(5.66) 
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By merging the scale factor of the coarse stage into the fine stage, we thus 
replace multiplications by the scale factor by additions. The final architecture is 
shown in FIG. 39, where the size of the multipliers is shown in FIG. 40. 

FIG. 39 illustrates an angle rotator 3900 according to embodiments of the 
invention. The angle rotator 3900 includes a ROM 3902, a fine adjustment circuit 
3904, a first butterfly circuit 3908, and a second butterfly circuit 3910. The angle 
rotator 3900 rotates an input complex signal 3906 according to angle 6 to 
produce a rotated complex output signal 3912. The angle 6 can be broken down 
into a most significant portion (or word) 6^ and least significant portion (word) 
9^. Note that normalized angle values are shown in FIG. 39, as represented by the 
9 nomenclature. However normalized angle values are not required, as will be 
understood by those skilled in the arts. 

The ROM 3902 stores the following for each corresponding 9 : sin 
[cos ej 5. , , and Gm - 9^, where all of these values have been exactly 

[cos^ ] 

defined in preceding sections. To summarize, the sin Oj and [cos 9i] values are 
MSBs ofsin 9m and cos 9i, respectively. The ^j^^^ j error value represents the 

difference between the cos 9^ and the [cos 9i] value.( The exact definition for 
^1 is given in (5.47)) Likewise, the (9^ - 9 J error value represents the 

difference between sin 9^ and the sin 9^ value. (The exact definition for 9„ is 
given in equation (5 .30)) 

The butterfly circuit 3908 includes multiple multipliers and adders as 
shoAvn. The implementation of these multipliers and adders is well known to those 
skilled in the arts. In embodiments, the sizes of the multipliers and adders are as 
shown in FIG. 40. Note that savings are obtained on the size of the multipliers 
because of bit truncated approximations that are described above. This produces 
a faster and more efficient angle rotator compared to other angle rotator schemes. 

The operation of the angle rotator 3900 is fiirther described in reference 
to the flowchart 4100. As with all flowcharts herein, the order of the steps is not 
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limiting, as one or more steps can be performed simultaneously (or in a different 
order) as will be understood by those skilled in the arts. 

In step 4102, the input complex signal is received. 

In step 4104, the sin 0i, cos [OJ, Sr and 0^ - 6„ values are 

retrieved from the ROM 3902, based on the rotation angle 6. 

In step 4106, the butterfly circuit 3908 multiplies the input complex signal 
by the sin O^ and [cos 6J values to perform a coarse rotation of the input 
complex signal, resuhing in an intermediate complex signal at the output of the 
butterfly circuit 3908. 

In step 4108, the adder 3914 adds the 6l value to the error value - 6„ 
to produce a 0; angle. 

In step 41 10, a fine adjustment circuit 3904 generates a fine adjust value 
( ^ 1 - 6/0 based on the 0, angle and S, , . 

[cos(^] / & [cos6l] 

12^e butterfly circuit 3910 multiplies the intermediate complex 
signal by the 0/ angle, and tHS-ftie^justment value {Sr , - 0/) to perform a 

fine rotation of the intermediate complex number, resuttiagin the output complex 
signal 3912. 

In embodiments, the ROM 3902 storage space is 2^^ words, where A/^ is 
the bit size of the real or imaginary input complex number 3906. Therefore, the 
overall size of the ROM 3902 can be quite small compared with other techniques. 
This occurs because of the two-stage coarse/fine rotation configuration of the 
angle rotator 3900, and saving of storing sin 0^, [cos OJ instead of sin 0 and cos 
0. Also, there is another advantage to having a small ROM: in certain 
technologies it is awkward to implement a ROM. Thus, if only a small ROM is 
needed, it is possible to implement the ROM's input/output relationship by 
combinatorial logic circuits instead of employing a ROM. Such circuits will not 
consume an unreasonable amount of chip area if they need only be equivalent to 
a small ROM. 
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5. 5 Computational Accuracy and Wordlength 

In this section we study the effect of quantization errors on the final 
output's computational accuracy and the most efficient way to quantize the data 
for a given accuracy. 

In our algorithm, the errors can be classified into three categories. The first 
category is the quantization of the values in the ROM table. The second category 
is the error due to the truncation of data before muhiplications, to reduce 
multiplier size. The third type of error is that resulting fi-om approximating sin 6, 
by 6,. Quantization errors are marked in FIG. 40 with an e marker as shown. The 
total error can be obtained by combining the errors propagated fi-om each source. 
To calculate the propagated error at the output with a given error at the source, 
we can first identify all paths by which the error reaches the output and then use 
the approach discussed in Section 5. 1.3. Let us first examine all the error sources 
and determine their effects on X, which is the real component of the output 
complex signal 3912. Table 5. 1 displays this information. (Sinular results apply 
tor.) 

The values stored in the ROM are sin 01, [cos G J, 6^ - Q„ and <5[^^^^ j , 
where sin 0, and [cos 0,] are MSBs of sin G^^and cos Bj, respectively. A loss of 
precision due to ROM quantization error depends only on the number of bits used 
in representing 0^, - 0„ and , . 
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Table 5.1 EfTect of errors at the X output 





Error source 


Error at the 
output 




quantizing (6^ -6„) in ROM 




^2 


quantizing ( ^5^^^^^ ^ ) in ROM 




^3 


truncating Yy for YyQ^ 






truncating for 










^2 


^5 


truncating 0; for 0/ 


Q 


^6 


truncating 0/ for 














r 

quantizing ^z, ~ ^ 






C^ 


quantizing X at the output 






approximating sin 0^ by 0; 




^,0 


neglecting d^^^^^O^m (5.62) 


^1 (-4-<i/'<) 



The total error in A" can be obtained by combining all the terms in the third column 
of Table 5 .1: 



6 Ml 



(5.67) 
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Since in Table 5. 1 is a truncation error, we have ^0. If we quantize 
^[cos^,] rounding it upward before storing it in ROM, then ^0. This way 
such errors tend to cancel each other. Cancelling errors are grouped together in 
(5.67) since the magnitude of their combined error is no greater than the larger of 
the two. This yields seven terms in (5.67), each contributing a maximum possible 
error of 2'^. If the multiplier sizes are as indicated in FIG. 40, the total error in 
JTis bounded by 7 x 2"^. 

From the above analysis it can be seen that the computation errors 
resulting from hardware reduction have similar magnitudes and no particular 
source dominates. This seems to provide the best trade-ofFbetween the hardware 
complexity and the accuracy of the entire system. 

According to (5.67), the total output error can be reduced by increasing 
the internal data wordlength and the wordlength of each sample in the ROM. For 
each bit increase, we get one more bit of precision at the output. Therefore, we 
can design the processor to have the minimum hardware for a given precision 
requirement. Next, we give a simulation example to illustrate this method. 

Example. A cosine waveform with an error less than 2"*^ is specified. 

According to (5.67), we chose 15, as indicated in FIG. 40. We obtained the 

maximum error to be approximately 5 x 10"^, which is considerably smaller than 
2-12 

In FIG. 40, the ROM is shown as having 2 ^ words to achieve no more 
than a total error of 7 x 2"^ in the X output. If is not a multiple of 3, we can 
choose the smallest A^' > AT that is a muhiple of 3. Having 2 ^ words in ROM, 
of course, suffices to achieve the required precision. As discussed before, the total 
output error is a combination of errors from various sources, such as from 
quantizing the data before multiplications and from approximating sin 6; by 0;, etc. 
However, our error bound estimation is rather conservative. Hence, the ROM 
size can be perturbed to determine the minimum size to satisfy a specific precision 
requirement. Our experience in designing the angle-rotation processor has shown 



1904.0140002 



-83- 



N 

that, even by rounding — down to determine the ROM size, the total error is still 
less than 7x2'^. 



5. 6 Comparison with the Single-Stage Mixer 

As mentioned earlier, the main advantage of the two-stage angle rotator 
5 is that it requires only a small ROM 3902. For the single-stage angle rotation, the 

ROM size grows exponentially with the precision of the angle. Thus, our two- 
stage method is well-suited for applications where more than 14 bits in the input 
angle are required. In this case, the sine lookup table for the single-stage angle- 
rotator, even with compression, is too large for high-speed operations (Vankka, 
10 J. , "Methods of mapping from phase to sine amplitude in direct digital synthesis," 

IEEE Trans. Ultrasonics, Ferroelectronics and Freq. Control 526-534 
(1997)). However, the following comparison of our method to a well-known 
single-stage method with 14-bit input angle shows that even in this case our 
method has advantages, and this is true even when the single-stage method is 
1 5 optimized for that particular precision requirement. 

To compare, we use the quadrature direct digital frequency synthesizer/ 
mixer (QDDFSM) with 14-bit input angle and 12-bit input data that is reported 
in (Tan, L. and Samueli, H., "A 200-MHz quadrature frequency synthesizer/mixer 
in 0.8-^m CMOS," /£E£:y, Solid-State Circuits 30. \9^'2mi\99'S)^. It achieves 
20 84.3 dB spurious free dynamic range (SFDR). According to this method, the sine 

and cosine values are generated using a DDFS, which employes lookup tables for 
these values. To reduce the ROM size, ROM compression techniques are used. 
The DDFS is followed by four 12x12 real multiplications. 




^^i^^ ^Fet-oyr^struc^^ chose the internal wordlengths and multiplier sizes 
25 as indicated in FIG. 42?"TTtefljase-accumulator that generates 6, as well as the 

circuit that maps an angle in the rMgS^fO,.^] into [0, ti/4], are described in 
(Madisetti, A., "VLSI architectures and IC imptemontations for bandwidth 
efficient communications," Ph.D. dissertation, University oF'^Salifornia, Los 
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Vj^^^-An gclos ( 1 99 6))-Jliese^tructo^ also employed here in our test. Truncating 
X the 32-bit phase word to 14 bitsTSusltnietnre^ws-achiej^ SFDR of 90.36 dB, 
as shown in FIG. 43. This is 6 dB better than the single-stagemeft 

The integrated circuit that implements this structure is currently being 
built. A preliminary estimation of its hardware complexity yields a similar 
transistor count as that of (Tan, L. and Samueli, H., IEEE J. Solid-State Circuits 
30: 193-200 (1995)). Thus, using approximately the same number of transistors, 
our structure achieves a 6 dB performance gain. Our structure requires a much 
smaller ROM (17 x 25 = 425 bits) in comparison to the single-stage method, 
which needs a 3072-bit ROM when the ROM compression technique is employed. 
Since the ROM access is hard to pipeline, it is usuaUy the bottleneck in the data 
path, thereby limiting the achievable data rate. Hence, one pronounced benefit of 
having a much smaller ROM would be the much faster ROM access. Also, there 
is another advantage to having a small ROM: in certain technologies it is awkward 
to implement a ROM. Thus, if only a small ROM is needed, it is possible to 
implement the ROM's input/output relationship by combinatorial logic circuits 
instead of employing a ROM. Such circuits will not consume an unreasonable 
amount of chip area if they need only be equivalent to a small ROM. 

5.7 A Modified Structure When Only One Output is Needed 



In some applications, such as the implementation of the trigonometric 
interpolator discussed in the previous sections, only one output, say A' is needed. 
In such cases, obviously, we can eUminate certain computations used to generate 
Y. However, using the angle rotator 3900, only those generatmg Y in the fine 
stage are subject to deletion, while the coarse stage must remain the same, since 
we need both^j and Y^ to generate the X output. Let us seek to fiirther simplify 
the coarse stage by attempting to eliminate one multiplication by cos 
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5. 7. 1 Modifications to the Coarse Stage 



10 



15 



If we factor out the cos term of the coarse stage in (5 .16), we can 
then apply the factor cos 6^ to the output of the second stage in (5. 1 7), because 
the two operations (scaling and rotation) are permutable, to obtain 



(5.68) 




7, = 



r 



X = cos (9 



^1 



V 

f f 



2 ) 



J 



Y = cos0^ 



M 



V 



2 J 



(5.69) 



In this case, we have only two multipliers in the coarse stage (5.68), and the 
multiplications by the scale factor Oj^ are applied to the output of the fine stage 
(5.69). Unlike the situation in (5.16) and (5. 17), if only one output fi-om the angle 
rotator, say ^, is needed, we can also eliminate one more multiplier — the one that 
multiplies the coarse stage output with the cos dj^ factor. As in Section 5.2, we 
now seek to simplify the coarse stage in (5.68). 

Let tan e„ be tan 0^^^ rounded upward at the (N/syrd bit. In other 
words, writing 6^^^ as the binary number 



tan 0^ = 0.^1 *iv/3* 



N/3 



where b„ € {0, 1 }, tan 0^ is obtained from tan 0^ according to 
tan0„ = 0.6,. 6^,3 + 2-^'^ 

Obviously, 



0< tane„- tanG^ < 2-^'^ 



(5.70) 
(5.71) 
(5.72) 



20 



The M3-bit number tan6„ decreases the number of partial products needed in 
computing JTotane^ and Kotan6„ to at most a third of those needed for A'otanOj^ 
and 7otan0;i^. 
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The resulting fine-stage angle is 0/ = 9-6„. Thus, as in Section 5.2, a 
modified fine-stage angle compensates for a simplified coarse-stage angle. If 0; 
satisfies (5.40), we then have |sin0/-0J<2"^/6. That is, the approximations 

sin0; = 0; and cosG ^ = 1 - 0 f / 2 can be applied. The proof that (5 .40) holds is 
as follows: 

Proof, According to the mean value theorem 

tane^- tan0 



9m - 



^ = tan' ^ (5.73) 



where ^=0a/ + (0„ - 0^/)^, 0 ^ h <, 1 . The derivative tan'^ satisfies 

tan' ^ = 1 + (tan^ f > 1, for every ^ . (5.74) 

10 Re-arranging (5.73), and using (5.74), we have 

tane„-tan0^ 

0«-eA/ = tan'^ ^ tan9,-tan9^. (5.75) 

Hence, according to (5.72), 

0< e^- 0^ < 2-^^^ (5.76) 

By definition, 

15 O^0^<2-^^. (5.77) 

Therefore, subtracting (5.76) firom (5.77) yields 

-2-^'^<0,-(0„-0^)<2-^'^ (5.78) 

which is exactly (5.40) because 

0;=0-0„ = 0^+0^.0^ = 0^.(0„-0^). (5.79) 

20 This concludes our proof 

This indicates that, instead of storing the tanO^^^ values in the ROM , we may store 
tan0„, which has N/3 bits for each sample, and we may store 0„ - dj^. This results 
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in a reduction of the multiplier size in the coarse stage. The difference between 
e„ and Qj^ can be compensated in the following fine rotation stage. Furthermore, 
the approximations (5.6) and (5.7) still apply to 0,, in view of (5.40). 
We can now implement the coarse rotation stage as follows: 

y, = 7o + Xotaii0„. (5^°) 

Accordingly, the scale-factor at the output of the fine stage is cos6„ instead of 
cosQj^. Since 0; satisfies (5.40), the fine stage simplification is similar to the 
method described in Section 5.3. Next we examine how the multiplications of the 
fine-stage output by cos0^ can be simplified. 

5. 7. 2 Scaling Multiplier Simplification 

A straightforward implementation would use the fiiU wordlength of cos0^ 
in the product X = JfjCOsO^, which would require a multiplier of size NxN. But 
this multiplier's size can be reduced as follows: By defining [cos0„] as the 2M3+1 
MSBs of cos0„ the scale factor can be written as 



cosG^ = [cos0^]+ A^3e. = [cosG^] 1 + 



[cose„]y 



(5.81) 



A. 



Let us define S^^^^ = ^ """^^ ^ and, since 0 ^ 0„ ^/4, we surely have [cos 0J > 
- [cos^„] 



0.5, and hence 



0<<J,„3^^ <2x2-'^''-V (5.82) 



Moving the factor 1 + ^cosc^ stage, we have 

X, = 0f /2)(1+ YA{1+ <Jeos.J (5 83) 

= X, + X,{S^,, - ef 1 2) - YA- (5.84) 



1904.0140002 



-88- 



The only significant ejTorin^proximating (5.83) by (5.84) is the absence ofthe 9,5 

term in the factor multiplying 7,. But this is tolerable since, according to (5.40) 
and (5 .82), 

l^/^co..J< 2-^. (5.85) 
In view of (5.40) we have 0 < Of < 2~^^'^ which, combined with (5.82), 

yields 

l^a>s.„-^/2|<2-^" (5.86) 
Thus, if we truncate S^^^g^ - Of / 2 to Ambits, only the least significant M3 bits 

in the truncated result will be non-sign bits. Therefore, in our computation of 
^\ (^cos/?„ ~ ^ / 2) in (5.84), if we truncate to N/3 bits, we can use an (M3) x 

(A73) multiplier, with the product's error bound being 

\S^^,^-efl2\2-^^'<2~r (5.87) 

The factorization of cosG^ in (5.81) allows a reduction ofthe multiplier to 
approximately 2/3 its original size. In this case, the values of [cosO^] and S^^^ 

are stored in the ROM instead of cos0„. 

The modified structure for one output is illustrated as angle rotator 4400 
in Figure 44.The angle rotator 4400 includes a ROM 4402, a fine adjustment 
circuit 4404, a first butterfly circuit 4408, and a second butterfly circuit 4410. The 
angle rotator 4400 rotates an input complex signal 4406 according to angle 6 to 
produce a rotated complex output signal 4412. As with the rotator 3900, the angle 
0 can be broken down into a most significant portion (or word) 6^^^ and least 
significant portion (word) 0^. Note that normalized angle values are shown in 
FIGs. 39, 40, 42, and 44, as represented by the 9 nomenclature. However 

normalized angle values are not required, as will be understood by those skilled 
in the arts. 
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The ROM 4402 stores the following for each corresponding normalized 
9: tan [cos 8„] S^^^^^ , and - e„, where all of these values have been 

exactly defined in preceding sections. 

In the butterfly circuit 4410, the arithmetic units that are encircled by the 
5 line 44 1 8 can be eliminated when only the X output is needed in the output signal 

44 1 2. This may be desirable for applications where only one output from the angle 
rotator 4400 is needed, such as when implementing a trigonometric interpolator, 
such as interpolator 1000 in FIG. 10 or interpolator 1400 in FIG. 14. 

The operation of the angle rotator 4400 is further described in reference 
10 to the flowchart 4500 in FIG. 45. As vnth all flowcharts herein, the order of the 

steps is not limiting, as one or more steps can be performed simultaneously (or in 
a different order) as vnll be understood by those skilled in the arts. 

In step 4502, the input complex signal 4406 is received. 

In step 4504, the tan 0„ [cos 0„], ^5^^^^^ , and 6^ - d„ values are 

1 5 retrieved from the ROM 4402, based on the rotation angle 0 (or the normalized 

value 0 ). 

In step 4506, the butterfly circuit 4408 multiplies the input complex signal 
4406 by tan 0^ to perform a coarse rotation of the input complex number, 
resulting in an intermediate complex signal at the output of the butterfly circuit 
20 4408. 

In step 4508, the adder 4414 adds the 0^ value to the error value 0^ - 0„ 
V to produce a 0/ angle. 

lHs§tep^ 10, a fine adjustment circuit 4404 generates a fine adjust value 
( ^cose„ - ^f) ba§ed^ the 0; angle and S^^^^^ . 

25 In step 45 1 2, the butt^rflxdrcuit 44 1 0 niultiplies the intermediate complex 

signal by the 0/ angle, and the fine ad^itstqient value ( S^^^^^ - Of) to perform a 

fine rotation of the intermediate complex signal, resulting in the output complex 
signal. 
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In step 4514, the X value for the output complex signal is scaled by the 
[cos 0„] value, resulting in the output complex number 44 12. As discussed above, 
the elements inside the outline 44 1 8 can be eliminated if only the X value of signal 
4412 is desired. Alternatively, similar elements could be eliminated from the 
butterfly circuit 4410 if only the Y value of signal 4412 was desired. 

5. 8 Application of Angle Rotation Processors 

This subsection describes exemplary applications for angle rotator 
processors. These applications are provided for example purposes only and are not 
meant to be limiting, as those skilled in the arts will recognize other applications 
based on the discussions given herein. These other applications are within the 
scope and spirit of the present invention. 

One application for the angle rotation processor is the Quadrature Direct 
Digital Frequency Synthesizer/Mixer (QDDFSM), including a few special cases 
that are candidates for the angle rotator algorithm. One is the case when only one 
of the outputs (XoT Y) is desired, as shown by angle rotator 4400 (FIG. 44). As 
shown in FIG. 44, this is accomplished by simply deleting the hardware required 
for the computation of the unused output. Yet another special case of QDDFSM 
is the Direct Digital Frequency Synthesizer (DDFS). In DDFS configuration we 
simply fix the input vector (Xq, Yq) to be (1,0). This enables the complete 
elimination of the coarse stage by taking advantage of the fact that 1 X A = A and 
0 X A = 0. In the following section we will concentrate our discussion on the 
QDDFSM, since it is the general case, while keeping in mind the special cases and 
the associated hardware reductions mentioned above. 

5. 8. 1 Using the Angle Rotation Processor in a Quadrature 
Direct Digital Frequency Synthesizer/Mixer 

The frequency synthesis and mixing operation can be described with the 
following pair of equations, which relate an input vnthx-y coordinates (Xq, Yq) and 
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a frequency control word (few) for the synthesizer, to an output with new x-y 
coordinates (X, Y). The following pair of equations establishes the relationship 
between (Xo, YqX few, and (X, Y). 

X=[XqX cos(fcw X w)] - X sin(fcw x «)] 

(5.88) 

Y=[YqX cos(fcw X w)] + [Xq X sin(fcw x n)] 
where n is the time index 



Per (5 .88), since the sine and cosine functions are periodic with period 2n 
(i.e., few X w = <fcw X w >2„ = <|), where o is a modulo operator) an overflowing 
adder is used as a phase accumulator to compute ^ from the input few, as shown 
by the the phase accumulator 4600 in FIG. 46. 

Now, for any given time instance w, we have a corresponding angle ^ from 
the phase accumulator, hence the original pair of equations (5.88) for QDDFSM 
can be rewritten in terms of the angle <t) as follows. 

X=[XoXcos<^] - [7oXsin(|)] 

(5.89) 

Y=[YqXcos^] + [A^oXsincI)] 

Note that the expressions (5.89) are exactly those of an angle rotator 
expressed by equations (5. 1). By applying a phase accumulator fed by an few, we 
have converted the QDDFSM into an angle rotation application. The only conflict 
between the above expressions and the angle rotation processor is that the angle 
rotation processor takes an angle 0 in the range [0, 7c/4], while the angle ^ in the 
above expressions is in the interval [0,27i;). 
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5. 8.1.1 A General Angle Rotator for Arbitrary Input 
Angles 

Let us consider the changes necessary to make the angle rotation processor 
use an input angle <|> that may lie outside the [0, n/4) range. Fortunately, a simple 
interchange operation at the input of the coarse stage, and an interchange/negate 
operation at the output of the fine stage is all we need in order to map (j) into an 
angle 9 in the range [0, n/4] and use it as the input to the angle rotator. Even 
though the input angle 0 is in the range [0, n/4], the rotation by 9 along with the 
interchange and interchange/negate operations make the overall rotation of the 
input (Xq, Yq) equivalent to a rotation by the original angle <|) in the full range 
[0,2n), The latter is possible because of the convenient symmetry properties of 
sine and cosine functions over the range [0, In]. 

For example, sin (|> = -sin(<|) - ti;) and cos <|) = -cos(<|) - 7c), while sin ^ = 
cos(<|) - 7c/2) and cos <|) = -sin((t) - n/l), and finally, for 0 < (|) < n/4, if we write 
n/4 + (|) for (j) then sin(7c/4 + <!))= cos(7i;/4 - (j)) and cos(7c/4 + <|)) = sin(7c/4 - (j)). 
Using the first pair of trigonometric identities, we can map 4) into the range [0, n]) 
by simply performing a negate operation at the output of the angle rotator. Using 
the second pair of identities along with the first pair enables one to map (|) into 
the range [0, n/2) by performing negate and interchange operations at the output 
of the angle rotator. Finally, using all three pairs of identities, the angle ^ can be 
mapped into the range [0, n/4) by performing an interchange operation at the 
input of the angle rotator, along v^th interchange and negate operations at the 
output of the angle rotator. Note that all of these interchange and negate 
operations are conditioned only on the issue of which octant <|) is in. This means 
that if (f) is a normalized angle, then the interchange and negate decisions depend 
only on the top three MSB bits of <|). The following tables show the interchange 
and negate operations required for all eight octants (specified by the three MSB 
bits of (|)). It is evident, as well, that other interchange and negate criteria for the 
input and output would also be suitable. 
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This table indicates when an interchange operation is required at the input 
and when an interchange operation is required at the output of the angle rotator. 



Octant of <t> (3 MSBs of i>) 


Output Interchange 


Input Interchange 


1-st octant (000) 






2-nd octant (001) 




Interchange inputs 


3-rd octant (010) 


Interchange outputs 




4-th octant (Oil) 


Interchange outputs 


Interchange inputs 


5-th octant (100) 






6-th octant (101) 




Interchange inputs 


7-th octant (110) 


Interchange outputs 




8-th octant (111) 


Interchange outputs 


Interchange inputs 



The following table indicates when a negation operation is required at the 
output of the angle rotator. 



Octant of (J) (3 MSBs of 6) 


Negation of output X 


Negation of output Y 


1-st octant (000) 






2-nd octant (001) 


Negate output X 




3-rd octant (010) 




Negate output Y 


4-th octant (Oil) 


Negate output X 


Negate output Y 


5-th octant (100) 


Negate output X 


Negate output Y 


6-th octant (101) 






7-th octant (110) 


Negate output X 




8-th octant (111) 







Note that the flag for input interchange is simply the 3rd MSB bit of <j), while the 
flag for output interchange is just the 2nd MSB bit of <|). Finally, to produce the 
remapped angle 0 in the range [0, 7i/4) for the angle rotation processor, we 
simply take the remaining bits of (|> after stripping the top two MSBs and 
performing a conditional subtract operation to produce 0. More specifically, if the 
MSB bit (after stripping the two MSB bits) is low, i.e., the angle is in an even 
octant (numbering them 0, . . ., 7), we pass the angle unchanged, otherwise we 
perform a "two's-complement type" inversion of the angle. Note here that after 
such remapping operation, the MSB bit of 0 is set to one only in the case when 0 
= 7t/4. This fact is usefial in determining the required amount of lookup table in the 
angle rotation processor. In other words, even though the MSB bit of 0 is an 
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address to the lookup table, since we know that when it is ' 1' the remaining bits 
have to all be '0' we only need to allocate a single address for that case (as 
opposed to increasing the size of the lookup table by an entire factor of two). 



5. 8, L 2 Adapting the General Angle Rotator to Make a 
QDDFSM 



The structure of the QDDFSM using an angle rotation processor 3900 is 
depicted in the figure 4700. It simply requires the employment of a phase 
accumulator 4702 and a conditional subtract 4704 to provide an input angle fi"om 
the input fi"equency control word few. We refer to the system of Fig. 47 with the 
10 phase accumulator excluded as a General Angle Rotator. It has the capability to 

receive an angle in the interval [0, 2n) and to perform an angle rotation of the 
input data {X^, by that angle. We show a general angle rotator in Fig. 48, but 
one in which further structural simplification has been made. The method of 
performing these simplifications will be discussed next. 



15 5,8.2 How to Use the Conditionally Negating Multipliers in 

the General Angle Rotator 



For a moment assume we have a powerful technique for making 
conditionally negating multipliers. What we mean by that is a multiplier which 
takes a negate flag to produce an output depending on that negate flag as follows: 

20 The output is simply the product of the input signals if the flag is low (0) and the 

output is the negative of the product of the input signals if the flag is high (1). 

Each one of the two outputs in the coarse and fine stages is computed with 
two multipliers and one adder as shown in Fig. 47. These multipUers and the 
adder are implemented in a single Carry-Save Adder (CSA) tree, with the partial 

25 products being generated fi-om Booth decode modules corresponding to the two 

multipliers. This technique of employing a single tree eliminates the need for 
intermediate carry propagation fi-om each multiplier and makes the propagation 
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delay of each rotation stage very short. Note that the single CSA tree 
implementation is possible since the multipliers are operating in parallel. 
Furthermore, because the structure that is needed to compute one output of a 
rotation stage is identical to the structure required by the other output (with the 
exception of the minus sign), a single CSA tree can easily be interleaved between 
the two outputs for a significant amount of hardware savings. The minus sign at 
the output of the multiplier can be implemented very efficiently by the technique 
described in the following sections (using the conditionally negating multiplier). 
The negation or non-negation of the multiplier output can be controlled with a flag 
that changes between the two cycles of the interleave operation. 

The angle at the output of the conditional subtract module 4704 in Fig. 47 
is in the range [0, n/4]. As already discussed, the outputs for the angles outside 
this range are constructed by mapping the angle into the range [0, n/4] while 
conditionally interchanging the inputs (inputs to the coarse stage) and 
conditionally interchanging and negating the outputs (outputs of the fine stage) of 
the angle rotator. A negation at the output of the fine stage simply means changing 
the output signs of the multipliers and negating the input of the adder coming fi*om 
the input of the fine rotation stage. Changing the output signs of the multipliers 
is once again accomplished by using conditionally negating multipliers. The 
negation of the input to the fine rotation stage can easily be implemented with 
XOR gates and a conditional high or low bit insertion into the CSA tree at the 
position corresponding to the LSB location of the input. Since this conditional 
high or low bit is inserted in the CSA tree, there is no additional carry propagation 
introduced for the negation of the input. Note that the latter technique eliminates 
any circuitry required to implement the conditional negation of the outputs, and 
hence eliminates any carry propagations associated with two's complement 
numbers. 

Furthermore, the conditional interchange of the outputs can be 
implemented by conditionally interchanging the inputs of the fine rotation stage 
and appropriately controlling the signs of the multiplier outputs in the fine stage, 
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The conditional interchange of the fine stage inputs can be propagated to the 
inputs of the coarse stage with the same line of reasoning. Remember that the 
inputs to the coarse stage were conditionally interchanged according to the three 
MSBs of the input angle anyway. In conclusion, the conditional interchange and 
negation operations of the outputs can be implemented by modifying only the 
condition of the interchange at the inputs of the coarse stage and appropriately 
controlling the multiplier output signs by using conditionally negating multipUers 
(which we had to do for interleaving anyway). This eliminates the conditional 
negate and interchange block at the output of the fine stage entirely (i.e., it 
eliminates muxes and two's complement negators), and also eliminates the need 
for storing and pipelining control signals (i.e., it eliminates registers) to perform 
the conditional interchange and negation operations at the output. The resulting 
General Angle Rotator 4800 is now depicted in the following FIG. 48. 

5.8.2.1 Booth Multiplier 

There are many algorithms for digital multiplication. One of the most 
popular is the Booth multiplier. The essence of the Booth muhiplier is in the 
decoding scheme performed on the multiplicand to reduce the number of partial 
products which, when added together, produce the desired product. For an JVx M 
Booth muhiplier, where N is the wordlength of the multiplier, and M is the 
wordlength of the multiplicand, there will be ceiling(A'^/2) Booth decoders. Each 
Booth decoder will take three bits fi-om the multiplier (with one bit overlapping 
the decoders on both sides) and will manipulate the multiplicand SLCCordingto the 
Booth decoding table 5000 shovra in Fig. 50. Some relevant details for a 10 x M 
Booth muhiplier are depicted in Fig. 49, especially how the muhiplier bits feed 
into the Booth decoders to produce the five partial products which, when added, 
compute the resuk (the product of the multiplier and the multiplicand). 



1904.0140002 



-97- 



5. 8. 2. 2 Haw to Make a Negating Booth Multiplier 

Suppose we wish to make a multiplier that produces the negative of the 
product. More specifically, suppose we wish to multiply two signals and A/ and 
get -C = -(NxM). The latter can be accomplished in a number of different ways. 
The most obvious is perhaps to use a regular multiplier to produce the product C 
= (NxM) and then negate C to achieve -C = -(N x M). In case of two's 
complement representation, this approach requires an additional carry propagation 
chain through the negator, which is costly in terms of speed and additional 
hardware associated with a negating circuit. Another approach, described below, 
is more favorable in a few key aspects. 

The product C is essentially the result of adding a number of partial 
products, which are generated by the Booth decode blocks as described in the 
previous section. Therefore, we can write the following sum expression for C: 

n 

C=Y.Pi (5 90) 

7 = 1 

where /7, are the n (in the 10 xM example above ai = 5) partial products generated 
from the n Booth decoders. Note that, in order to negate C, we can negate all of 
the partial products and proceed with the summation of the negated partial 
products to produce -C. The expression for -C is then the following: 

n 

-C=Z-A > (5 91) 

1=1 

where are the negated n partial products generated from the n Booth decoders. 
Let us investigate how the Booth decoders need to change to produce the desired 
negated partial products. All we need to do is to change the decoding table 5000 
fi-om that of Fig. 50, to the decoding table 5100 in Fig. 51. Note that the 
difference between the tables is only in the partial product columns and, more 
specifically, the partial product colunrn 5102 of table 5100 is the negative of the 
partial product column 5002 of table 5000. This means that by simply modifying 
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the Booth decode table to the negating Booth decode table shown in Fig. 5 1, the 
result will be the negative of the product, as desired, with absolutely no additional 
hardware and absolutely no speed penalty. An example for a 10 xM negating 
Booth multiplier 5200 is shown in Fig. 52. 



A particularly interesting case arises when one wishes the multiplier 
product to be negated sometimes, and normal (non-negated) the other times. One 
can extend the idea presented in the previous section to craft the following 
powerful technique. Let us investigate the original Booth decode table 5000 
depicted in Fig. 50 and the negating Booth decode table 5100 of Fig. 51 a bit 
more closely. Note the horizontal line of symmetry that runs through the midline 
of both decoding tables. This line of symmetry suggests that we can create the 
negating Booth decode table 5100 from the original Booth decode table 5000 by 
simply inverting the three bits (b2 bl bO). For example, if the three bits (b2 bl bO) 
are (0 1 0), then, according to the original Booth decode table, the corresponding 
partial product is A, where A is the multiplicand. If we invert the three bits (b2 
bl bO) as suggested above, we will have (10 1) and the corresponding partial 
product will be -A, exactly what is needed for a negated partial product. 

Given a signal F which specifies when the output of the multiplier should 
be negated and when not (F = 0 implies regular multiplication, F = 1 implies 
negating muhiplication), F can simply be XORd with the three bits (b2 bl bO) at 
the input of the regular Booth decoders to make a new conditionally negating 
Booth decoder, hence a conditionally negating multiplier. The details of a 
conditionally negating Booth decoder 5300 are captured in Fig. 53 . Note that with 
a minimal amount of hardware {N XOR gates for an A^x A/ multiplier, which is 
insignificant compared to the hardware cost of the entire multiplier), we have the 
means to control the sign of the multiplier product. Also note that the overall 
latency of the multiplier is increased insignificantly since the latency through a 
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single XOR gate is much smaller than the latency through the entire multiplier. 
Furthermore, the latency of a single XOR gate is much smaller than the latency 
associated with a carry propagation chain that would be necessary if one built such 
a circuit with a two's complement negator. A 10 xM conditionally negating 
multiplier 5400 is shown in Fig. 54. 

5. 8.3 Using the Angle Rotation Processor in a Quadrature 
Direct Digital Frequency Synthesizer 

As mentioned above, the angle rotator is useful in implementing various 
forms of direct digital synthesizers. In this case, all starting points for the angle 
rotations are = 1, ^ 0 (with, of course, the various usual 
interchange/negation requirements). Fig. 55 shows a quadrature direct digital 
synthesizer (QDDS) 5500, a system having two outputs, one being samples of a 
cosine waveform and the other being samples of a sine waveform. An exact 90- 
degree phase offset between the two waveforms is obtained by the QDDS, and 
numerous applications for such a device are well knovm. No and input 
samples are shown in the Fig. 55 system. These fixed values have been "built in" 
and used to greatly simplify the coarse rotation stage. 

Notice that the angle rotator 5502 is preceded by a system 5504 that 
generates a data stream of input rotation angles, a so-called overflov^ng phase 
accumulator 5506, and its input is a single fixed data word that precisely controls 
the frequency of the output sample waveforms. The three MSBs of each phase 
accumulator output word, of course, assess the approximate size of the angle that 
is being used as a rotation angle (i.e., these three bits show how many octants the 
rotation angle spans), and they are stripped off to control the interchange/negation 
operations that are appropriate for obtaining the desired output samples. Also, the 
third MSB is used, as described previously, to determine whether or not to 
perform a "two's complement type" inversion of the LSBs. One other operation 
is required by the "Conditional Subtract" module 5508 shown in Fig. 55; in 
addition to stripping off the three MSBs, it appends one MSB having the value 
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zero except in the case where a rotation angle of exactly 7c/4 is required. In that 
case, the appended MSB is one and all other ROM-address bits are zero. 

A special case of the QDDS system, one having only a single output data 
stream, which could be either of the two, but which we call the "cosine-only" case, 
5 is also useful for various well-known applications. Fig. 56 and Fig. 57 show two 

specializations of the angle rotator circuits previously discussed to implement the 
cosine-only DDS. The system 5600 in FIG. 56 results from specializing the 
angle-rotation system 3900 in FIG. 39. The system 5700 in Fig, 57 is a 
specialization of the angle rotator 4400 in FIG. 44. 

10 5.9 Conclusion 



Based on the design method discussed, for a given accuracy requirement, 
an architecture with the least amount of hardware is produced by balancing the 
precision of intermediate computations and the complexity of each arithmetic 
block, while keeping the output error within the specified bound. Furthermore, our 

15 architecture consolidates all operations into a small number of reduced-size 

multipliers. This permits us to take advantage of many efficient techniques that 
have been developed for multiplier implementation, such as Booth encoding, 
thereby yielding a smaller and faster circuit than those previously proposed. 

Simulations and preliminary complexity estimation show that, even 

20 comparing to the method of (Tan, L. and Samueli, H., IEEE J. Solid-State 

Circuits 30: 193-200 (1995)) that is optimized for a 14-bit input angle, our method 
achieved 6 dB more SFDR while using approximately the same number of 
transistors as those needed by (Tan, L. and Samueli, H., IEEE J, Solid-State 
Circuits 30: 193-200 (1995)). In addition, since our structure employs only a small 

25 ROM, it overcomes the problem of slow access time that occurs when large 

ROMs are used, thereby facilitating a higher data rate. Using the two-stage 
method, when a higher precision is needed, it is very straightforward to satisfy 
such a requirement, since more accurate results can be attained simply by 
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increasing the wordlength and the multiplier size. For the single-stage method, 
however, when high precision is desired, the required lookup table is likely to be 
too large to be practical, particularly for high-speed operation. 

6. Symbol Synchronization for Bursty Transmissions 

We have thus far discussed methods that provide efficient implementations 
of the resampler for symbol synchronization in a digital receiver using 
trigonometric interpolation as well as the phase rotator for carrier recovery. To 
produce the correct samples, a timing recovery circuit must supply the resampler 
with symbol timing information, as shown in Figure ID . We will now consider 
how this can be accomplished. 

6. 1 Initial Parameter Estimations for Burst Modems 

There are many methods to derive timing information from the received 
signal. According to their topologies, synchronization circuits can be divided into 
two categories: there are feedback and feedforward schemes. Feedback structures 
usually have very good tracking performance, and they work quite well in 
continuous mode transmissions. For packet data systems used by third-generation 
mobile communications, where the transmission is bursty, it is essential to acquire 
initial synchronization parameters rapidly from the observation of a short signal- 
segment. 

A typical packet format is shown in Figure 58. It includes a short preamble 
5802 followed by user data 5804. The preamble 5802 is a set of known 
modulation symbols added to the user data packet at the transmitter with the 
intention of assisting the receiver in acquisition. 

There are many approaches to burst demodulation, depending on the 
specific system requirements. In one approach (S. Gardner, "Burst modem design 
techniques, part 1," Electron. Eng. 77:85-92 (Sept. 1999); Gardner, S., "Burst 
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modem design techniques, part 2," Electron, Eng. 77:75-83 (Dec. 1999)) , the 
receiver first detects the presence of the preamble, using a correlator, whose 
output should produce a large magnitude when the preamble is present. It then 
estimates the symbol timing. If the sampling fi-equency error is small, the total 
change of the timing phase fi*om the start of the short preamble to the end is 
negligible. Next, it estimates the initial carrier frequency and phase. The above 
steps assume that the impairment caused by the channel is small enough that the 
modem can successfully track the timing carrier phase prior to equalization. 
Otherwise, equalizer training prior to the timing and carrier recovery is needed. 

With a typical preamble of 8 to 32 symbols, depending on the required 
system performance, for QPSK modulation, rapid acquisition is desired. 
Feedforward timing estimation is known to have rapid acquisition, since it 
produces a one-shot estimate instead of tracking the initial timing through a 
feedback loop. 

A well-known method, digital square timing recovery (Oerder M,, and 
Meyr, H., IEEE Trans, Comm. 55:605-612 (1988)), has shown rapid acquisition, 
but it requires oversampling of the signal at, typically, four times the symbol rate, 
which imposes a demand for higher processing speed on the subsequent digital 
operations. Moreover, it does not work well for signals employing small excess 
bandwidth. However, pulses with small excess bandwidth are of interest for 
bandwidth-efficient modulation. 

For applications where low power and low complexity are the major 
requirements, such as in personal communications, it is desirable to sample the 
signal at the lowest possible rate and to have the synchronizer be as simple as 
possible. In this section, a synchronizer is proposed that needs just two samples 
per symbol period. In addition, it has been shown to work well for small excess 
bandwidth, which is important for spectral efficiency. Using this method, the 
estimations of the symbol timing and the carrier phase can be carried out 
independently of each other. Hence, they can be carried out in parallel. Using the 
proposed structure, the timing and carrier-phase estimators can be implemented 
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efficiently by means of direct computation (instead of a search, as is employed, for 
example, by (Sabel, L., and Cowley, W., "A recursive algorithm for the estimation 
of symbol timing in PSK burst modems," in Proc. Globecom 1992, vol 1 (1 992), 
pp. 360-364) using an efficient rectangular-to-polar converter (to be discussed in 
Section 7). This yields a very small computation load. Thus, this structure is well 
suited for low-power, low-complexity and high-data-rate applications, such as 
those in multimedia mobile communications. 

6, 2 Background Information 

The system model 5900 used in developing the symbol timing and carrier 
phase recovery algorithm described in this section is shown in Figure 59. 

Here h(t) \s?i real-valued, unit-energy square-root Nyquist pulse and w(t) 
is complex white Gaussian noise with independent real and imaginary components, 
each having power spectral density NJ2. 

As mentioned in Section 6.1, a typical data packet for a burst modem 
consists of a short preamble 5802followed by user data 5804. According to the 
approach of (Gardner, S., Electron, Eng. 77:85-92 (Sept. 1999)), the matched 
filter output is sampled every 7; = 772 seconds, i.e., at twice the symbol rate. The 
receiver then detects the presence of the preamble in the received signal by 
correlating the conjugate of the known preamble sequence whose length is L, 
with the sampled data x(nTJ as 

L-\ 

r^(n) = Z ^>inT, + 2mTJ. (6.1) 

m=0 

The correlator output r^(n) should produce a large magnitude \r^(n)\ when the 
preamble is encountered. It then estimates the initial synchronization parameters, 
namely the symbol timing and the carrier phase, assuming the transmitter/receiver 
fi-equency mismatches are insignificant. 

The complex envelope x(t) of the received signal, after the matched 

filter, is 
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^(0 = S ^kg(t -kT-r)^ v(0 (6.2) 

>t=-oo 

where {a^) is a sequence of independent equally-probable symbols with 
E[\a^ 1^ ] = 1 . We also have that v(t) = w(t)(» h(-t) and that g(t) = hft)t» hf-t) 

is a Nyquist pulse. The time delay x and the carrier phase 0 are both unknown. 
5 To estimate the data sequence we want sample values of x(t) at 

t = mT-\-T^ with m an integer, whereas only the samples x(nTJ are available after 
sampling x(t) by a fixed clock. 

Now let us examine how the correlator output relates to symbol timing and 
carrier phase. Inserting (6.2) into (6. 1) yields 

L-\ 00 

r^iri) = X S^^^a^C^T; + 2mT, - kT - r)eJ' 
10 (6.3) 

#n=0 

Since the data are independent, and they are independent of the noise, we have 

fl k=m 
E{a„aA=\^ , (6.4) 

E[alv{nT^ + 2mT,)-\=Q. (6.5) 
According to (6.4) and (6.5), and because T = 2T^ the expectation of r^(n) with 
1 5 respect to the data and the noise is (for simplicity, we omit the constant real scale 

factor L) 

E[r^{ny\=eJ'ginT^-T). (6.6) 
Thus, the mean value of the complex preamble correlator output actually 
equals the sample of the delayed signaling pulse g(t), with delay being x, rotated 
20 by the angle 0. This is shown in Figure 60 for 0 = 0, where g(t) is a raised cosine 

Nyquist pulse with a = 0.35. The total timing delay r can be expressed as 

t=«,7;+n (6.7) 
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where the integer Wq represents the portion of x that corresponds to an integer 
multiple of the sampling interval 7; and 0 ^ |i < is the sampling time mismatch. 

Most practical signaling pulses g{t) are symmetrical and their peak value 
occurs at g(0). If 6 is known, using these properties, v/e can estimate the 
sampling time mismatch from the correlator output r^(n) . In the next section 
we will discuss such an algorithm. We will derive this algorithm by first assuming 
that 0 = 0, Then we will discuss how the method can be carried out independently 
of the carrier phase. Simultaneously, we also derive a phase estimation algorithm 
that is independent of the symbol timing. 

6.3 Symbol Timing Estimation Assuming 0=0 
From (6.6), with 6 = 0, we have 



According to (6.7) and (6.8), if the transmission delay t is exactly an integer 
multiple of we must have ^ = 0, and thus r^(nf) must correspond to the peak 
g(0). Otherwise, we have ^i^O, v^th r^(nQ) and rJ[no^l) being the two correlator 
output values nearest the peak value g<0), as shown in Figure 60. That is, rj[n() 
and r^{no+l) must be the two largest correlator outputs. Therefore, once the 
largest correlator output is located, we can obtain n^, the integer part of x. 

We now turn to finding n. Without loss of generality, let us assume T^=l. 
Replacing n by n^+n we have, according to (6.8) and (6.7), 



For simplicity in our discussion on finding the fractional delay fi, we henceforth 
drop the index /Iq, which corresponds to an integer multiple of sample delays, from 
our notation. Next we define R(^'^) as the Fourier transform of r^(n) \ 



(6.8) 



(6.9) 



00 




(6.10) 
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The expectation of R(e'") can be expressed as 

E[R{e^'^)\ = EiFTir^in))] = FTiE[r^m). (6.11) 
According to (6.9), and (C.4) in Appendix C, we have 

EiRie^")^ = FT(g(n-M)) = e^'"" G(eJ'" ) (6.12) 
where G(e"^ is the Fourier transform of g(ji). Since g(n) is symmetrical, G(e"^ 
must have zero phase. Thus, according to (6. 12), 

arg(E[R(eJ'")]) = arg(e^'"^G(e^" )) = com. (613) 
Evaluating (6. 13) at co = tz/2, we can obtain an estimate of \i as 

2 

//= -arg(i?(e^''2)). (6.14) 

Therefore, the unknown sampling mismatch can be obtained by taking the 
Fourier transform of r^(n) and evaluating the phase corresponding to co = n/2. 

To make the implementation of (6.14) realistic, we should truncate the 
sequence r^(n) before taking its Fourier transform. For example, using only the 
four samples r^(-l), r^(0), /"^(l), and r^(2), we have 

^r(^^"'') = [/-^(O)- r^(2)]^j[r^(~l)-r^(l)l (615) 
Using the correlator output, the ^ value can be obtained by first computing 

Rj-ie^"'^) according to (6. 15), and then from the following: 

2 

Ai=-aig(R^ie''"^)). (6.16) 
n 

For low precision requirements, this operation can be accomplished using a small 
CORDIC processor (Chen, A., et al, "Modified CORDIC demodulator 
implementation for digital IF-sampled receiver," in Proc. Globecom 1995, vol. 2 
(Nov. 1995), pp. 1450-1454) or a ROM lookup table (Boutin, N., IEEE Trans. 
Consumer Electron. CE-38:5-9 (1992)). With high accuracy requirements, 
however, the CORDIC processor will have long delays, while the table-lookup 
method vnll certainly require a very large ROM. In this case, we propose to use 
the rectangular-to-polar converter which v^U be discussed in Section 7. This 
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rectangular-to-polar converter requires two small ROMs and it consolidates the 
operations into small array-multipliers, which can yield a smaller and faster circuit 
using well-known efficient multiplier implementation techniques. 

A synchronizer 6100 for implementing the synchronization scheme 
5 described above is illustrated in Figure 61. The synchronizer 6100 includes a 

correlator 6102, a Fourier Transform module 6104, and a rectangular-to-polar 
converter 6 1 06. The Fourier transform module 6 1 04 includes various delays and 
adders that are known to those skilled in the arts. The rectangular-to-polar 
converter is described further in Section 7. 
1 0 The synchronizer 6 1 00 receives data samples associated with sampling one 

or more received symbols and determines an offset 7rji/2, where \i represents a 
synchronization offset of the data samples relative to the incoming symbols. The 
operation of synchronizer 6100 is described in reference to the flowchart 6200, as 
follows. 

15 In step 6202, a set of complex data samples is received. 

In step 6204, the correlator 6 1 02 correlates the complex data samples with 
a complex conjugate of a preamble data set (a„*), resulting in correlated complex 
data samples. 

In step 6206, the Fourier transform module 6104 determines the Fourier 
20 transform of the correlated data samples signal, according to equations (6.10) - 

(6. 13) and related equations; 

In step 6208, the Fourier transform module 6104 evaluates the Fourier 

transform of the correlated data samples at 7r/2, generating a complex signal 

representing a complex number; 
25 In step 62 1 0, the rectangular-to-polar converter 6 1 06 determines an angle 

in a complex plane associated with the complex number of step 6210, where the 

angle represents synchronization between the data samples and the incoming 

symbols. 

In step 6212, the angle from step 6210 is scaled by 2/n to determine the 
30 synchronization offset. 
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6,4 Bias in Symbol Timing Estimation due to Truncating the 
Sequence 



By truncating the sequence r^(ri) before taking the Fourier transform, we 
have produced a very simple structure to compute \i. However, since Rj(e^) 
5 differs from we must determine how this difference would affect the 

estimated value. The truncated sequence rj(ri) is related to the original 
sequence r^(w) as 

^ri^) = rjn)win) (6.17) 
where w(n) is a rectangular function whose Fourier transform Wie''^ is a sine 
10 function. Thus, 

Rj.{e''')= R{e''')® w{e'''), (6.18) 
Taking the expectation of (6.18) we have 



(6.19) 



Obviously, the \i value obtained using Rj(el'^ in (6.16) would be different from 
1 5 that obtained usingi?(e''^. This will introduce a non-zero timing-jitter mean (bias) 

to the \i value obtained using /J^O instead of /?(e"^. But the phase difference 
of the expected values of /?7<e"^ and R(e^"^ can be computed for a given ^0 
The procedure is as follows: 

1. Given the pulse waveform g(/), obtain, for each value n, the 
20 samples g(w - ^i), w = -1, 2. 

2. Compute iJX^'^ using these samples^/? - x) according to (6. 15). 

3 . Find the value \i according to (6. 16). The difference between the 
desired value ii and the value jl computed using finite samples 
g(n - |i), w = -1, 2, is the bias. 

25 This bias is illustrated in Figure 63, where g{t) is a raised cosine Nyquist pulse 

with roUoff factor oc == 0. 1 . 
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From Figure 63, the bias is a function of and it can be precalculated and 
stored in a ROM in the receiver. Although the size of the ROM depends on the 
precision of ji, for typical precision requirements on |i the ROM can be quite 
small Let us illustrate this point using an example: If an 8-bit accuracy is desired 
for the bias, the bias value corresponding to the three most significant bits (MSBs) 
in li is indistinguishable from that corresponding to the full-precision \i value. 
Hence, we can use only the 3 MSBs in to determine the bias, thereby needing 
only 8 words in the bias lookup table. 

Thus, for each of our symbol timing detector output samples, we can 
obtain the corresponding bias value fi-om the ROM table, then subtract this bias 
fi-om the original timing estimate to obtain an unbiased estimate. 

We have thus far restricted our discussion to the timing recovery algorithm 
for 0 = 0. We now consider how this algorithm can be made to accommodate an 
arbitrary carrier phase 0. 

6. 5 Carrier-Independent Symbol Timing Recovery 

According to (6.6), with the T^=l normalization, the complex correlator 
output rj(n) is dependent on 0. Although the expectation of its magnitude 



does not depend on 0, it is non-trivial to compute the magnitude of rj(ji) from its 
real and imaginary components. Expressing r^(n) in terms of its real and 
imaginary components, according to (6.6), we have 




(6.20) 




(6.21) 



Thus, 




(6.22) 



(6.23) 
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20 



Since the carrier phase 0 does not depend on ^ we can treat it as a 
constant scale factor in Re[r^(ri)] md Im[r^(n)] when we are only concerned with 
extracting the timing information. 

Clearly, therefore, instead of using the magnitude of the complex 
value, we can use one of its real and imaginary parts, which are available at the 
output of the preamble correlator. 

We, of course, must decide which of Re[rj(ri)] and Im[r^(n)] to use. If 
the unknown phase 8 is such that cosG » 0 it is certainly desirable to use Im[r^(n)] 
instead of Re[r^(ny], and vise versa. But we don't know the 0 value thus far. 
How do we decide which one to use? 

From (6.22) and (6.23) we can see that the relative magnitudes of cos0 
and sin0 can be obtained from the real and imaginary components of rj(ri). For 
example, i{\Re[r^(ri)]\ > \Im[r^(n)]\ we certainly have that |cos0| > |sin0|, thus we 
should use the real part of the correlator output to find |i. Henceforth we denote 
the appropriate (real or imaginary) part of r^(n) by r^(rt). 

6. 6 Carrier Phase Computation 

Next, let us examine the problem of extracting the carrier phase. From 
(6.6) we can see that the phase of the complex number £'[r^(w)] does not depend 
on n. Moreover, the carrier phase can simply be obtained by extracting the phase 
of rjiri). In order to achieve the best precision, it is desirable to choose the r^{ri) 
value with the largest magnitude for carrier phase estimation. For example, if 
r^(«o) is the correlator output with largest squared-magnitude, we choose rj[n^ 
to compute 



One advantage of this approach is that the symbol timing and carrier phase 
estimations are independent of each other. They can thus be carried out in parallel 




(6.24) 
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As for symbol timing estimation in (6.16), the computation in (6.24) can 
be accomplished efficiently using the rectangular-to-polar converter to be 
discussed in Section 7. 

A synchronizer 6400 for determining timing and phase offsets is shown in 
FIG. 64. Similar to synchronizer 6100, the synchronizer 6400 receives data 
samples associated vsath sampling one or more received symbols and determines 
a timing oflFset 7rji/2, where fi represents a synchronization offset between the 
data samples and the incoming symbols. Additionally, the synchronizer 6400 
determines a carrier phase offset represented by 0. The synchronizer 6400 includes 
the correlator 6102, sample selectors 6404 and 6406, the Fourier transform 
module 6104, and two rectangular-to-polar converters 6106. The operation of 
synchronizer 6400 is described in reference to the flowchart 6500 in FIGs. 65 A-B, 
as follows. The order of the steps in flowchart 6500 is not limiting, as one or more 
steps can be performed simultaneously or in a different order, as will be 
understood by those skilled in the relevant arts. 

In step 6502, a set of complex data samples is received. 

In step 6504, the correlator 61 02correlates the complex data samples with 
a complex conjugate of a preamble data set (a^*), resulting in correlated complex 
data samples. Each correlated complex data sample includes a real sample and an 
imaginary sample. 

In step 6506, the sample set selector 6404 selects either the set of real 
correlated samples or the set of imaginary correlated samples. In embodiments, the 
set v^th the larger magnitude is selected. 

In step 6508, the Fourier transform module 6104 determines the Fourier 
transform of the selected real or imaginary data samples, according to equations 
(6. 10) -(6. 13) and related equations; 

In step 6510, the Fourier transform module 6104 evaluates the Fourier 
transform at nil, generating a complex signal representing a complex number; 

In step 6512, the rectangular-to-polar converter 6 1 06a determines an angle 
in a complex plane associated with the complex number of step 6510, where the 
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angle represents synchronization between the data samples and the incoming 
symbols. 

In step 6514, the angle from step 6512 is scaled by 2/n to determine the 
synchronization offset. 

5 In step 6516, the selector 6406 selects the largest correlator complex 

output. This selection can be based on an examination of one of the parts (real, 
imaginary) of the data sequence. 

In step 65 1 8, the rectangular-to-polar converter 6 1 06b determines an angle 
in a complex plane associated with complex output of step 65 16, where the angle 
10 represents the carrier phase offset 6. 

6, 7 Simulation Result 

We have used the above procedures to estimate the timing delay and the 
carrier phase of binary PAM symbols. The pulse shape was raised cosine with 
roUoff factor a = 0.4. The block size was L = 32 preamble symbols. To 
15 demonstrate its performance for signals with small excess bandwidth, we also 

tested this method with a = 0. 1 . For a carrier phase offset 0 = 45 °, we ran the 
simulation repeatedly using the synchronizer 6400, each time using a \i value 
randomly chosen between 0 and 1 . 

In addition to synchronizer 6400, we have also used the following two 
20 well-known methods to estimate the sampling mismatch: 

1 ) the DFT-based square-timing recovery (Oerder M. , and Meyr, H. , 
IEEE Tram. Comm, 36:605-612 (1988)), 

2) the method of (Gardner, S., Electron. Eng. 77:75-83 (Dec. 
1999)) that maps rJrtQ + l)/r^(wo)— the ratio of the two 

25 correlation values nearest the peak (see Figure 60) — to the 

sampling mismatch value |a. 
The variances of the timing jitter using these estimation methods for 
a = 0.4 and a = 0.1 are plotted in Figure 66 and Figure 67, respectively. The 
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corresponding Cramer-Rao bounds (CRB) — the theoretical lower bounds of 
estimation errors (Meyr, H., ei al.. Digital Communication Receivers: 
Synchronization, Channel Estimation and Signal Processing, Wiley, New York, 
NY (1998)) — are also shown. We can see that, in both cases, the timing- jitter 
5 variance using the proposed synchronizer is quite close to the theoretical bound. 

It clearly outperforms the other two methods, even for signals employing small 
excess bandwidth, as seen in Figure 67, 

The variance of the phase estimation error is depicted in Figure 68. It 
shows that, using the proposed method, the phase estimation error agrees quite 
10 well with the theoretical bound. 

6.8 Conclusion 

A synchronizer for initial symbol timing and carrier phase estimation using 
preambles has been presented. This synchronizer requires just two samples per 
symbol. Since the two estimations are independent of each other, they can be 

15 carried out simultaneously. These characteristics would ease the demand for 

computational speed for high-data-rate applications. Moreover, this synchronizer 
has demonstrated very good timing estimation performance even for signals with 
small excess bandwidth, which is essential for bandwidth eflBcient communications. 
The parameter estimations can be implemented very eflBciently using the 

20 synchronizer 6400. Due to its simplicity, this method is attractive for applications 

where low power and low complexity are desired, such as in a hand-held 
transceiver. 

7. A High-Speed Processor for Rectangular-to-Polar Conversion 

25 As discussed previously, the rapid acquisition characteristic of feedforward 

S5mibol synchronizers is essential to symbol synchronization for burst modems. 
Many feedforward structures require the evaluation of the phase of a complex 
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number. That is, an efBcient implementation of the phase extraction process is 
crucial. In order to handle a wide range of communications problems (Section 8), 
a general rectangular-to-polar conversion problem is considered. 

There are several well-known implementations for a rectangular to polar 
5 coordinate conversion, i.e. obtaining the magnitude and phase of a complex 

number. One method uses a ROM lookup table with both the real and imaginary 
components as input. This is practical only for low bit-accuracy requirements, as 
the ROM size grows exponentially with an increasing number of input bits. To 
reduce the ROM size, we can first divide the imaginary by the real component, 

1 0 then use the quotient to index the lookup table. But the hardware for a full-speed 

divider is very complicated and power consuming. An iterative divider 
implemented using shifting and subtraction requires less hardware, but it is usually 
quite slow. Recently, CORDIC has been applied in this coordinate conversion 
(Chen, A., and Yang, S., "Reduced complexity CORDIC demodulator 

1 5 implementation for D- AMPS and digital EF- sampled receiver," in Proc, Globecom 

J 998, vol 3 (1998), pp. 1491-1496). However, due to the sequential nature of 
CORDIC, it is diflBcult to pipeline, thus limiting the throughput rate. 

In burst-mode communication systems, rapid carrier and clock 
synchronization is crucial (Andronico, M., et aL, "A new algorithm for fast 

20 synchronization in a burst mode PSK demodulator," in Proc. 1995 IEEE Int. 

Conf, Comm,, vol, 3 (June 1995), pp. 1641-1646). Therefore, a fast rectangular- 
to-polar conversion is desired. In this section, we present an apparatus and method 
that implements the angle computation for rectangular-to-polar conversion with 
low latency and low hardware cost. This processor and the polar-to-rectangular 

25 processor presented in Section 5 (See rotator 3900 in FIG. 39), together, can 

perform the M-ary PSK modulation devised in (Critchlow, D., "The design and 
simulation of a modulatable direct digital synthesizer with non-iterative coordinate 
transformation and noise shaping filter," M.S. thesis, University of California, San 
Diego (1989)). 
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7. 1 Partitioning the Angle 

Figure 69 displays a point in the Cartesian X-Y plane having coordinates 
(Xq, ^0)3 wherein Xq and Yq represent the real and imaginary parts of an input 
complex signal. The angle 9 can be computed as 

5 <p=tan-'{YJX^). (7.1) 

In deriving the core of our algorithm, we assume the dividend and divisor satisfy 

X^>Y^> 0. (7.2) 
We will discuss how to extend the result to arbitrary values in Section 7.4. To 
achieve the highest precision for given hardware, the inputs and Yq should be 
10 scaled such that 

1 < < 2. (7.3) 
A straightforward method for fast implementation of (7. 1) can be devised 
as follows: 

1) Obtain the reciprocal of Xq from a lookup table. 
15 2) Compute Y^ x (IIX^ with a fast multiplier. 

3) Use this product to index an arctangent table for (p. 
However, the size of the two tables grows exponentially with increased precision 
requirements on (p, and rather large tables would be required to achieve accurate 
results. Therefore, for high-precision applications, such an implementation seems 
20 impractical. 

If we approximate MX^ by the reciprocal of the most significant bits 
(MSBs) of Xq, denoted by [Xq], then the required reciprocal table is much smaller. 
We can then multiply the table output by Yq to yield YJ\X^, which is an 
approximation of Y^l X^. This quotient can then be used to index an arctangent 
25 table. Similar to the reciprocal table, a much smaller arctangent table is needed 

if we use only the MSBs of VTO. denoted by [lV[-^o]]> to address the table, 
which returns <pi = tan**( [Kq / V^i^] ) Obviously, this result is just an 



1904.0140002 



-116- 



approximation to (p. We will subsequently refer to the computation of (pi as the 
coarse computation stage. 

Let be the diflFerence between cp and (p^ Using the trigonometric 

identity 

tan (p2 = tan ((p - (pj) = (tan (p - tan cpj) / (1 + tan cp x tan (pj (7.4) 



and the definitions tan (p == YJX^ and tan cpi = [IVLA^oU. we have 



7o/jro-[yo/[^o]] 


Yq~ Xq)(. 




u{YJX,y[Yj[x,Y 







(7.5) 



Using this relationship, (pj can be determined fi"om [lV[-^t^o]]> the coarse 
computation results. Therefore, the desired result (p can be obtained by adding the 
fine correction angle (p2 to the coarse approximation (pi . This procedure of finding 
(P2 will subsequently be referred to as the fine computation stage. 

By partitioning the computation of (7. 1) into two stages, the table size in 
the coarse stage can be reduced significantly at the expense of additional 
computations, which are handled by the fine stage. Let us now examine the 
complexity of the fine stage. To find (p2> we can first compute 



y, = y,-x,x[yj[x,\^ 

and then find (pj as 

(p^ = tan"^(7i/Xi). (7.7) 
The computation in (7.6) involves only adders and multipliers, while (7.7) 
requires lookup tables. Moreover, it seems we can't use the same coarse-stage 
tables because they have low resolution and thus can't satisfy the high precision 
requirements for the fine angle (pj. Now let us analyze cpj to see if there is any 
property that can help in this situation. 
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If (pi is a good approximation of cp, then cpj = (p - cpi is close to zero * In 
view of (7.7), YJX^^ should be very small too. This property helps us in two 
respects: 1) The difference between Y^IX^ and YJ[X^^ is much smaller than that 
between lAYj and \I[X{\. This suggests that if we use the same low resolution 
reciprocal table as in the coarse stage, the error contributed to the final result will 
be very small. We will demonstrate this in the next section. 2) If YJX^ is 
sufficiently small to satisfy 

|?;/A^J= |tan^2|< 2"^'" (7.8) 
where A'^ denotes the desired number of bits in (p, then 

(p^ = tan-'(y; I X,)^YJ X, (7.9) 
and we can compute (p2 without using an arctangent table. This is explained as 
follows: 

From the Taylor expansion of XdLn^{Y^IXy) near Y^IX^ = 0, we obtain 

tm~'{Y, /X,)=Y, '/ X, - {y, I X^ / 3 + o{^, I J^J'). (7.10) 
Since OifJxIX^^) is negligible in comparison to {JylX^IZ, it can be omitted. 
Therefore, if Y^IX^ is used to approximate tan'^(7i/Xi), an error 

A ^ = tan-^ (l^ / X, ) - 7, / X, = - ( / X, )' / 3 (7.11) 
will occur. However, according to (7.8), ts.^ is bounded by 

Atan|< 2-^/3 (7.12) 

which is very small. This indicates that the approximation (7.9) is quite accurate 
if (7.8) is satisfied. 

From the above analysis, no additional tables are needed for the fine stage 
if (Pi is sufficiently close to (p. On the other hand, the better that cpi approximates 
<p, the larger the tables required for its computation become. As mentioned 
previously, table size grows exponentially as the precision increases. A good 
trade-off is obtained when the result cpi of the coarse stage is just close enough to 
<p that (7.8) is satisfied, thereby eliminating the additional tables in the fine stage. 
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A detailed description of a rectangular-to-polar converter that implements the 
algorithm follows. 

FIG. 71 illustrates a rectangular-to-polar converter 7100 that implements 
the coarse and fine rotation described in section 7 herein, including equation (7. 1)- 
(7.53). The converter 7100 receives a complex input signal 7102 (that represents 
a complex number having Xq and Yo components) and determines the angle (p, 
which represents the position of the complex signal 7102 in the complex plane. In 
doing so, the converter 7100 determines a coarse angle computation that is 
represented by the angle <Pi, and performs a fine angle computation represented 
by the angle (pj. Once <Pi is determined, the input complex number 7102 is 
conceptually rotated back toward the X-axis to an intermediate complex signal 
71 15 as represented in FIG. 72, and (pj is determined fi-om intermediate complex 
signal 7115. The angles (p^ and (p2 are added together to determine cp. 

The converter 7100 includes: an input mux 7104, reciprocal ROM 7106, 
output demux 7108, an arctan ROM 71 10, a multiplier 71 12, a butterfly circuit 
71 14, a scaling shifter 71 16, a fine angle computation stage 7124, and an adder 
7126. The fine angle computation includes a multiplier set 7118, a one's 
complementer 7120, and a multiplier 7122. 

The ROM 7106 stores reciprocal values of [X^], wherein [Xq] is defined 
as the most significant bits (MSB) of Xq of the input signal 7102. The reciprocal 
of [Xq] is represented as Zq, for ease of reference. As will be shown, the ROM 
7 1 06 is re-used to determine the reciprocal of [XJ, where X^ is the real part of the 
intermediate complex number 7115 shown in FIG. 7 1 and FIG. 72. The reciprocal 
of [XJ is represented as Z^, for ease of reference. In embodiments, the ROM 7 106 
has 2^^'^^ storage spaces, where A^is the number of bits that represents Xq (and Yq) 
of the input signal 7102. 

The input mux 7104 chooses between [Xq] and [XJ as an input to the 
reciprocal ROM 7106, according to the control 7128. The output demux 7108 
couples an output of the ROM 7106 to Zq or Z^ according to the control 7128. 
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The control 7128 assures that Zq receives the stored reciprocal value for [Xq], 
and that Zj receives the stored reciprocal value for [X^]. 

The arctan ROM 7110 stores the coarse approximation angle (Pi based on 
a [YqZo] input. Therefore, a coarse stage can be described as including the ROM 
ROM 7110, the ROM 7106, and the multiplier 7112, as they are used in the 
coarse angle computation. 

The operation of the converter 7100 is described further with reference to 
the flowchart 7300, as follows. The order of the steps in the flowchart 7300 is not 
limiting as one or more of the steps can be performed simultaneously, or in 
different order. 

In step 7302, the input complex signal 7102 having a Xq component and 
a Yq component is received. In embodiments, the Xo and Yq components are N-bit 
binary numbers. 

In step 7304, the control 7128 causes Zq to be retrieved from the ROM 
7106, where Zq represents l/[Xo], and wherein [Xq] is the MSBs of X^, 

In step 7306, the multiplier 7112 multiplies Yq of the input complex 
number 7102 by Zq, resulting in a [ZoYq] component. The [ZqYo] component is an 
approximation of Yq/Xq 

In step 7308, the coarse angle <Pi is retrieved from the ROM 7110 based 
on [ZqYo], and is sent to the adder 7126. Note that the coarse stage can be 
described as including the ROM 7110, the ROM 7106, and the multiplier 7112, 
as they are used in the coarse angle computation. 

In step 73 10, the butterfly circuit 71 14 multiplies the input complex signal 
7102 by [ZqYo]. This causes the input complex signal 7102 to be rotated in the 
complex plane toward the real axis to produce the intermediate complex signal 
7115 (representing an intermediate complex number), having a real X^ component 
and an imaginary component. 

In step 73 12, the scaler 7116 scales the X^ component of the intermediate 
complex signal so that it is compatible with the reciprocal values stored in the 
ROM 7106. The scaler also scales the Y^ component by the same amount. 
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In step 7314, the control 7128 causes to be retrieved from the ROM 
7106 based on [XJ, where Zi represents l/[Xi], and wherein [XJ is the MSBs of 
Xp Note, that the ROM 7106 is efficiently used twice to calculate two different 
reciprocals Zq and Zj, thereby reducing overall memory size. 

In step 7316, the fine angle computation stage 7124 determines the fine 
angle (p2 based on Z^ and the scaled intermediate complex number 71 1 5. In doing 
so the Newton- Raphson method is emulated in hardware to estimate cp2, which 
is the arctan of ZiYi. More specifically, multiplier set 71 1 8 multiples X^,Y^ by Z^. 
The ones' (approximating two's) complement 7120 is then determined for X^Z^. 
After which, the multiplier 7127 multiplies (2-XiZi) by YiZj, to determine tan 
cp2. Since cpj is a small angle, the value tan is used as an approximation of (pj 

In step 7316, the cpi and (pj are added together to get (p. 

A more detailed description of the algorithm follows. 

7. 2 The Two-Stage Algorithm 

In this section we first analyze how the coarse approximation error (pj = 
q) - (Pi depends upon the precision of the tables 7106 and 7110, in order to 
determine the amount of hardware that must be allocated to the coarse stage. 
Next we explore ways to simplify the computations in the fine stage. 

7,2.1 Simplification in the Coarse Computation Stage 

The main concern in the coarse stage design is how the lookup table values 
are generated to produce as precise results as possible for a given table size. As 
mentioned previously, there are two lookup tables: 
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7.2. LI The Reciprocal Table 7106 
The input to this table, 1 <. Xq<2^ can be expressed as 

where only bits through x„ are used to index the table. To generate the table 
5 value, if we merely truncate as 

{X^\ = \,X^X^...X^ (7.14) 
then the quantization error A = - \Xq ] is bounded by 

0< A;^^ < 2"'". (7.15) 
Thus, the difference between the table value and I/Xq, 

10 \/X,-\l[x,]= ([Xo]- Xo)/([Xo]Xo)« -t^xJXl (7.16) 
is bounded by 

-2"'" < l/[^o]^ 0- (7 17) 
But if we generate the table value corresponding to 

[X^]= \,X,X^...X^\ (7.18) 

15 with a bit "1" appended as the LSB, then the quantization error in (7.15) is 
centered around zero: 



-2"'""' < A;^^ < 2"^""' (7.19) 
hence, the error in the reciprocal is also centered around zero: 



_2— 1 <\l X^- l/[Xo]< 2-"-\ (7.20) 
20 Comparing (7.20) to (7. 17), the maximum absolute error is reduced. This is the 

technique introduced in (Fowler, D.L., and Smith, J.E., "An accurate high speed 
implementation of division by reciprocal approximation," in Proc. 9th Symp. on 
Computer Arithmetic (19S9X pp. 60-67). 
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Since the output of the table will be muhiplied by Yq, the fewer the bits in 
the table value, the smaller the required multiplier hardware. Let the table value 
Zq be generated by rounding l/[Xo] to n bits: 

Zo = 0.1Z2Z3...Z„. (7.21) 
The quantization error A = 1 / [x^ ] - is then bounded by 

A < 2—\ (7.22) 
Once we have obtained Zq from the reciprocal table, we can get an 
approximation to the quotient Yq/Xq by computing Yq^q, This result is then used 
to address the arctangent table for tpj. 



7.ZL2 The Arctangent Table 7110 

In order to use a very small table, Y^ Zq is rounded to k bits to the right of 
the radix point to become [Y^ Z^], with the rounding error bounded by 

-2"*-^ < ^r,z, = Y,Z, - [y,Z,] < 2-'-\ (7.23) 
Then, [Y^ Zq] is used to index the arctangent table, which returns the coarse angle 
(Pi = tan'^([7oZo]). 

Now we must determine the minimum /w, n and k values such that (7.8) is 
satisfied. First, let us examine and Y^ which are computed using [Yq Zq] as 

= ^0 ^ YoIYoZo] (7.24) 

yi = ' ^o[yoZo] (7-25) 
Dividing (7.25) by (7,24), and then dividing both the numerator and 
denominator by Xq, we have 
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(7.27) 



\y,/X,\= (Fo / X, - [y,Z,])/(i + (Fo / FoZo ) 

I r 1 

<|7o/;^o-[>^o^o]- 

The inequality is true because A'q ^ Yq^O and [ro^o] ^ 0 Taking into account all 

the quantization errors in (7.20), (7.22) and (7.23), we can express Y^/Xq in terms 

of [Ko^o] as 

= Yo{{z, + A,yA^JX^) 
= Y,Z, + Y,A,^-Y,{A^JXi) 

= [YoZo] + A r,z, + Yo^ - Yoi^ I )■ 
Substituting this result into (7.26), we have 

k / ^il ^ h ^0^0 + ^'oA z„ - J^o(a I (7.28) 
Since 7o(a / X^') = {y^ I X„)(a / x)j , from (7.2) and (7.19), 

-2-"'-'<7„(a^^/^o^)<2-'"-'. (7.29) 
Also, according to (7.2) and (7.22), we have 

-2-"<7oA^^ <2-". (7.30) 
Applying (7.23), (7.29) and (7.30) to (7.28), we obtain \Y,/X,\ < T"-' + T" + T^^. 
If we choose /M^iV^/3 + l.«^A^/3+2andA:^7//3 + l, then 

l^/Jr,|<0.75x2-^'^ (7.31) 
Therefore, since the inputs A', and to the fine stage satisfy (7.8), no additional 
tables are needed for the fine stage. Henceforth we choose m = N I 3 + 1, 
« = iV/ 3 + 2 and yt = iV^/ 3 + 1. 
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7. 2, 2 Hardware Reduction in the Fine Computation Stage 
7124 



Since (7.8) is satisfied, we can obtain the fine angle ^2 by computing the 
quotient Y^IXy. From (7.24), we haveXp <. X^<.Xq + Y^, hence \ <. Xy<A. In 
5 order to use the same reciprocal table as in the coarse stage, X^ should be scaled 

such that 

\^X^<2, (732) 
This can be satisfied by shifting ^1 to the right ifXi ^ 2. Of course should also 

be shifted accordingly so that Y^/X^ remains unchanged. 

10 As in the coarse stage, the reciprocal table accepts A^/ 3 + 1 MSBs ofX^ 

and returns Z^. We define the reciprocal error 6^ = l/Xy Z^. Since the same 

reciprocal table is used as in the coarse stage, 6^ and 60 must have the same 

bound. Since 

S, = S, = \/X,-Z, = l/X,-l/[X,-\+^,^ (7.33) 
15 we can use (7.20) and (7.22) to obtain 

-0.75 X 2-"^'^-' <S,< 0.75 x 2-"^'^'. (7.34) 
The bound on Y^ can be found using (7.31) and (7.32): 

|7J< 0.75 X 2"'^''^^ (7.35) 
Now we can obtain the final error bound in approximating X^/Y^ by Y^Z^^ 
20 according to (7.34) and (7.35), as 

\Y,/X, - Y,X,\=\Y,S,\< (0.75)' x 2"'^''. (7.36) 
Clearly, this approximation error is too large. To reduce the maximum 
error below 2'^, the bound on IdJ should be approximately 2"^^, which would 
require the reciprocal table to accept 2M3 bits as input. That is, the table needed 
25 for such a high-resolution input would be significantly larger than the one already 

employed by the coarse stage. 

To overcome this difficulty, we can apply the Newton-Raphson iteration 
method (Koren, L, Computer Arithmetic Algorithms, Prentice Hall, Englewood 
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ClifFs, NJ (1993)) to reduce the initial approximation error 61. We now briefly 
explain how this method works. First, let us define the following function: 

/(Z)=l/Z-Jr,. (7.37) 
Obviously, we can obtain Z=yx^ by solvingy(Z) = 0. In other words, we can find 
\IXi by searching for the Z value such thaty(Z) intersects the Z-axis, as shown in 
Figure 70A. ^ 

Shifting the Z-axis down by X^, we obtain a new fiinction /i(Z) = — , 
shown in Figure 70B. At = \I[X{\ , being the initial guess, the slope of 
A(Z)-\IZ is 

The tangent, shown as the dashed line 7102, intersects the/,(Z) = A', line at a new 
point Z2. From Figure 70B, Z2 is much closer to the desired value IIX^ than the 
initial guess Z,. Let us now find Z2. According to Figure 7-3, we must have 

^■-'^^-.--L. (7 39) 

Expressing in terms of and we have 

Z2 = Zi(2-XiZi). (7.40) 
Thus, we can obtain Zj, a more accurate approximation of l/Xj, from Z^. One may 
wonder how accurate Zj is in approximating \IX^, Let us now examine the 
approximation error bound. 

Substituting Z^ = \IX^-b^ into (7.40), we have 

Z, = (1 / X, - S,){2 - X,{\IX, - S,)) 
= 1/X,-X,S^, 

According to (7.32), (7.34) and (7.35), after one Newton-Raphson iteration, the 
error in K1Z2 is reduced to 

\Y,/X, - Y,Z^\ = \Y,X,S^\< (0.75)' x 2"^. (7,42) 
Thus, a rather accurate result is obtained with just one iteration. 
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Finally, the fine angle can be computed by multiplying by Y^ : 

q>^ « Y,Z^ = Y,Z,(2- X,Z,\ (7.43) 
Although there are three multipliers involved in (7.43), the size of these 
multipliers can be reduced with just a slight accuracy loss by truncating the data 
before feeding them to the multipliers. The computational procedure of (7.43) is 
as follows: 

1 ) The inputs to the fine stage, and 7^, are truncated to 2A73+2 and 
N+3 bits to the right of their radix points, respectively. Since the M3-1 MSBs in 
7i are just sign bits, as indicated by (7.35), they do not influence the complexity 
of the muhiplier that produces Y^Z^, The corresponding quantization errors are 
bounded by 

0< A;,^<2-^^^^-^ (7.44) 
0<A^^<2-^-^ (7.45) 

2) Both quantized and Y^ are multiplied by Zj. 

3) To form 2- X^Z^, instead of generating the two's complement of 
XjZi, we can use the ones complement v^th only an insignificant error. Since this 
error is much smaller, in comparison to the truncation error in the next step, we 
can neglect it. 

4) The product Y^Z^ is truncated to N+3 bits. We would also truncate 
the ones complement of X^Z^, But since the inverted LSBs of J^^Z^ will be 
discarded, we can truncate X^Z^ to 2M3+2 bits and then take its ones 
vcomplement. The corresponding quantization errors, as discussed above, are: 

0<^^^,^<2-'^"-' (7.46) 

0< A^,^, <2-^-^ (7.47) 
After including all the error sources due to simplification, we now analyze 
the effects of these errors on the final result cpj. Taking the errors into account, we 
can rewrite (7.43) as: 
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^((I'l-Aj, )(l/jri -'5i) + Aj.Zj)(2-(>'^i-A;^.^)(l/jrj-<yj) + A;j.^^^) (7.48) 

Expanding this product and neglecting terms whose magnitudes are insignificant, 
we have 

-{\IX, )A + (y; / X, )(A ) + A ,,,, . ^^^^^ 
As mentioned in Section 7. 1, FjAYi is an approximation of tan'^CFj/A^j). Its 
approximation error, defined in (7.8), is bounded by 

|A..„| = \(Y^/X{f/3\ < (0.75)3 X 2-^/3. (7.50) 
Replacing Y^/Xi by tan"' (K,/^,) + (y,/Jr,)'/3 in (7.49), we have 



« tan-'(7, / J,) + (7, /X,f/3- + (7, / jf )A;,, 
-(l/X,)AK, + (i;/^,)A^.^. + A,,,. 
10 The total error, e = cpj - tan * (FiAY,), is 



(7.51) 



^=(Y,/ xM(Y, /X,f/3- {XAf + A X, + A } 

V i I 1 ^ (7.52) 



\2 



All terms in the subtotal (7^ / X^f 1 3 - (^i^i)^ ^ I ^\ + A 2 ^® non- 
negative. Thus, the lower bound of this subtotal is the minimum value of 
-iX,S,y, which is -0.75^ X 2-^'' = -0.56 x 2-^'^, according to (7.34). 
1 5 Correspondingly, its upper bound is the sum of the maximum values of the other 

three terms: (0.75^/3 + X'^ + 2"^) x 2'^'^= 0.68 x 2'^'^. 
Finally, we can obtain the total error bound as: 



|£|< 0.75 X 2"^'^ X 0.68 X 2-"^'" + 2"^-' + 2"^-=' 
= 0.76 X 2 



N (7-53) 
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7.3 Magnitude Calculation 

Once the angle of the vector {X^^Y^ is known, its magnitude can be 
obtained by multiplying Xq by 1/cos <p, whose values can be pre-calculated and 
stored in a ROM, thereby requiring only a single multiplication. However, if we 
use all the available bits to index the ROM table, it is likely that a very large ROM 
will be needed. 

As we know from the preceding discussion, the coarse angle <Pi is an 
approximation of cp. Similarly 1 /cos g)^ approximates 1/cos (p. Therefore, we can 
expand the coarse-stage ROM 71 10 to include also the 1/cos (p^ values. That is, 
for each input [Y^^], the coarse-stage ROM would output both q)y^ = tan"^ ([^^o^o]) 
and 1/cos q)^. Since and Yq satisfy (7.2) and (7.3), the 1/cos Rvalue is within 
the interval [1,V2]. 

For many applications, the magnitude value is used only to adjust the 
scaling of some signal level, and high precision is not necessary. For applications 
where a higher precision is desired, we propose the following approach: 

First, instead of using the above-mentioned table of 1/cos q)^ values, 
we pre-calculate and store in ROM the 1/cos <Pi^ values, where contains only 
the m MSBs of <p. Obviously a small table, one of comparable size to the l/coscp^ 
table, is needed. Then, we can look up the table entries for the two nearest values 
to q), namely and <Pm^ <Pm^ 2"'" . Then a better approximation of 1/cos <p 
can be obtained by interpolating between the table values 1/cos and 1 / cos^j^ 
as 



Let <Pj^^ (p- q)^. Obviously, simply contains the LSBs of q). We can now 
rewrite (7.54) as 



1/ cos^» Icos^^ + 



9m- (Pm 



^{(P^<Pm)' (7.54) 



1904.0140002 



-129- 



l/cos^« l/cos^2?^+(l/cos^;^-l/cos^2>^)x X 2^" (7.55) 

which involves only a multiplication and a shift operation, in addition to two 
adders. 

7. 4 Converting Arbitrary Inputs 

In previous sections we have restricted the input values to lie within the 
bounds of (7.2) and (7.3). However, if the coordinates of (X^J^ do not satisfy 
that condition, we must map the given point to one whose coordinates do. Of 
course, the resulting angle must be modified accordingly. To do that, we replace 
JSTo and7o by their absolute values. This maps {X^J^ into the first quadrant. Next, 
the larger of lA'ol and \Yq\ is used as the denominator in (7. 1) and the other as the 
numerator. This places the corresponding angle in the interval [0, 7i/4]. We can 
now use the procedure discussed previously to obtain q>. Once we get ^, we can 
find the angle (|) that corresponds to the original coordinates fi-om q>. First, if 
originally [A'ol l^ol we should map <p to [ttM, 7r/2] using <f>' = n 12 - (p . 
Otherwise (f>' - cp . We then map this result to the original quadrant according to 
Table 7.1. 

Next, let us examine the affect of the above-mentioned mapping on 
the magnitude calculation. Since the negation and exchange of the original 
and To values do not change the magnitude, whose value is (^Xl + Y'^ , the 
result obtained using the X^ and values after the mapping needs no correction. 
However, if the input values were scaled to satisfy (7.3), we then need to scale the 
computed magnitude to the original scale of J^o and Y^. 



Table 7.1 Converting Arbitrary Inputs 


Original coordinates 


Modification to the angle 


Xo<o,yo>o 


(j>- TV - (f>^ 
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Xo>o,yo<o 


^ = 2;r - ^' 



7.5 Test Result 

We have verified our error bound estimation by a bit-level simulation of 
the rectangular-to-polar converter 7100. To test the core algorithm described in 
Section 7.2, we generated the pair of inputs Xq and randomly within the range 
described by (7.2) and (7.3). This test was run repeatedly over many different 
values of Xq and Yq, and the maximum error value was recorded. Choosing N = 
9 for this simulation, the error bound estimate according to (7.53) is 0.0015. Our 
test results yielded the error bounds [-0.00014, 0.00051], well within the 
calculated bound. 

7. 6 Conclusion 

An efficient rectangular-to-polar converter is described. The angle 
computation of a complex number is partitioned into coarse and fine 
computational stages. Very small arctangent and reciprocal tables are used to 
obtain a coarse angle. These tables should provide just enough precision such that 
the remaining fine angle is small enough to approximately equal its tangent value. 
Therefore the fine angle can be obtained without a look-up table. The 
computations are consolidated into a few small multipliers, given a precision 
requirement. While a low-precision magnitude can be obtained quite simply, a 
high-precision result can be achieved by combining the angle computation with the 
angle rotation processor 3900 of Section 5. 

The applications of this converter include the implementation of the 
converter 6106 in the symbol synchronizer 6100 and the synchronizer 6400. 
However, the converter is not limited to sjonbol synchronization. It also provides 
eflScient implementation of computational tasks for many communication systems. 
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such as constant-amplitude FSK and PSK modems (Chen, A., and Yang, S., 
"Reduced complexity CORDIC demodulator implementation for D-AMPS and 
digital IF-sampled receiver," in Proc. Globecom 1998, vol. 3 (1998), pp. 1491- 
1496; Boutin, N., IEEE Trans. Consumer Electron. 38:5-9 (1992)), DMT 
modems (Arivoli, T., etal, "A single chip DMT modem for high-speed WLANs," 
in Proc. 1998 Custom Integrated Circuits Conf. (May 1998), pp. 9-11), as well 
as carrier synchronization (Andronico, M., et al, "A new algorithm for fast 
synchronization in a burst mode PSK demodulator," in Proc. 1995 IEEE Int. 
Conf. Comm.. vol. 3 (June 1995), pp. 1641-1646; Fitz, M P., andLindsey, W.C., 
IEEE Trans. Comm. ¥0:1644-1653 (1992) where the computation of phase and 
magnitude from the rectangular coordinates is essential. 

8. Exemplary Conq)uter System 

Embodiments of invention may be implemented using hardware, 
software or a combination thereof and may be implemented in a computer 
system or other processing system. In fact, in one embodiment, the invention 
is directed toward a software and/or hardware embodiment in a computer 
system. An example computer system 7702 is shown in FIG. 77. The 
computer system 7702 includes one or more processors, such as processor 
7704. The processor 7704 is connected to a communication bus 7706. The 
invention can be implemented in various software embodiments that can operate 
in this example computer system. After reading this description, it will become 
apparent to a person skilled in the relevant art how to implement the invention 
using other computer systems and/or computer architectures. 

Computer system 7702 also includes a main memory 7708, preferably 
a random access memory (RAM), and can also include a secondary memory or 
secondary storage 7710. The secondary memory 7710 can include, for 
example, a hard disk drive 7712 and a removable storage drive 7714, 
representing a floppy disk drive, a magnetic tape drive, an optical disk drive, 
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etc. The removable storage drive 7714 reads from and/or writes to a 
removable storage unit 7716 in a well known manner. Removable storage unit 
7716, represents a floppy disk, magnetic tape, optical disk, etc. which is read 
by and written to by removable storage drive 7714. As will be appreciated, 
5 the removable storage unit 7716 includes a computer usable storage medium 

having stored therein computer software and/or data. 

In alternative embodiments, secondary memory 7710 may include other 
similar means for allowing computer software and data to be loaded into 
computer system 7702. Such means can include, for example, a removable 

10 storage unit 7720 and an storage interface 7718. Examples of such can include 

a program cartridge and cartridge interface (such as that found in video game 
devices), a removable memory chip (such as an EPROM, or PROM) and 
associated socket, and other removable storage units 7720 and interfaces 7718 
which allow software and data to be transferred from the removable storage unit 

15 7720 to the computer system 7702. 

Computer system 7702 can also include a communications interface 
7722. Conmiunications interface 7722 allows software and data to be 
transferred between computer system 7702 and external devices 7726. 
Examples of communications interface 7722 can include a modem, a network 

20 interface (such as an Ethernet card), a communications port, a PCMCIA slot 

and card, etc. Software and data transferred via communications interface 7722 
are in the form of signals, which can be electronic, electromagnetic, optical or 
other signals capable of being received by the communications interface 7722. 
These signals are provided to the communications interface 7722 via a channel 

25 7724. This channel 7724 can be implemented using wire or cable, fiber optics, 

a phone line, a cellular phone link, an RF link and other communications 
channels. 

Computer system 7702 may also include well known peripherals 7703 
including a display monitor, a keyboard, a printers and facsimile, and a 
30 pointing device such a computer mouse, track ball, etc. 
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In this document, the terms "computer program medium" and "computer 
usable medium" are used to generally refer to media such as the removable 
storage devices 7716 and 7718, a hard disk installed in hard disk drive 7712, 
semiconductor memory devices including RAM and ROM, and associated 
5 signals. These computer program products are means for providing software 

(including computer programs that embody the invention) and/or data to 
computer system 7702. 

Computer programs (also called computer control logic or computer 
program logic) are generally stored in main memory 7708 and/or secondary 
5 10 memory 7710 and executed therefrom. Computer programs can also be 

J received via communications interface 7722. Such computer programs, when 

^ executed, enable the computer system 7702 to perform the features of the 

=P present invention as discussed herein. In particular, the computer programs, 

a when executed, enable the processor 7704 to perform the features of the present 

15 invention. Accordingly, such computer programs represent controllers of the 

W computer system 7702. 



^_ 'o.^ Aii-itt*,6mbodiment where the invention is implement using software, the 
^^software may hectored in a computer program product and loaded into 



computer system 7702 usingl-emovable storage drive 7714, hard drive 7712 or 
20 communications interface 7722. In&^cmitrol logic (software), when executed 

by the processor 7704, causes the processoK^04 to perform the ftinctions of 

the invention as described herein. 

In another embodiment, the invention is implemented primarily in 

hardware using, for example, hardware components such as application specific 
25 integrated circuits (ASICs), stand alone processors, and/or digital signal 

processors (DSPs). Implementation of the hardware state machine so as to 

perform the functions described herein will be apparent to persons skilled in the 

relevant art(s). In embodiments, the invention can exist as software operating 

on these hardware platforms. 
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In yet another embodiment, the invention is implemented using a 
combination of both hardware and software. 

9. Appendices 

The following Appendices are included. 

9. 1 Appendix A: Proof of the Zero ISI Condition. 

In Section 3.4 we have examined the condition for a real-valued function 
f(t) to have zero crossings at integer multiples of the sampling period, i.e., for 
satisfying (3 .16). We have stated that / (0 satisfies (3 .16) if and only if F{k) , 
defined as samples of F (Q) in (3.8) (F is the fi-equency response off), satisfy 
(3.18). Here, < A: < oo is an integer. We now provide the proof. 

Proof. First let us define a periodic, of period N, extension ofXO as 

/c(0= Y.f(t-Nny (A.1) 

n=-oo 

Its Fourier transform consists of a sequence of impulses 

^^(^)= z^(^K"-^4 ^^^^ 

^ = — CO 

Next, consider an impulse chain. 

c(t)= f^S{t-rj). (A.3) 

n=-oo 

whose Fourier transform is 

CO 

C(n) = ^ ^(Q - iTTk) . (A.4) 
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The convolution F^^C can be expressed as 

F,(Q)0C(n) 

QO 

m=-oo 

Substituting (A.2) into (A.5) yields: 
F,(fi)®C(n) 

Therefore, we have the following relationships 

(3 16) o fJt) = I ^ t = Nm,m an integer 
[O all other integers. 

m=—co 

00 

<» J] F{k-Nm)=\ 

m=-oo 

(3.18) 



(A.5) 



(A.6) 



This concludes the proof 
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9- 2 Appendix B: Impulse Response of the Simplified 
Interpolators 



In Section 2, after introducing a preliminary interpolation method, we have 
shown that we can trade one angle rotator for a multiplier by conceptually 
modifying the input samples, then by "correcting" the interpolated value obtained 
from the "modified" samples. A simpler implementation structure as well as better 
performance in interpolating most practical signals have been facilitated. We now 
derive the impulse response of this interpolation filter. As discussed in Section 2, 
the interpolated sample is computed as 



^ N(2-\ 

y{^>-h (B.i) 

where K is defined in (2.30), and 

N/2 N/2 
m^-N/2+l OT=-7//2+l 

where ^-0, ..,,N/2-l. 

Substituting (2.30) into (B.2), we have 
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Replacing K and c^. in (B. 1) by (2.30) and (B.3), respectively, we have 

1 Nfl-\ Nn f N/2 2 ^ 

2 



///2 / 1 N/2-1 



^ N/2 „ A//2-1 ^^ 



2 

- — c- ly 



Z ^ Z 

Vn=-A^/2+l k=-N/2^l > 



N/2 

I 

m=-A^/2+l 



where /(l^^, the impulse response of the simplified interpolation filter discussed in 
Section 2.6. 1, is now defined as 



1 N/2-\ 



^ k=- N/2+1 



9 f N/2-l N/2-\ \ C^^) 

-4(-in Z ^ Z 



for -/w ^ / ^ 1 - m, -M2+1,...,M2. Otherwise,y(0 = 0- 

The fi-equency response, of course, can be obtained by taking the Fourier 
transform of f[t). To modify the fi-equency response of f[tX we can multiply the 

coefficient in (B.3) by a value denoted by F/i(k) . In designing an optimal 

10 interpolation filter as discussed in Section 4, we search for the F/i(k) value that 

minimizes (4.4), i.e., FJu>) most accurately approximates the desired 
response (4.3). 
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9.3 Appendix C: Fourier Transform of g(nT^-fi) 

Since, in Section 6, the g(riT^ are samples of the continuous-time pulse 
g(0, assuming, without loss of generality, that = 1, it is well-known (Freeman, 
H., Discrete-Time Systems, Wiley, New York, NY (1965)) that the Fourier 
transforms are related as 

G(e''')^ X GicD-^lTik) (C.l) 
where G(e^'") and G(o}) are the Fourier transforms of g(«) and g(t\ 
respectively. Since ^0 is bandlimited, i.e., |G(ey)|= 0 for \co\> tt, we have 
G(e^" ) = G(o}X -7r<co<7r, (C.2) 

Using the Fourier transform' s time-shifting property, the Fourier transform of g(t- 

e'^'^Gio)), (C.3) 

Since the gin-ix) are samples of git'\i) , for the same reason as the above, their 
Fourier transforms are the same in the interval -7t^a)<7i, as in (C.2). Thus, 
according to (C.2) and (C.3) we have 

FT(gin - ju)) = e^"^G(e^" ), -7r<co<7r, (C.4) 

9. 4 Appendix D: Interpolation on Non-Center Intervals 

When we first discussed the interpolation problem in Section 2, we 
focused on interpolating between the two samples in the middle of the set of 
samples used to generate a synchronized sample. What is the impact on 
interpolation performance when we interpolate in an interval not at the center of 
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the samples being used? Figure 74shows such an example for N= 4, where the 
interpolation is performed between XO) and>'(l) using y(-2Xy('lXy(0) andXO 
(as opposed to using yiO)y andX2), as seen in Figure 2-1). 

Using the procedure described in Section 2, given samples y(n), 
« = -A^+ 2,..., 1, we first compute the Fourier coefficients as 

Comparing (D. 1) to (2.9), their only difference is the range of the summation. As 
in Section 2, for a given offset 0 ^ ji < 1, the synchronized sample y{\i) can be 
computed as: 

f N/2-i ^ 

(D.2) 



We can express ^^i) in terms of yiri) by substituting (D.2) into (D. 1), as 
1 Ink 



y^l^^^lTj Z X«) 1+2^ cos— r (//-«)+ cos ^(//-w) 
= 2 y(f^)fif^-n) 



where 



f(0 = 



1+ 2 > cos— — / + COSTT/- 1< t < ~ N 2 

T^i N (D.4) 

0 otherwise 



is the impulse response of the corresponding interpolation filter. For 4,/{t) is 
plotted in Figure 75A. Taking the Fourier transform of J^t), we obtain the 
corresponding fi*equency response, which is shovm in Figure 76A. 

Comparing the Figure 76 A fi-equency response to Figure 7A, both for 
we can see that the interpolation performance degraded significantly, as 
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shown by the ripples in the passband and large sidelobes in the stopband in Figure 
76A, when the interpolation is not done in the center interval. 

However, using the optimization method discussed in Section 4, we can 
"reshape" the impulse response M such that the corresponding frequency 
response is a better approximation to the ideal interpolation frequency response. 

The impulse response of an optimized interpolation filter for a non-center 
interval is illustrated in Figure 75B. The corresponding frequency response is 
shown in Figure 76B, which is clearly better than Figure 76A, since it has less 
ripple in the passband and more attenuation in the stopband. 

Using samples X-^ + 2), ...,y{-\), and^O), to interpolate obviously 
reduces the latency in generating synchronized samples, as compared to using 
samplesX-A^/2+ 1), . ,y(0), ...,y(N/2), since the interpolator does not have to 
wait until samples XI), ■■.,y(N/2), become available. In applications where low 
latency takes a higher priority than interpolation accuracy, this approach will be 
useful. 

9. 5 Appendix E 

The following documents are incorporated by reference in their entireties: 

1. Buchanan, K., etal.. lEEEPers. Comm. ¥:8-13 (1997); 

2. Reimers, U., IEEE Comm. Magazine i<J: 104-1 10 (1998); 

3. Cho, K., "A frequency-agile single-chip QAM modulator with 
beamforming diversity," Ph.D. dissertation. University of California, Los 
Angeles (1999); 

4. Oerder M., and Meyr, H., IEEE Trans. Comm. 3<J:605-612 (1988); 

5. PoUet T., and Peeters, M., IEEE Comm. Magazine 57:80-86 (1999); 

6. Gardner, F.M., IEEE Trans. Comm. ¥7:502-508 (1993); 

7. Erup, L., et al, IEEE Trans. Comm. ¥7:998-1008 (1993); 



1904.0140002 



-141- 

Farrow, C, "Acontinuously variable digital delay element," inProc. IEEE 

Int. Symp. Circuits Syst. (June 1988), pp. 2641-2645; 

Proakis, J.G., Digital Communications, McGraw-Hill, New York, NY 

(1993); 

Proakis, J.G., Digital Signal Processing, Macmillan, New York, NY 
(1992); 

Vesma, J., and Saramaki, T., "Interpolation filters with arbitrary 
fi-equency response for all-digital receivers," in Proc. 1996 IEEE Int. 
Symp. Circuits Syst. (May 1996), pp. 568-571; 

Ahn, Y., et al., "VLSI design of a CORDIC-based derotator," in Proc. 
1998 IEEE Int. Symp. Circuits Syst., Vol. //(May 1998), pp. 449-452; 
Wang, S., etal, IEEE Trans. Comp. 4&.Vl<!n.-\2'!yi (1997); 
Madisetti, A., etal, IEEE J. Solid-State Circuits 34:1034-1043 (1999); 
Tan, L. and Samueli, H., IEEE J. Solid-State Circuits 50:193-200 (1995; 
Madisetti, A., "VLSI architectures and IC implementations for bandwidth 
efficient communications," Ph.D. dissertation. University of California, 
Los Angeles (1996); 

VankksL, J., IEEE Trans. Ultrasonics, Ferroelectronics and Freq. Control 
'^^:526-534 (1997); 

Daneshrad, B., "System design of a 1 .6 Mbps all-digital QAM transceiver 
for digital subscriber line applications," Ph.D. dissertation. University of 
California, Los Angeles (1993); 

Tan, L., "High performance architectures and circuits for QAM 
transceivers," Ph.D. dissertation. University of California, Los Angeles 
(1995); 

Gw[<iner,S., Electron. Eng. 7/: 85-92 (1999); 
Gardner, S., Electron. Eng. 7/:75-83 (1999); 
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Sabel, L., and Cowley, W., "A recursive algorithm for the estimation of 
symbol timing in PSK burst modems," in Proc, Globecom 1992, vol 1 
(1992), pp. 360-364; 

Chen, A., et al, "Modified CORDIC demodulator implementation for 
digital IF-sampled receiver," in Proc. Globecom 1995, vol 2 (1995), pp. 
1450-1454 ; 

Boutin, N., IEEE Trans. Consumer Electron, 38:5-9 (1992); 

Meyr, H., et al. Digital Communication Receivers: Synchronization, 

Channel Estimation and Signal Processings Wiley, New York, NY 

(1998); 

Chen, A., and Yang, S,, "Reduced complexity CORDIC demodulator 
implementation for D-AMPS and digital IF-sampled receiver," in Proc. 
Globecom 1998, vol. 3 (1998), pp. 1491-1496; 

Andronico, M., et al. , "A new algorithm for fast synchronization in a burst 
mode PSK demodulator," in Proc. 1995 IEEE Int. Conf. Comm., vol. 3 
(June 1995), pp. 1641-1646; 

Critchlow, D., "The design and simulation of a modulatable direct digital 
synthesizer with non-iterative coordinate transformation and noise shaping 
filter," M.S. thesis. University of California, San Diego (1989); 

Arivoli, T., etal, "A single chip DMT modem for high-speed WLANs," 
in Proc. 1998 Custom Integrated Circuits Conf. (May 1998), pp. 9-11; 
Fitz, M P,, and Lindsey, W.C., IEEE Trans. Comm. ^0:1644-1653 
(1992); 

Fowler, D.L., and Smith, J.E., "An accurate high speed implementation of 
division by reciprocal approximation," in Proc. 9th Symp. on Computer 
Arithmetic (1989), pp. 60-67; 

Koren, I., Computer Arithmetic Algorithms, Prentice Hall, Englewood 
Cliffs, NJ (1993); 
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33. Fu, D., and Willson, Jr., A.N., "Interpolation in timing recovery using a 
trigonometric polynomial and its implementation," in Globecom 1998 
Comm, TheoryMmiConf Record (Nov. 1998), pp. 173-178; 

34. Lang, T., and Antelo, E., IEEE Trans, Comp, ^7:736-749 (1998); 

35. Moeneclaey, M., IEEE Trans, Comm. 37:1029-1032 (1983); 

36. Vuori, J., IEEE Trans. Consum, Electron. ^5:1098-1 103 (1999); and 

37. Freeman, H., Discrete-Time Systems, Wiley, New York, NY (1965). 

10. Conclusion 

Example implementations of the methods, systems and components of the 
invention have been described herein. As noted elsewhere, these example 
implementations have been described for illustrative purposes only, and are not 
limiting. Other implementation embodiments are possible and covered by the 
invention, such as but not limited to software and software/hardware 
implementations of the systems and components of the invention. Such 
implementation embodiments will be apparent to persons skilled in the relevant 
art(s) based on the teachings contained herein. 

While various application embodiments of the present invention have been 
described above, it should be understood that they have been presented by way of 
example only, and not limitation. Thus, the breadth and scope of the present 
invention should not be limited by any of the above-described exemplary 
embodiments. 
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